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L  SUMMARY 


A.  Introduction 

This  Annual  Report  details  a  'v-ogram  of  research  on  composite  spinel 
ceramic  materials  with  enormous  specific  heat  maxima  in  the  temperature  range 
5-20  K  carried  out  by  CeramPhysics,  Inc.,  Ohio  State  University,  Pennsylvania 
State  University,  and  West  Virginia  University.  In  this  section  the  main 
advances  and  accomplishments  of  the  year  are  summarized  and  the  organization 
of  the  rest  of  the  report  outlined. 

We  start  with  a  brief  review  of  the  interpretation  of  these  spinel  materi¬ 
als  developed  in  our  previous  research.  The  results  obtained  so  far  have 
suggested  an  attractive  picture  for  the  ordering  phenomena  in  the  B-site 
spinels,  CdCr204  and  2nCr20^.  The  results  have  shown  that  at  least  two  types 
of  magnetic  correlations  are  present,  antiferromagnetic  and  paraaiag"®^ ic.  and 
that  frustration  and  the  presence  of  strong  spin-lattice  coupling  play  an 
important  role  in  the  anomalously  large  specific  heats  and  thermal  conductivi¬ 
ties. 

A  basic  examination^  of  the  structure  of  the  spinel  phases  of  CdCr204  and 
ZnCr204  revealed  an  interesting  pattern  among  the  B-slte  spinels,  suggesting 
that  additional  systems  of  considerable  interest  might  be  made  by  filling  in  a 
table  of  materials  constructed  by  replacing  the  A-site  atom  (Zn  or  Cd)  by 
isoelectronic  atoms  or  mixtures  of  atoms.  In  addition  th"'  transitions  in 
these  materials  were  seen  to  have  a  peculiar  nature  in  which  the  spins  order 
weakly  in  a  lattice  which  has  a  high  degree  of  frustration.  This  has  great 
importance  for  our  understanding  of  these  systems,  since  it  means  that  large 
numbers  of  spins  can  remain  unordered  below  the  transition,  resulting  in 
anomalously  lar'-  ■'pecific  heats,  and  furthermore,  that  distortions  of  the 


_ ..iiarstri 

pri .  J  IK  ^ 

I 


Rr/iS^#Ai 


EPORT  DOCUMENTATION  PAGE 


ta.  REPORT  SECURITY  CLASSIFICATION 
Unclassified 


2a  SECURITY  CLASSIFICATION  AUTHORITY 


b  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


A  PERFORMING  ORGANIZATION  REPORT  NUMBER{S) 


6b  OFFICE  SYMBOL 
(If  applicable) 


6a.  NAME  OF  PERFORMING  ?  ■  C aNIZATION 
CeramPhysics ,  Inc. 


6c  ADDRESS  (Gty,  State,  and  ZIP  Code) 

921  Ea&twind  Drive,  Suite  110 
Westerville,  Ohio  43031 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

AFOSR 


8c.  ADDRESS  (City,  State,  and  ZIP  Code) 

Bldg.  410 

Bolling  AFB  DC  20332-6448 


1 1  title  (Include  Security  Classification) 


Research  on  High  Specific  Heat  Dielectrics 


lb  RESTRICTIVE  MARKINGS 
None 


3  DiSTRlBUTiON/AVAILASlIITY  Of.  REPORT 

lppi*owd  tov  puWlfl 

d  i  St  T  Ibut  1  on  ual  l«it  ed 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

IlFOSR-TR.  8  7  -  145  2 


7a  NAME  OF  MONITORING  ORGANIZATION 

Air  Force  Office  of  Scientific  Research 


7b  A.NDRESS  (C/ty,  Stare,  and  ZIP  Code) 
Bldg.  410 

Bolling  AFB  DC  20332-6448 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F49620-86-C-0049 


10  SOURCE  OF  FUNDING  NUMBERS 


PROJECT 

NO. 


WORK  UNIT 
ACCESSION  NO 


12  PERSONAL  AUTHOR(S)  Lawless ,  W.N.  and  Clark,  C.F.,  C^eramPhysics ,  Inc.);  Patton,  B.R.,  (Ohio 

d  Dalai.  N . S . . (W. Virainia  Univ . ) i Snrout .  T.  (Pennsylvania  State  Uni 


14.  DATE  OF  REPORT  (Year,  Month.  Day)  15  PAGE  COUNT 

1987,  May  10  188 


16  SUPPLEMENTARY  NOTATION  01/30 


WgiimuaoigaKigEiiCT 


17 

COSATI  COOES  1 

FIELD 

1  GROUP 

SUB-GROUP  1 

mtmmam 

18  SUBJECT  terms  {Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Sl2|  i\  ,  ;  c  , 


19  ABSTRm 

T 

c  e  r  am i 
range 
p  i  c  t  u  r 
Z nC  r ^O 

of  the 
well  a 
M 

that  t 
and  Cd 
z  a  t  i  o  n 
p  a  r  a  m  a 


(Cont,r,i>€ 
his  An 
c  mate 
5-2  0  K 
e  f  o  r 

4  ‘  ..  . 

1  g  n  1  f  1 

s  p  i  n  e 

s  the 

a  g  n  e  t  o 

he  p  ho 

spine 

c  o  o  1  i 

g  n  e  c  i  0 


on  reverse  if  necessary  and  identify  by  block  number) 

nual  Report  details  a  program  of  research  on  composite 
rials  with  enormous  specific  heat  maxima  in  the  temper 
.  The  results  obtained  so  far  have  suggested  an  attra 
the  ordering  phenomena  in  the  B-site  spinels,  CdCrl^O, 

cant  progress  in  characterizing  and  controlling  the  fa 
1  materials  has  been  made.  Pure  powders,  compacted  di 
9/1  columbite  composite  have  been  made  and  characteriz 
caloric  experiments  reveal  no  trace  of  hysteresis,  whi 
noraena  involved  in  the  low  temperature  transitions  in 
Is  must  involve  second  order  transitions.  Adiabatic  d 
ng  above  T’^  and  also  at  low  temperatures  is  indicative 
spins,  while  a  region  of  demagnetization  heating  .just 


spinel 
a  t  u  r  e 
c  t  i  ve 
and 

b  r  i  c  a  t  i  o  n 
s  k  s  ,  as 
ed  . 

c  h  means 
the  Z  n 
e  m  ;i  g  n  e  t  i  - 
of 

below  r . . 

.'l 


20  OlSTRiBuTlON/AVAlLABiLiTv  OF  A^TTRACT 
OUNCLASSIFIED/UNLIMITED  ©5AME  AS  RPT 


DO  FORM  1473,  34  mar 


3  3  A-iR  eUit  on  m<iy  be  used  urni  i?»r.nistcd 
All  Otbpr  editions  J't.‘  ntisoletp 


_ srcj  miTs'^Cl  ASSiMCATiON  Of  Jjs  PA'Sf 

ilNCLASSIHED 


19.  ABSTRACT  (continued) 


seems  to,  indicate  ant  if e c romagne 1 1  c  ordering.  Measurements  have  also 
discovered  an  amazing  anomaly  in  the  dielectric  constant  at  the  anti¬ 
ferromagnetic  transition.  >. 

The  susceptibility  mea^sur  ament  s  on  the  pure  spinel  powder  reveal  no 
low  temperature  paramagnetic Curie  tail.  Implying  that  all  the  spins  are 
well  ordered  by  T  =  0.2  T^.  On  the  other  hand  the  grains  reacted  with  10% 
columbite  reveal  a  low  temperature  paramagnetic  tail.  Comparison  of  the 
measured  s.u  s  c  e  [i  t^i.|)  i  1  it  y  with  a  mean-field  Currie-Weiss  law  above  T^,  and  a 
Curie  law  at  1  o'w’ tdmphratur e  seems  to  imply  that,  at  low  -temperatures  at 
least,  the  number  of  paramagnetic  spins  N  is  a  small  fraction  of  the 
anti -ferromagnetic  spins  of  order  l%^in  the  Cd  spir--l  and  0.1%  in  the 

Zn  spinel. 

Theoretical  Monte  Carlo  calculations  have  found  that  the  spinel 
lattice  with  spins  only  at  the  B  sites  has  a  magnetic  transition  at  a 
relatively  low  temperature.  The  magnetization  is  linear  in  field  in  good 
agreement  with  experiment.  The  calculated  low  field  susceptibility  shows  a 
modest  cusp  at  an  ordering  temperature  T^  =  J,  where  J  is  the  nearest 
neighbor  an t i -f e r r omag ne t i c  exchange  interaction. 


lattice,  which  remove  the  frustration,  can  couple  strongly  to  the  spins,  thus 
leading  to  dielectric  anomalies  and  large  thermal  conductivities  due  to  spin 
energy  being  transported  through  the  spin-phonon  interaction. 

R  Current  Work 

Recent  and  interesting  new  work  has  been  done  under  this  proposal  on  the 
properties  of  the  B-site  spinels,  CdCo204  and  ZnCo204.  In  this  summary,  we  we 
focus  on  1)  recent  advances  in  the  fabrication  and  structural  characterisation 
of  various  ceramic  phases  of  CdC0204  and  ZnC02O4  (section  III),  2) 
experimental  research  on  the  magnetocaloric  properties  of  the  B-site  spinels 
and  the  experimental  discovery  of  a  significant  dielectric  anomaly  at  T]^, 
which  reveals  a  novel  coupling  of  the  spin  and  lattice  (section  IV),  3) 
experimental  work  on  the  magnetization  and  susceptibility  which  reveal  an 
antiferromagnetic  transition  at  Tjj  «  6-llK  (section  V),  4)  theoretical  results 
for  a)  the  properties  of  the  ideal  spinel  lattice  and  b)  the  properties  of  the 
tetragonally  distorted  lattice  obtained  from  Monte  Carlo  computer  simulations 
on  the  frustrated  spinel  spin  lattice  (section  VI),  and  5)  an  overall  picture 
for  the  phenomena  occurring  in  these  fascinating  materials  which  integrates 
the  numerous  results  listed  above  (section  VII).  We  describe  each  of  these 
accomplishments  in  more  detail  in  the  following: 

1)  Significant  progress  in  characterizing  and  controlling  the  fabrication 
of  the  spinel  materials  CdCo204  and  ZnCo204  has  been  made.  Pure  powders, 
compacted  disks,  as  well  as  the  9/1  columbite  composite  have  been  made  and 
characterized.  The  fabrication  experiments  in  section  III  give  clear  evidence 
that  the  spinel  powder  of  nominal  grain  size  10  microns  reacts  with  the 
columbite  component  to  produce  paramagnetic  spins  at  low  temperatures.  The 
same  ceramic  material  can  be  made  in  several  different  ways,  with  about  lox 


addition  of  the  columbite  phase.  SEM  evidence  indicates  that  the  spinel  and 
columbite  phases  remain  distinct,  forming  a  composite  of  the  two  different 
substances,  with  grain  sizes  of  order  lOp,  consistent  with  earlier  thermody¬ 
namic  determinations  at  CeramPhysics. 

2)  The  magnetocaloric  experiments  reported  in  section  IV  reveal  no  trace 
of  hysteresis,  which  means  that  the  phenomena  involved  in  the  low  temperature 
transitions  in  the  Zn  and  Cd  spinels  must  involve  second  order  transitions. 
Adiabatic  demagnetization  cooling  above  and  also  at  low  temperatures  is 
indicative  of  paramagnetic  spins,  while  a  region  of  demagnetization  heating 
just  below  Tjj  seems  to  indicate  antiferromagnetic  ordering.  These  results  are 
consistent  with  the  thermodynamic  relation  which  requires  that  the  sign  of  the 
temperature  change  (heating  or  cooling)  upon  adiabatic  demagnetization  depends 
on  the  sign  of  the  temperature  derivative  of  the  susceptibility.  Measurements 
have  also  discovered  an  amazing  anomaly  in  the  dielectric  constant  at  the 
anti-ferromagnetic  transition.  This  data  clearly  shows  that  the  ordering  of 
the  frustrated  spin  system  has  an  enormous  effect  on  the  structure  of  the 
lattice.  We  recall  that  the  pattern  in  the  B-slte  spinel  single  crystals  is 
toward  a  slight  tetragonal  distortion  of  the  high  temperature  cubic  phase 
below  the  magnetic  transition.^  We  know  of  no  other  case  of  such  a  strong 
coupling  between  the  magnetic  and  the  lattice  degrees  of  freedom  in  any 
material.  The  present  measurements  are  the  first  to  find  an  anomaly  at  low 
temperatures.  A  theory  was  developed  (section  VI)  which  gives  both  a  qualita¬ 
tive  explanation  as  well  as  quantitative  agreement  with  this  anomaly. 

3)  The  susceptibility  measurements  described  in  section  V  on  the  pure 
spinel  powder  reveal  no  low  temperature  paramagnetic  Curie  tail,  Implying  that 
all  the  spins  are  well  ordered  by  T  s  0.2  Tjj.  On  the  other  hand  the  grains 
reacted  with  lOZ  columbite  reveal  a  low  temperature  paramagnetic  tail.  Compar¬ 
ison  of  the  measured  susceptibility  with  a  mean-field  Curie-Welss  law  above 


Tjj,  and  a  Curie  law  at  low  temperature  seems  to  imply  that,  at  low  tempera¬ 
tures  at  least,  the  number  of  paramagnetic  spins  Np  is  a  small  fraction  of  the 


anti-ferromagnetic  spins  Naf,  of  order  ly-  in  the  Cd  spinel  and  O.ix  in  the  Zn 
spinel. 

4)  Finally,  the  theoretical  Monte  Carlo  calculations  reported  in  section 


VI  have  found  that  the  spinel  lattice  with  spins  only  at  the  B  sites  has  a 
magnetic  transition  at  a  relatively  low  temperature.  The  magnetization 
is  linear  in  field  in  good  agreement  with  experiment.  The  calculated  low 


field  susceptibility  shows  a  modest  cusp  at  an  ordering  temperature  Tjj  ~  J, 
where  J  is  the  nearest  neighbor  anti-ferromagnetic  exchange  interaction. 

In  conclusion,  it  is  clear  that  a  wide  variety  of  experimental  and  theoretical 
work  on  the  frustrated  B-site  spinels,  CdCr204  and  ZnCr204  has  begun  to 
produce  a  coherent  picture  for  the  very  interesting  behavior  of  these  systems. 
K  magnetic  transition  of  anti-ferromagnetic  character,  Incomplete  ordering  of 
the  spins,  and  strong  coupling  to  the  lattice  degrees  of  freedom  are  all  well 
established.  Further  work  is  required  to  pin  down  the  origins  and  parameters 
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II. 


INTRODUCTION  AND  BACKGROUND 


In  the  late  70's,  Wright  Patterson  AFB  began  funding  a 
series  of  applied  programs  aimed  at  developing  dielectric 
insulation  systems  incorporating  enthalpy  stabilization  for  the 
superconductors  NbTi  and  Nb^Sn  (WPAFB  Contracts  #F33615-80-C- 
2202,  -82-C-2227,  -82-C-2229,  and  -84-C-2418).  These  programs 
were  based  on  new  materials  which  are  two-phase,  spinel  + 
columbite  ceramics  and  which  sinter  at  ~  1300“C.  The  columbite 
phase,  although  minor,  must  be  present  for  the  spinel  phase  to 
density,  and  these  ceramics  have  huge  specific  heats  at  low 
temperatures . 

In  the  Summer  of  '83,  research  on  the  physics  of  these  new 
materials  was  begun  by  CeramPhysics ,  Inc.  and  the  Physics  Dept, 
at  the  Ohio  State  Univ.  under  a  subcontract  to  AFOSR  Contract 
#F49620-82-C-0129 .  This  research  is  summarized  in  Eckels  e^  al . 
(1985),  and  this  Annual  Report  documents  the  continuation  of  this 
research  by  CeramPhysics,  Inc.,  Ohio  State  Univ.,  Pennsylvania 
State  Univ.,  and  W.  Virginia  Univ. 

The  columbite  mineralizers  CdNb20g  and  ZnNb20g  must  be 
present  for  the  spinels  CdCr20^  and  ZnCr20^  to  density  in  the 
ceramic,  and  the  optimum  spinel rcolumbite  molar  ratio  is  9:1. 

The  thrust  of  the  previous  research  has  been  to  understand  the 
enormous  specific  heat  maxima  in  these  spinels  at  low  tempera¬ 
tures  rather  than  the  ceramic-formation  phenomena. 

Specific  heat  data,  1.5-40  K,  were  measured  on  these  9:1 
ceramics  and  on  ceramic  samples  of  the  CdNb20g  and  ZnNb20g 
mineralizers  over  broad  temperature  ranges.  Next,  the  specific 
heats  of  the  columbites  were  fitted  to  Schottky  and  Einstein 
functionals,  and  these  fitted  data  were  then  used  to  correct  the 
9:1  ceramics  specific  heat  data  to  determine  the  specific  heats 
of  the  spinel  phases. 

The  separated  specific  heat  data  for  the  CdCr20^  and  ZnCr20^ 
spinels  are  shown  in  Fig.  2-l(a)  where  the  data  are  plotted  as 
C/T^  to  illustrate  the  non-Debye  contributions  at  the  lowest 
temperatures.  We  point  out  that  these  specific  heat  maxima  are 
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C/T  (ergo 


equivalent  to  the  specific  heat  of  water  at  room  temperature; 
hence,  the  applied  interest  in  enthalpy  stabilization. 

These  specific  heat  data  were  decomposed  into  Debye, 
Schottky,  and  magnetic  contributions,  and  the  entropies  of  these 
contributions  in  CdCr20^  are  shown  in  Fig.  2-l(b).  These  contri- 


T  (K) 


(a) 


(b) 


Figure  2-1.  (a)  Specific  heats  of  the  new  spinels;  (b) 

Contributions  to  the  specific  heat  of  CdCr20^. 


butions  were  determined  by  the  following  procedure  based  on  the 
data  in  Fig.  2-l(a): 

1.  The  high-temperature  data  (i.e.,  T>2Tj^,  where  T^^  is  the 
peak  temperature  in  Fig.  2-l(a)]  were  fitted  to  the 
general  Schottky  term,  and  these  fittings  yielded  the 
Debye  temperatures,  =  420  and  463  K  for  the  cadmium 
and  zinc  spinels,  respectively.  These  0^^  -  values  are  in 
excellent  agreement  with  predictions  of  the  Lindemann 
relation . 

2.  The  low-temperature  data  (i.e.,  T<Tj^/2)  follow  a  Schottky 
term  very  well  with  two-level  splittings  of  0.33  and  0.37 
K  for  the  cadmium  and  zinc  spinels.  However,  these  fits 
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K  (rrf^cm  K  ) 


yielded  unrealistically  small  0jj's  (e.g.,  ~  100  K),  and 
it  was  concluded  that  antiferrimaqnetic  spin  waves  are 
contributing  a  term  which  is  indistinguishable  from 
the  Debye  acoustic  background. 

3.  Adopting  the  more  realistic  Debye  temperatures  allows  a 
determination  of  the  spin-wave  parameters, 

J's/c^^/^  =  5.94x10"^®  erg  (CdCr204) 

=  9.65x10  erg  (ZnCr20^) 

4.  Finally,  the  determination  of  the  parametric  data  above 
allows  the  separation  of  the  specific  heat  contributions 
and,  by  integration,  the  entropy  of  these  contributions, 
as  illustrated  in  Fig.  2-l(b).  The  magnetic  entropies  of 
the  two  spinels  determined  from  these  analyses  are 

=  2.340  (CdCr20^) 

=  1.434  (ZnCr204) 


(a)  (b) 

Figure  2-2.  Thermal  Conductivity  data  on  the  spinels,  (a) 
As  a  function  of  temperature;  (b)  In  magnetic  fields. 

Here  C(9/l)  and  D(9/l)  refer  to  the  Cd-and  Zn-spinels, 
respectively. 


The  thermal  conductivities  of  the  spinels  display  "jumps" 
associated  with  the  (antiferrimagnetic)  transitions  as  shown  in 
Fig.  2-2(a).  The  magnetic-field  dependence  of  these  anomalies 
was  measured  just  below  T^^,  and  these  data  are  shown  in  Fig. 
2-2(b).  Surprisingly,  the  anomaly  in  ZnCr20^  is  unaffected,  but 
in  CdCr20^  a  magnetic  field  "quenches"  the  anomaly  as  shown 
schematically  by  the  dashed  curve  in  Fig.  2-2(a).  This  excess 
thermal  conductivity  in  CdCr20^  is  =  0.23  mW  cm  ^  K  ^ ,  or  roughly 
half  the  zero-field  thermal  conductivity  at  5  K. 

By  far  the  most  intriguing  measurements  previously  made  on 
the  spinels  have  been  the  maqnetocaloric  measurements  at  T<Tj^/2. 
Both  spinels  display  adiabatic-demagnetization  cooling  and 
adiabatic-magnetization  heating,  and  these  caloric  effects  are 
perfectly  reversible .  The  measured  AT  data  up  to  10  T  are  shown 
in  Fig.  2-3,  and  we  point  out  that  these  are  very  large  effects 
because  the  specific  heats  of  these  spinels  are  large  at  3-5  K. 


Figure  2-3.  Reversible  magnetocaloric  effects  in  the  new 
spinels  at  T<T^/2.  Here  CCN(9/1)  and  ZCN(9/1)  refer  to  the 
9:1  spinel  +  columbite  ceramics. 

The  measurements  in  Fig.  2-3  led  naturally  to  measurements 
of  the  magnetization  of  the  spinels  at  4.2  K,  and  these  data  are 
shown  in  Fig.  2-4.  Here  it  was  found  that  both  spinels  act  as 
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perfect  paramaqnets  and  display  no  hysteretic  effects. 

A  self-consistent  picture  began  to  emerge  from  these 
measgrements :  Namely,  the  specific  heat  peaks  in  Fig.  2-l(a) 

appear  due  to  an  antiferrimagnetic  spin  ordering  as  evidenced  by 
the  antiferrimagnetic  spin-wave  contribution  below  about  3  K. 
These  ordered  spins  do  not  affect 
the  magnetization.  Fig.  2-4.  How¬ 
ever,  several  degrees  of  freedom 
remain  below  and  these  para¬ 

magnetic  spins  cause  the  reversi¬ 
ble  magnetocaloric  effects. 

Actually,  there  is  competition 
between  the  ordered  spins  and  the 
disordered  spins  in  the  magneto¬ 
caloric  effects  because  the  former 
spin  system  would  show  demagnet¬ 
ization  heating,  the  latter  system, 
demagnetization  cooling.  There  may 
be  different  spin-phonon  relaxation 
rates  involved  which  cause  the  do¬ 
minance  of  the  latter  effects  assoc¬ 
iated  with  the  paramagnetic  spin 
system . 

It  is  puzzling  that  both  the  specific  heat  maxima  at  T^^  and 
the  magnetocaloric  effects  below  Tj^  are  so  large  (i.e.,  the 
energetics  associated  with  these  phenomena  are  very  large,  ~ 

0.1-1  J/g) . 

A  very  simple  approach  was  tried  theoretically  here;  namely, 
a  two-spin  model  where  some  fraction  of  the  Cr^^  spins  partake  in 
the  antiferrimagnetic  ordering  at  T^^,  the  balance  of  the  spins 
remain  unordered,  and  coupling  between  the  spin  systems  is 
ignored.  Using  s  =  3/2,  it  was  found  from  the  magnetic  entropy 
Sjn  above  that  84  and  52%  of  the  spins  order  at  Tjg  in  CdCr204 
ZnCr20^,  respectively.  From  the  spin-wave  specific  heat  contri¬ 
butions  above,  we  have  for  the  exchange  constants. 


H  (T) 

Figure  2-4.  Magnetiza¬ 
tion  data  at  4.2  K  on 
the  new  spinels. 


J/k  =  3.42,  CdCr20^, 

=  4.72,  ZnCV20/^‘ 

From  the  Schottky  specific  heat  terms  well  below  we  find  for 
the  two-level  splitting, 

6  =  0.59  K,  CdCr20^, 

=  0.37  K,  ZnCr20^, 

where  these  terms  are  ascribed  to  the  unordered  spins. 

This  simple  model  yields  reasonable  parameter  values  from 
the  specific  heat  data,  but  the  model  cannot  explain  the  magnetic 
properties.  From  the  magnetization  data.  Fig.  2-4,  the  model 
indicates  that  the  spin  densities  in  CdCr20^  are  three  times 
larger  than  in  ZnCr20^,  but  the  magnetocaloric  data  indicate  just 
the  opposite  from  this  simple  model. 

The  more  basic  theoretical  approach  started  with  a  basic 
examination  of  the  structure  of  the  spinel  phases  of  CdCr20^. 

This  study  revealed  an  interesting  pattern  among  the  B-site 
spinels,  which  suggests  that  additional  systems  of  considerable 
interest  might  be  made  by  filling  in  a  table  of  materials  con¬ 
structed  by  replacing  the  A-site  atom  (Zn  or  Cd)  by  isoelectronic 
mixtures  of  Cu  and  In,  or  Ag  and  In.  In  addition,  the  transi¬ 
tions  in  these  materials  were  seen  to  have  a  peculiar  nature  in 
which  the  spins  order  weakly  in  a  lattice  which  has  a  high  degree 
of  frustration .  This  has  great  importance  for  our  understanding 
of  these  systems,  since  it  means  that  large  numbers  of  spins  can 
remain  unordered  below  the  transition,  resulting  in  anomolously 
large  specific  heats  and  magnetocaloric  effects.  Furthermore, 
distortions  of  the  lattice,  which  remove  the  frustration,  can 
couple  strongly  to  the  spins,  thus  leading  to  large  thermal 
conductivities  due  to  spin  energy  being  transported  through  the 
spin-phonon  interaction. 

The  first  attempt  to  understand  the  theoretical  properties 
of  these  systems  was  to  derive  the  Hamiltonian  describing  the 
interaction  of  the  spins  on  the  spinel  lattice  sites.  From  this 
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we  extracted  a  simple  model  in  which  we  kept  only  two  possible 
sublattice  magnetizations,  an  antiferromagnetic  one  which  orders 
at  the  transition,  and  a  ferromagnetic  one  coupled  to  it,  which 
is  polarized  in  a  magnetic  field,  giving  rise  to  spin-flop  like 
effects.  In  the  mean-field  approximation,  this  model  already 
gives  interesting  effects,  including  a  temperature  dependent 
magnetization  linear  in  field  and  a  magnetic  field  dependent 
specific  heat  which  increases  with  field.  These  results  are  in 
good  qualitative  agreement  with  the  experimental  data  and  confirm 
the  basic  correctness  of  our  approach.  The  next  step  was  to 
examine  the  renormalization  effects  due  to  fluctuations,  which 
led  to  the  novel  feature  that  the  non-ordering  spins  give  rise  to 
a  large  contribution  to  the  specific  heat  even  well  below  the 
ordering  temperature  of  the  antiferromagnetic  spins.  This  is  not 
possible  in  a  simpler  system  with  only  one  order  parameter,  again 


(a)  and  thermal  conductivity  (b)  data  for  ZnCr20^. 

suggesting  our  model  has  correct  general  properties.  Detailed 
results  were  obtained  in  the  gaussian  fluctuation  approximation 
for  the  specific  heat  as  a  function  of  temperature  and  magnetic 
field  both  above  and  below  the  transition.  Finally,  a  theory  for 
the  thermal  conductivity  in  this  class  of  materials  was  de¬ 
veloped.  The  thermal  conductivity  depends  on  both  the  specific 


-11- 


‘aJ- ‘Ata v«»-' 


heat  and  the  scattering  rate,  which  we  showed  involved,  in  addi¬ 
tion  to  the  usual  phonon  and  antiferromagnetic  spin  wave  contri¬ 
butions,  a  spin  fluctuation  term  which  limits  the  thermal  con¬ 
ductivity  in  the  vicinity  of  the  transition.  Fits  were  made  to 
the  specific  heat  and  thermal  conductivity  data  using  the  same 
set  of  fitting  parameters,  which  revealed  that  the  peak  in  the 
thermal  conductivity  at  the  transition  arises  from  the  peak  in 
the  specific  heat,  but  that  it  is  typically  a  factor  of  about 
three  smaller  due  to  the  spin  scattering  at  the  transition.  An 
example  of  the  theoretical  fits  for  the  ZnCr20^  spinel  are  shown 
in  Fig.  2-5. 

The  research  reported  in  this  Annual  Report  was  based  on  the 
above  experimental  and  theoretical  findings  and  represents  a 
continuation  of  this  previous  research. 

Our  broad  goals  in  the  present  research  are  as  follows: 

1.  To  investigate  the  ceramic  formation  of  these  new  spinel 
+  columbite  materials  (i.e..  Why  must  the  columbite 
phase  be  present  for  densif ication? )  and  to  explore 
ceramic  fabrication  of  the  pure  spinels  (Pennsylvania 
State  Univ . ) . 

2.  To  extend  the  magnetocaloric  measurements  of  Pig.  2-3 
over  broad  temperature  ranges  for  both  the  spinel  + 
columbite  ceramics  and  for  the  pure  spinel  ceramics 

( CeramPhysics ) . 

3.  To  extend  the  specific  heat  measurements  of  Fig.  2-1  to 
the  pure  spinel  ceramics  (CeramPhysics). 

4.  To  measure  magnetic  susceptibility  data  on  these 
materials  over  broad  temperature  ranges  and  in  magnetic 
fields  (W.  Va.  Univ.). 

5.  To  pursue  EPR  measurements,  including  electron  nuclear 
triple  resonance  measurements,  on  these  spinel  systems 
( W.  Va .  Univ . )  . 

6.  To  pursue  Monte  Carlo  simulations  of  the  spinel  lattice, 
including  frustratin  and  variable  boundary  conditions, 
and  to  couple  these  calculations  to  Ginzberg-Landau 
theory.  Renormalization  group  theory  will  be  employed 
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III.  CERAMIC  PREPARATION  STUDIES 


3. A  Introduction 


A  key  to  the  understanding  of  the  thermal,  dielectric,  magnetocaloric, 
and  magnetic  susceptibility  measurements  performed  In  this  program  Is  the 
crystal  chemistry  and  ceramic  fabrication  of  the  CdCr204  and  ZnCr204  spinel 
materials.  A  portion  of  this  program  was  devoted  to  this  topic.  In  earlier 
CeramPhysIcs  studies,  It  was  found  that  dense  ceramics  of  the  CdCr204  or  ZnCr204 
spinels  could  only  be  fabricated  by  the  addition  of  at  least  10  mole  percent 
of  the  respective  (CdNb20g  or  ZnNb20g)  columbites.  Ceramic  compositions  with 
various  spinel /col umbite  ratios  were  the  basis  of  earlier  studies.  Several 
of  these  existing  ceramic  samples  were  thoroughly  characterized  by  x-ray 
diffraction  and  scanning  electron  microscopy.  In  order  to  better  understand 
the  crystal  chemistry  of  these  materials.  It  was  hoped  that  If  an  understanding 
of  the  crystal  chemistry  and  densificatlon  mechanisms  of  the  spinel -col umbite 
ceramics  could  be  obtained,  fabrication  methods  could  be  established  to  allow 
for  the  preparation  of  more  phase-pure  spinel  ceramics,  with  Improved  thermal 
properties.  Ceramic  preparation  studies  were  also  performed  In  order  to  prepare 
dense  samples  of  the  pure  spinel  ceramics.  These  efforts  were  aimed  at 
supplying  Improved  ceramic  samples  for  thermal  and  magnetic  measurements. 

3.B  Spinel -Col umbite  Ceramics  (CeramPhysIcs  Samples 


3.B.1  X-Ray  Diffraction  Studies 


X-ray  diffraction  and  SEM  studies  were  performed  on  previously 
synthesized  spinel -col umbite  ceramics  (both  calcined  powders  and  sintered 
disks)  with  spinel /col umbite  ratios  varying  from  1/1  to  11/1.  X-ray  diffraction 
results  are  presented  In  Table  3.1;  the  phases  found  for  each  of  the  samples 
are  listed  In  order  of  decreasing  relative  amount.  The  results  for  the 
CdCr204/CdNb20g  samples  (Table  3.1a)  confirmed  the  multi -phase  nature  of  these 
ceramics.  With  low  ratios  {<5/1)  of  spinel  to  columbite,  the  major  phase  In 
the  ceramic  was  pyrochlore  Cd2Nb207,  with  lesser  amounts  of  spinel  CdCr204, 
columbite  CdNb20g,  and  Nb205.  As  the  spinel /col umbite  ratio  Increased,  the 
amount  of  pyrochlore  decreased,  and  the  amounts  of  spinel  and  columbite 


Increased.  Spinel  CdCr204  was  the  dominant  phase  with  spinel /col umbite  ratios 
of  5/1  and  greater.  With  high  spinel/columbite  ratios  (9/1),  the  pyrochlore 
phase  was  not  observed. 

X-ray  diffraction  patterns  of  the  ZnCr204/ZnNb205  samples  (Table 
3.1b)  were  much  simpler.  Spinel  ZnCr204  was  the  dominant  phase  for  all  of 
the  ceramics  (ratios  of  3/1  to  9/1).  Columbite  ZnNb20g  phase  was  the  only 
other  phase  In  these  samples  (except  for  a  trace  of  excess  Nb205  detected  in 
the  5/1  sample.  The  amount  of  columbite  decreased  as  the  spinel/columbite 
ratio  was  Increased,  as  expected. 

Additional  ceramic  samples  were  studied  by  x-ray  diffraction,  and 
gave  Interesting  results.  The  compositions  were  Cd5Zn7Cr2oNb4052  and 
Cd7Zn5Cr2oNb4052;  the  Intent  here  was  to  form  a  50/50  solid  solution  of  CdCr204 
and  ZnCr204  spinels,  witha  10  percent  mineralizer  of  the  ZnNb20g  and  CdNb20g 
columbites,  respectively.  As  shown  by  the  data  In  Table  3.1c,  the  spinel  solid 
solution  was  not  formed.  For  each  composition,  the  same  three  phases  were 
formed;  spinel  ZnCr204,  columbite  CdNb205,  and  CrNb04.  This  suggests  that 
the  ceramic  solid  solution  of  these  two  spinels  will  be  difficult  (If  not 
Impossible)  to  fabricate. 

3.B.2  Microstructure 

X-ray  diffraction  of  the  CeramPhysIcs  samples  revealed  that  most  of 
the  ceramics  consisted  primarily  of  two  phases,  the  respective  Cd  or  Zn  spinels 
and  columbites.  SEM  backscattered  and  x-ray  fluorescence  analyses  of  the 
sintered  ceramic  sample  of  the  5/1  CdCr204/CdNb205  composition  are  presented 
In  Figure  3.1.  The  SEM  (backscatter  mode)  micrograph  (Figure  3.1a)  reveals 
the  coexistence  of  two  phases,  with  the  major  phase  being  light  In  contrast 
and  a  minor  phase  evident  as  round,  greyish  spots,  which  appear  to  be 
precipitates.  X-ray  fluorescence  analyses  (Figures  3.1b,  c,  and  d)  Indicate 
that  the  major  phase  Is  rich  In  Cd  and  Cr,  and  Is  probably  the  CdCro4  spinel 
phase.  The  grey  areas  are  rich  In  Nb  and  deficient  In  Cr,  suggesting  that 
this  phase  Is  CdNb20g. 

The  microstructure  of  the  ZnCr204/ZnNb20g  samples  were  studied  by 
both  optical  and  scanning  electron  microscopy.  Optical  microstructures  of 
polished  faces  of  three  sintered  ceramic  disks  with  spinel/columbite  ratios 
of  3/1,  5/1,  and  7/1  are  presented  In  Figure  3.2.  Two  types  of  microstructures 


are  clearly  indicated,  with  a  core-like  structure  at  the  center  of  the  disks. 

It  Is  Important  to  note  that  the  core  size  Increases  with  decreasing 
spinel /col umbite  ratio.  Higher  magnifications  of  the  7/1  sample  were  studied 
by  SEM,  and  are  shown  In  Figure  3.3.  The  SEM  (backscatter  mode)  micrograph 
(Figure  3.3a)  revealed  that  both  the  outer  ring  and  inner  core  are  composed 
of  two  phases,  evident  as  grey  (low  contrast)  and  light  (bright  contrast) 
areas.  The  outer  ring  Is  composed  primarily  of  the  grey  areas  with  Isolated 
regions  rich  In  the  light  phase.  The  Inner  core  Is  also  composed  of  a  mixture 
of  the  of  two  phases,  but  with  a  finer  distribution  than  the  outer  ring.  X- 
ray  fluorescence  revealed  that  the  samples  are  uniform  In  Zn  (Figure  3.3b), 
the  light  areas  are  rich  In  Nb  and  deficient  In  Cr  (Figure  3.3c),  and  the 
grey  areas  are  rich  In  Cr  and  deficient  In  Nb  (Figure  3.3d).  These  results 
suggest  that  the  Inner  core  Is  composed  of  a  fine  mixture  of  the  ZnCr204  spinel 
and  ZnNb205  col umbite  phases,  whereas  the  outer  ring  Is  rich  In  spinel  with 
Isolated  columbite-rich  regions. 

3.B.3  Initial  Conclusions 

The  XRD  and  SEM  characterization  studies  on  the  CeramPhysIcs 
spinel /col umbite  ceramic  samples  revealed  contrasting  phase  constitutions  and 
microstructures  between  the  Cd  and  Zn  spinel /col umbite  compositions.  Hoewever, 
the  results  these  studies  did  not  offer  any  conclusive  evidence  as  to  how  the 
columbite  affects  the  sintering  behavior  of  the  spinel  ceramics,  or  why  It  Is 
necessary  to  enhance  denslflcatlon.  Further  studies  were  thus  warranted. 

3.C  Ceramic  Studies  (Penn  State 


3.C.1  Pure  CdCroO^  and  ZnCroO^  Spinels 


In  order  to  understand  how  the  additions  of  the  columbites  (CdNb205 
and  ZnNb205)  were  necessary  to  density  the  respective  CdCr204  and  ZnCr204 
spinel  ceramics,  the  single  phase  spinel  ceramic  powders  were  fabricated  and 
sintering  studies  were  performed.  Reagent  grade  oxides  (CdO,  ZnO,  and  Cr203) 
were  batched  to  the  stoichiometric  spinel  compositions,  vibratory  milled  In 
nalgene  Jars  with  alcohol  and  zirconia  grinding  media,  and  dried.  Samples  of 
the  two  spinel  compositions  were  calcined  at  various  temperatures,  and  x-ray 
diffraction  was  used  to  monitor  the  formation  of  the  spinel  phase.  It  was 


found  that  the  CdCr20^  spinel  was  formed  by  calcination  at  950°C  for  four 
hours,  while  the  formation  of  the  ZnCr204  spinel  required  a  calcination 
temperature  of  1050®C. 

Disks  were  prepared  from  the  two  spinel  powders,  and  a  sintering 
study  was  performed.  The  disks  were  placed  on  zirconia  setters  and  sintered 
at  various  temperatures  for  one  hour.  The  results  of  this  sintering  study 
are  presented  In  Table  3.2.  It  was  found  that  neither  of  the  two  spinels 
could  be  densifled  to  near  their  theoretical  densities  of  5.79  and  5.30  g/cc, 
for  CdCr204  and  ZnCr204,  respectively.  The  maximum  densities  achieved  for  the 
the  spinel  samples  were  4.4  g/cc  (76  percent  of  theoretical)  for  the  CdCr204 
spinel  ceramic  sintered  at  1290®C,  and  3.8  g/cc  (72  percent  of  theoretical), 
for  the  ZnCr204  spinel  ceramic  sintered  at  1650®C.  It  was  Interesting  to  note 
that  below  1350®C,  the  fired  disks  experienced  a  weight  gain  of  between  1  and 
2  percent.  The  origin  of  this  weight  gain  Is  yet  unknown,  as  the  fired  disks 
remained  single  phase  spinel  by  XRD.  For  sintering  temperatures  above  1350®C, 
weight  loss  due  to  volatilization  of  CdO  or  ZnO  was  observed;  the  weight  loss 
was  very  pronounced  In  the  CdCr204  spinel  sample  sintered  at  1350®C  and  was 
accompanied  by  a  substantial  decrease  In  density.  Such  weight  loss  Is  typically 
expected,  as  both  ZnO  and  especially  CdO  are  volatile  at  these  high 
temperatures. 

The  inability  to  densify  the  two  spinel  ceramics  can  be  attributed 
to  the  volatilization  of  CdO  and  ZnO  at  the  high  sintering  temperatures  that 
are  required  for  these  spinel  ceramics.  If  the  volatilization  of  CdO  or  ZnO 
could  be  prevented,  then  densificatlon  of  the  ceramics  would  occur.  In  order 
to  prevent  weight  loss,  the  spinel  disks  were  sintered  In  a  CdO-  or  ZnO-rIch 
atmosphere,  by  burying  the  disks  In  a  sintered  powder  (or  sand)  of  their  own 
composition.  For  the  ZnCr204  samples,  this  technique  was  effective  In 
preventing  ZnO  loss,  and  dense  ceramics  (90  percent  of  theoretical)  were 
achieved,  although  extremely  high  sinterinc  temperatures  (1600®C)  were  required. 
The  SEM  microstructure  of  a  90  percent  dense  ZnCr204  sample,  sintered  at  1600®C 
for  four  hours.  Is  shown  In  Figure  3.4,  and  compared  with  the  SEM  microstructure 
of  the  analagous  Zn  spinel /col umbite  (9/1)  sample  sintered  at  1350°C.  The 
grain  structure  size  of  the  1600®C  sintered  ZnCr204  sample  was  non-uniform 
but  with  a  large  grain  size  (4  to  12  microns);  conversely,  the  Zn 
spinel/columbite  (9/1)  sample  had  a  much  more  uniform  microstructure  and  smaller 
grain  size  (about  1-2  microns). 


sintering  the  CdCr204  spinel  disks  using  the  "sand"  method  to  control 
CdO  loss  was  effective,  but  no  appreciable  difference  density  was  achieved. 
However,  with  a  sintering  temperature  of  1350°C,  the  weight  oss  was  reduced 
from  13  to  less  1  percent.  The  "sand"  technique  was  Ineffective  for  preventing 
CdO  loss  with  sintering  temperatures  above  1350®C,  and  resulted  In  lower 
densities.  Thus,  dense  samples  (90  percent  of  theoretical)  of  CdCr204  could 
not  be  prepared.  The  SEM  microstructure  of  the  CdCr204  sample  sintered  at 
1300°C  for  0.5  hours  In  a  CdO-rIch  atmosphere  Is  compared  with  that  of  the 
analagous  Cd  spinel /col umbite  (9/1)  sample.  In  Figure  3.5.  The  CdCr204  sample 
had  a  uniform  microstructure  with  a  grain  size  of  about  2  microns,  and  a 
significant  amount  of  porosity.  The  Cd  spinel /col umbite  (9/1)  sample  had  a 
non-uniform  microstructure  with  grain  sizes  ranging  from  4  to  8  microns,  but 
with  less  porosity. 


3.C.2  Modified  Spinel  Compositions 

Stoichiometry  and/or  other  compositional  variations  ofte.i  are 
Important  factors  In  determining  whether  a  particular  ceramic  powder  will 
density.  A  study  was  undertaken  to  determine  whether  compositional 
modifications  would  affect  the  sintering  behavior  of  the  spinel  ceramics.  It 
was  hoped  that  the  results  of  this  study  would  Improve  the  understanding  of 
how  the  columbite  additions  enhance  the  sintering  In  these  spinels. 

Various  compositional  modifications  were  made  to  the  pre-calcined 
spinel  powders.  After  additions  such  as  CdO,  ZnO,  or  Nb205  were  mixed  with 
the  appropriate  spinel  (by  vibratory  milling  followed  by  drying),  disks  were 
prepared  and  sintering  studies  were  carried  out.  Sintering  results  (weight 
loss  and  density)  and  x-ray  diffraction  data  are  presented  In  Table  3.3. 

Note  that  when  both  CdO  and  Nb205,  or  ZnO  and  Nb205,  were  added  simultaneously 
the  amounts  corresponded  to  an  addition  of  10  mole  percent  columbite.  As 
shown  by  the  data  In  Table  3.3,  excess  additions  of  ZnO  to  ZnCr204,  or  CdO  to 
CdCr204,  did  not  enhance  densi f Icatlon.  However,  enhanced  densi ficatlon  was 
achieved  In  all  cases  of  excess  Nb205,  regardless  of  whether  ZnO  or  CdO  was 
also  added.  For  the  case  of  ZnCr204,  a  1  weight  percent  addition  of  Nb205  was 
sufficient  to  Improve  the  densi ficatlon  of  the  spinel  ceramics,  whereas  a  10 
mole  percent  addition  of  columbite  corresponds  to  more  than  13  weight  percent. 
The  best  density  were  achieved  with  a  5  weight  percent  Nb205  addition  (4.95 


g/cc  (93%  of  theoretical).  For  the  case  of  CdCr204,  a  5  weight  percent  Nb205 
addition  resulted  in  a  sintered  spinel  ceramic  with  a  density  of  5.0  g/cc. 

It  was  found  that  if  CdO  loss  was  prevented,  by  using  a  closed  crucible  and  a 
CdO-rich  atmosphere  source  powder  (CdCr204  sand),  further  densification  to 
5.3  g/cc  (91%  theoretical)  could  be  achieved. 

X-ray  diffraction  data  of  the  sintered  samples  of  the  ZnCr204  spinel 
with  the  various  ZnO  and/or  Nb205  additions  Indicated  that  the  major  phase 
was  spinel,  and  a  trace  amount  of  columbite  ZnNb205  was  detected.  For  the  case 
of  CdCr204  spinel  with  additions  of  CdO  and/or  Nb205,  the  spinel  was  again 
the  major  phase,  but  minor  amounts  of  several  other  phases  were  detected: 
columbite  CdNb205,  pyrochlore  Cd2Nb207,  CrNb04,  Cr203,  and  Nb205.  The 
microstructures  of  sintered  disks  of  the  Cd  and  Zn  spinels  with  7  percent 
additions  of  Nb205,  are  presented  in  Figure  6.  The  microstructure  of  the 
CdCr204  spinel  sample  indicated  a  larger  grain  size  than  the  ZnCr204  sample, 
though  bith  samples  displayed  non-uniform  grain  structures,  with  grain  sizes 
ranging  from  1  to  5  microns. 

It  is  likely  that  the  mechanism  responsible  for  enhanced  densification 
of  the  spinel  ceramics  caused  by  excess  Nb205  is  identical  to  that  of  the 
columbite  additions.  Although  there  are  no  existing  phase  diagram  data  for 
the  CdCr204  and  ZnCr204  systems,  a  possible  explanation  for  the  enhanced 
sintering  can  be  inferred  from  existing  phase  diagrams  of  the  Cd0-Nb205  and 
Zn0-Nb205  systems.  These  phase  diagrams  are  shown  in  Figure  3.7.  Liquid 
phases  are  present  in  the  Cd0-Nb205  system  at  temperatures  as  low  as  1350®C, 
and  in  the  Zn0-Nb205  system  at  temperatures  below  1300®C.  It  is  possible 
that  the  presence  of  chromium  further  reduces  these  liquidus  temperatures. 

The  presence  of  a  liquid  phase  at  temperatures  in  range  where  a  liquid  phase 
is  present  would  suggest  that  the  Nb205  enhances  the  sintering  of  the  spinel 
ceramics  through  a  liquid-phase  sintering  mechanism.  It  is  unfortunate  that 
no  phase  diagram  data  exist  for  either  of  the  spinel -Nb205  systems.  To  verify 
the  existence  of  a  liquid-phase  sintering  mechanism,  differential  thermal 
analysis  (DTA)  was  carried  out  on  a  sample  of  CdCr204-9%Nb205  powder.  The 
DTA  pattern  is  presented  in  Figure  3.8.  An  endothermic  peak  was  observed  at 
a  temperature  of  1297®C,  confirming  the  presence  of  a  liquid  phase  at  this 
temperature.  This  endothermic  peak  was  not  observed  in  pure  CdCr204  powder. 

This  result  strongly  suggests  that  liquid-phase  sintering  is  the  mechanism 
whereby  the  Nb205  (or  columbite)  enhances  the  sintering  of  the  spinel  ceramics. 
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3.D  Summary  of  Results  of  Ceramic  Studies 


The  results  of  the  ceramic  studies  performed  In  this  Investigation 
are  summarized  below: 


1)  X-ray  diffraction  and  mlcrostructural  analyses  of  the  original 
CeramPhysIcs  ceramics  with  various  spinel /col umbite  ratios  were  performed. 
Different  microstructures  were  observed  between  the  Cd  and  Zn  analogues  of 
the  spinel /col umbite  ceramics. 

2)  Pure  CdCrpO^  could  not  be  densifled  (1n  air)  to  greater  than  80%  of 
theoretical,  because  of  excessive  CdO  loss.  The  denslfication  did  not 
improve  significantly  when  CdO  loss  was  limited  by  sintering  in  a  CdO-rIch 
atmosphere.  It  was  also  Impossible  to  density  pure  ZnCr20^  by  sintering 

In  air,  due  to  ZnO  volatilization.  However  90  percent  dense  ZnCr204  ceramics 
with  large  grain  size  (about  10  microns)  could  be  prepared  by  sintering  at 
1600®C  In  a  ZnO-rIch  atmosphere. 


3)  Denslfication  of  the  ZnCrpO^  ceramics  could  be  improved  by  a  one  weight 
percent.  Instead  of  the  13  weight  percent  columbite  (ZnNb20g)  addition 
that  has  previously  been  used.  Similarly,  CdCr204  ceramics  could  be 
densifled  by  a  5  weight  percent  addition  of  Nb205. 

4)  The  Improved  sintering  behavior  caused  by  the  columbite  or  Nb205 
additions  was  the  result  of  a  liquid  phase,  which  appears  to  be  present. 
This  conclusion  Is  supported  by  incomplete  phase  diagram  data  which  suggest 
the  presence  of  low-temperature  eutectics  In  the  Cd0-Nb205  and  Zn0-Nb205 
systems. 


Table  3.1a 


> 


1 


x-ray  Diffraction  Analysis  of  Spinel :Columbite  Composition. 
(CdCr20^) ;CdNb20g)  “  CeramPhysics 


Composition  Phases  Present 


CCN  (1/1) 

(#592) 

Cd2Nb202  (Pyrochlore) ,  Spinel 
Colofflbite, 

CCN  (3/1) 

(cal.  powder 
1050«C-2hr) 

Pyrochlore.  Spinel.  Nb20j 

CCN  (5/1) 

(#571) 

Spinel.  Pyrochlore,  Columbite 
^^2®3 

CCN  (7/1) 

(#563) 

Spinel.  Pyrochlore.  Columbite 

CCN  (9/1) 
(#312+6.75%  Fe) 

Spinel.  Columbite,  Pyrochlore 

CCN  (9/1) 

(#352) 

Spinel,  Columbite,  Cr203 

CCN  (11/1) 

(CPI-13  powder) 

Spinel,  Cr20^,  CdO,  Columbite 
(Note:  incompletely  reacted) 

Note:  All  samples  were  sintered  disks  unless 
otherwise  stated. 
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Table  3.1b 

X-ray  Diffraction  Analysis  of  Spinel:Columbite  Compositions. 
(ZnCr20^ :ZnNb205)  “  CeramPhysics 


Phases  Present 


Composition 


ZCN  (3/1) 

(mi) 

ZCN  (5/1) 
(#302) 


Table  3.1c 


■  X-ray  Diffraction  Analysis  of  Solid  Solution  Spinel :Columbite 

h 

k 

V 

Composition  Phases  Present 


ZnCr20^  Spinel,  CdNb20g  Columbite, 
CrNb04 

Cd2Zn^Cr2()Nb40j2  ZrCr204  Spinel,  CdNb20g  Columbite, 

CrNb04 


Sintering  Study  of  Zncr20^  and  CdCr20^  Disks 


Conpos it  ion 

Sintering  Condition 

Wt  Change  % 

Density  g/cc 

ZnCr20^ 

1290«C/1  hr 

+1-2% 

3 

1340»C/1  hr 

+  1-2% 

3 

1420»C/1  hr 

-0.6% 

3.4 

1550»C/1  hr 

-2% 

3.7 

1650*C/1  hr 

-3% 

3.8 

CdCrjO^ 


1290«C/1  hr 
1350«C/1  hr 


+1% 

-13% 


4.4 
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Figure  3.1c,  i'.-ray  fluorescence  of  Cr  in  CdCr.,0 


CCN  5 ■ 1 


2SKU  433X  seu  0029 


c.}'  f louresccrsce  of  :!'d  in  CcIC'.’oO 


Figure  3.  Ir 


Fioure  3.3b. 


;  j,  ..  . 

;  ■  ''v' 

f  ^ — ■  »  ■. 

i  ^  '  V  .  ' 

'  V  ■  ' »  ,  ' 

V--  ■  \  -  ^  ' 

.....  •  ( , 

■  '  ‘  'd 

r'-'  ’  ■  _ 

f  •_-  * 

t  ■  *  ^  '■■  ■'  s  ■■' 

.  .  . 

■  ’  '  ■••••'  »sV  . 

■  V  ..  ... 

2SKU  166X  leeii  eezs 

r--; 

’  'r r 
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Figure  3.3d. 
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Figure  3.5a.  Microstructure  of  pure  CdCr  0,  sample  sintered  in 
CdCr^O,  sand  at  1300°C  for  6.^  hrs. 
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Figure  3.5b.  Microstructure  of  CdCr  0  /CdNb  0  ( 9  / 1  )  s.irnple  sintort 
at  1300°C  for  one  hour. 
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Figure  6h.  Microstructure  of  7,nCr^O,/77'  Nb^O^  sample  sintered  at 
for  1 . a  hours . 
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Figure  3.7a.  System  Zn0-Nb20 » "Z^N 
‘^2,^1.^  ~  ^^2^^34^87* 
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Figure  3.7b.  Cd0-Nb20j  phase  diagram. 
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IV.  SPECIFIC  HEAT,  MAGNETOCALORIC,  AND  DIELECTRIC  CONSTANT 


MEASUREMENTS  AT  LOW  TEMPERATURES 

This  Section  is  devoted  to  experimental  research  on  the 

specific  heat,  magnetocaloric,  and  dielectric  constant  properties 

of  CdCr-,0.  -  and  ZnCr^O.  -  containing  ceramics  at  low 
2  4  2  4 

temperatures  and  in  intense  magnetic  fields.  All  measurements 
here  were  made  on  dense  or  semi-dense  ceramic  pellets,  and  the 
nomenclature  used  below  is  as  follows: 

CCN(9/1)  =  Mixture  of  10  mole%  CdNb20g,  90  mole%  CdCr20^, 
fully  dense,  fine  grain  size 

ZCN(9/1)  =  Zinc  analog  of  CCN(9/1) 

CdCr.O.  =  Pure  CdCr_0.,  74%  dense,  fine  grain  size 
2  4  2  4 

ZnCr20^  =  Pure  ZnCr20^,  72%  dense,  fine  grain  size 

This  Section  is  organized  along  the  following  lines:  First, 
the  dielectric-constant  data  will  be  presented,  followed  by  the 
specific  heat  data.  Next,  the  specific  heat  data  will  be 
analyzed.  Thirdly,  the  magnetocaloric  data  will  be  presented, 
and  lastly  these  data  will  be  analyzed. 

A.  Dielectric  Constant  Measurements 

Dielectric  measurements  of  CCN(9/1)  and  ZCN(9/1)  disks  have 
been  made  at  low  temperatures.  Each  sample  had  sputtered  gold 
electrodes  and  was  a  disk  approximately  one  centimeter  in 
diameter  and  one  millimeter  thick.  The  measurements  were  made  in 
a  dielectric  probe  where  the  lead  capacitance  was  on  the  order  of 
a  tenth  of  a  picofarad  (about  1%  of  the  measured  values). 

Dielectric  constant  versus  temperature  data  for  each  sample 
are  shown  in  Fig.  4-1.  In  each  case  the  dielectric  constant 
peaks  above  the  temperature  at  which  the  specific  peaks  occur  [8 
K  and  10.7  K  for  CCN(9/1)  and  ZCN(9/1),  respectively].  The 
specific-heat-peak  temperatures  are  shown  as  the  vertical  arrows. 

From  Fig.  4-1,  it  appears  that  the  slope  of  the  dielectric 
constant  increases  at  the  specific-heat-peak  temperatures.  To 


CCN(9/1) 

ZCNO/I) 


Temp  (K) 


test  this  observation,  the  average  slope  (Ae/AT)  was  calculated 
between  each  pair  of  data  points  and  is  plotted  in  Fig.  4-2(a) 
where  the  temperature  of  each  point  is  the  average  temperature  in 
that  interval.  Figure  4-2(a)  clearly  shows  the  dramatic  in¬ 
creases  in  the  slopes  very  near  the  specific-heat-peak  tempera¬ 
tures  . 

Figure  4-2(b)  shows  the  variation  in  the  dielectric  constant 
as  a  function  of  frequency  at  constant  temperature  near  the  peaks 
of  each  material  (8.01  K  and  10.55  K,  respectively).  Both  sets 
of  data  are  normalized  to  the  dielectric  constant  at  200  Hz. 
Although  both  changes  are  small,  the  ZCN(9/1)  changes  are  much 
larger  than  the  CCN(9/1)  changes. 

Following  the  discovery  of  sharp  maxima  in  Ae/AT  at  the 
temperature  of  the  specific  heat  maximum  for  both  CCN(9/1)  and 
ZCN(9/1),  Fig.  4-2(a),  a  decision  was  made  to  pursue  these 
measurements  on  the  pure  materials  and  in  intense  magnetic 
f  ields . 

Figure  4-3  shows  the  original  dielectric  data  for  CCN(9/1) 
and  the  new  dielectric  data  for  pure  CdCr20^  in  zero  field.  The 
latter  sample  was  a  disk  that  was  only  54%  of  theoretical  densi¬ 
ty,  and  the  data  have  been  corrected  for  this  low  density 
following  a  formula  given  by  Niesel  (1952).  The  pure  material 
has  a  smaller  dielectric  constant  than  the  CCN(9/1)  and  does  not 
have  the  rapid  rise  associated  with  the  spinel/columbi te  mixture. 
At  first  glance,  it  appears  there  might  be  a  small  structure  in 
the  pure  CdCr20^  data  at  8  K,  but  upon  repetition  of  the 
measurements  across  this  region  (not  shown)  it  became  clear  that 
this  apparent  structure  is  near  the  noise  level  and  is  probably 
not  present.  Except  for  the  broad  maximum,  the  dielectric 
constant  of  the  pure  material  is  distinctly  different  from  that 
of  CCN (9/1 ) . 

Figure  4-4  shows  similar  data  for  ZCN(9/1)  and  pure  ZnCr20^ 
where  again  the  pure  material  data  have  been  density-corrected 
(70%  of  theoretical).  Here  the  pure  material  does  show  slight 
structure;  i.e.,  a  small,  comparatively  broad  rise  in  the 
dielectric  constant  between  6  and  8  K.  There  is  no  broad  peak 
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•  CCN(9/1)  ■ZCNO/l) 


Figure  4-2.  (a)  Variation  of  tite  derivative  &e/AT  witli  temperature  in 

the  neighborhood  of  the  specific  heat  maxima.  (b)  Frequency 
dependence  of  e  near  the  temperature  of  the  specific  heat 
maxima . 
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Figure  4-A.  Comparison  of  dielectric  data  measured  on  ZCN(9/1)  and  ZnCr  0, 
The  structure  in  the  latter  data  is  small,  but  this  material 
exhibits  a  substantial  maximum  in  the  specific  heat  (see  Fig. 
4-  10) . 


above  this  rise  as  in  ZCN(9/1).  The  anomalous  rise  in  the 
dielectric  constant  coincides  with  the  peak  in  the  specific  heat 
for  the  spinel  +  columbite  samples.  There  is  no  correlation  in 
either  of  the  pure  materials  between  the  specific  heat  data  and 
dielectric  constant  data  (see  Figs.  4-9  and  -10  below). 

The  three  materials  which  show  dielectric  structure  [i.e., 
CCN(9/1),  ZCN(9/1)  and  ZnCr20^]  were  measured  at  the  National 
Magnet  Laboratory  at  MIT.  In  each  case,  the  following  procedure 
was  followed.  At  zero  field,  the  dielectric  constant  was 
measured  as  a  function  of  temperature  through  the  range  of  the 
dielectric  structure.  Then  at  a  constant  temperature  slightly 
below  the  maximum  in  3e/3T  the  dielectric  constant  was  measured 
as  a  function  of  H-field.  Holding  the  field  constant  at  15  T, 
another  temperature  sweep  was  made  over  the  same  temperature 
range  as  above. 

At  zero  field,  the  temperatures  were  held  constant  using  a 
capac i tance- temperature  controller  sensing  a  field-independent 
capacitance  thermometer  (Lawless,  1971)  in  the  probe.  This 
thermometer  was  calibrated  in  zero  field  against  a  germanium 
thermometer  during  the  first  sweep.  Then  at  H  /  0,  the  con¬ 
troller  could  be  set  to  the  same  known  temperature  points. 

In  the  following  data,  it  will  be  noted  that  the  absolute 
dielectric  constants  are  slightly  higher  in  each  case  compared  to 
the  original  data  of  Figs.  4-3  and  4-4.  This  is  due  to  a  slight¬ 
ly  larger  lead  capacitance  in  the  MIT  probe  (~  0.1  -  0.2  pF).  No 
correction  was  made  for  this  effect.  In  each  case,  the  samples 
were  identically  the  same,  and  the  dielectric  constants  of  the 
pure  spinels  were  corrected  as  above  for  density. 

Figures  4-5,  4-6,  and  4-7  show  the  results  of  these  measure¬ 
ments  where  in  each  figure  the  main  graph  shows  the  results  of 
the  two  temperature  sweeps  at  0  T  and  15  T  and  the  inset  shows 
the  field  dependence  of  the  change  in  dielectric  constant  at 
constant  temperature.  For  both  the  pure  ZnCr204  and  the 
CCN(9/1),  the  dielectric  constant  is  consistently  higher  across 
the  full  temperature  range,  while  the  dielectric  constant  of 
ZCN(9/1)  at  15  T  compared  to  0  T  is  first  above  at  low 
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Figure  4-3.  Magnetic  field  dependence  of  the  dielectric  constant  of  CCN(9/1) 
in  the  neighborhood  of  the  specific  tieut  maximum  at  8  K. 


temperatures,  then  below  at  high  temperatures.  The  inset  graphs 

2 

indicate  that  Ae  in  CCN(9/1)  is  proportional  to  H  whil-  Ae  in 
both  of  the  Zn  samples  is  proportional  to  H  with  the  slope 
approximately  the  same  in  both  cases  [1.09x10  ^  (1/T)  for 
2CN(9/1)  and  1.26xl0~^  (1/T)  for  pure  ZnCt^O^^]. 

The  derivatives  Ae/AT  for  the  temperature  sweeps  in  a  mag¬ 
netic  field  are  plotted  in  Fig.  4-8  as  a  function  of  temperature 
for  both  of  the  sp i ne 1/col umbi te  samples.  The  peak  in  Ae/aT  in 
ZCN(9/1)  has  not  shifted  in  field  while  the  peak  for  CCN(9/1)  has 
shifted  by  approximately  0.5  K.  From  ptevious  woik,  it  has  .oaen 
found  that  the  temperatures  of  the  specific  heat  peaks  of 
CCN(9/1)  and  ZCN(9/1)  do  not  shift  in  field. 

It  is  an  open  question  what  role  the  columbi te  phase  plays 
in  each  case;  it  is  possible  that  denser  samples  of  the  pure 
spinels  would  more  closely  mimic  the  sp i ne 1/co  1 umb i te  behavior. 
This  possibility  is  suggested  by  Figs.  4-3  and  4-4  where  the  more 
dense  (70%)  ZnCr20_^  showed  structure  closer  to  the  9/1  samples 
than  the  less  dense  (54%)  CdCr^O^. 

One  possibility  to  consider  is  a  magnetic  phase  transition 
occuring  in  the  spinel  component  accompanied  by  a  volume  change 
which  in  turn  applies  a  pressure  to  the  columbi te  phase  changing 
the  dielectric  constant  of  the  columbite  phase  by  e lec t ros t r i c- 
tion.  The  dielectric  changes  are  then  only  indirectly  probing 
the  phase  transition  in  the  spinel.  To  analyze  this  situation, 
consider  the  free  energy  of  the  columbite  phase  as  a  function  of 
polarization  P  and  hydrostatic  pressure,  p,  in  the  Ginzburq- 
Landau -De vonsh i re  formalism, 

A(P,p)  =  +  ^CP"*  ...  +  5(0^^  +  ^  ••• 

where  and  e  lec  tros  tr  ic  t  i  ve  coefficients.  For  P  =  0, 

Xp  =  4n/£p  =  8^A/9P^  =  ^  2C)^2)P  +  ••• 


Therefore , 


ZCN(9/1) 
o  H  =  0.0  T 
•  H  =  15.16  T 


where  q  =  0^^  +  ^*^12  subscripts  on  c  refer  to  pressure. 

Here,  p  in  Eq.(4-1)  is  the  pressure  exerted  on  the  columbite 
phase  by  the  spinel  phase. 

Now  consider  three  general  types  of  phase  transitions  in  the 
spinel  phase  and  use  Eq.(4-1)  to  predict  the  behavior  of  e(T). 

In  a  first-order  transition  the  pressure  will  be  a  step  function 
with  a  cut-off  at  the  transition  temperature  From  Eq.(4-1), 

£  will  also  change  in  a  step  pattern.  Although  there  is  a  sudden 
rise  in  the  dielectric  constant  of  both  (9/1)  mixtures,  this 
model  does  not  explain  either  the  maxima  in  c  in  the  (9/1)  sam¬ 
ples  nor  the  rise  in  e  in  the  pure  ZnCr^O^  sample. 

Next  consider  a  second-order  transition  with 

p  -  (T^  - 

Using  Eq . ( 4-1 ) , 

3e  /3T  <x  c^/(T^  - 
P  p  N 

Here  3g/3T  stays  finite  at  because  of  the  [.-isence  of  higher 
order  terms  in  the  free  energy  expansion.  Again,  this  case  does 
not  correspond  to  the  experimental  results. 

The  third  case  is  a  distributed  phase  transition  with 
mechanical  relaxation  where  approx ima tely 

p  «  exp[-a{Tj^  -  T)^] 

Then,  using  Eq.(4-1): 

3c/3T  -  (Tj^  -  T)  exp(-a(Tj^  -  T)^]. 

In  this  case  an  extremum  occurs  in  c  at  T  =  T^^,  and  a  change 
in  sign  is  predicted  for  Ae/AT  at  T  =  T^. 

None  of  these  cases  fits  the  experimental  observations,  and 
we  are  left  with  the  conclusion  that  the  columbite  phases  in 
CCN(9/1)  and  ZCN(9/1)  probably  do  not  play  a  role  in  these 


dielectric  anomalies.  Consequently,  these  dielectric  anomalies 
are  associated  with  a  coupling  between  the  optic  modes  and  spin 
systems  within  the  spinel  phase. 


B .  Specific  Heat  Measurements 

As  mentioned  in  Section  III  above,  compacted  disks  of  the 
pure  spinel  powders  were  successfully  formed  at  1290°C  for  1  hr 
in  the  case  of  CdCr^O^,  1650°C  for  1  hr  for  ZnCr20^.  Weight 
changes  were  negligible  (~  +1  and  -3%,  respectively). 

Specific  heat  samples  were  prepared  as  follows:  Two  disks 
of  CdCr20^  were  first  grooved  across  the  diameter,  and  these 
grooved  faces  were  then  bonded  together  to  form  the  basic  sample. 
Similarly  for  ZnCr^O^.  The  grooved  faces  thus  formed  a  slot  for 
the  insertion  of  a  carbon-chip  thermometer.  A  generous  amount  of 
copper  foil  was  used  in  the  sample  assembly  (e.g.,  between  bonded 
faces,  on  the  disk  periphery,  etc.)  to  provide  good  thermal 
diffusion,  as  later  proved  to  be  the  case  in  the  actual 
experiments.  Nonetheless,  the  specific  heats  of  the  disks  were 
large  enough  that  the  addenda  contributions  were  very  small 
( <  5% )  despite  the  copper- foi 1  additions. 

The  densities  of  the  starting  disks  were  as  follows: 

CdCr204:  4.359±0.009  (74.4±0.1%  of  theoretical) 

ZnCr^O^:  3.865±0.021  (72.0+0.4%  of  theoretical), 

using  the  NBS  XRD  densities  (3.862  and  5.367,  respectively).  It 
is  worth  noting  that  the  density  of  a  f ree  powder  is  ~  65-68% 
(sphere  close-packing),  so  the  dens i f ica t i on  of  these  disks  is 
only  marginally  larger  than  that  of  the  free  powders.  It  should 
also  be  pointed  out  that  these  disks  proved  surprisingly  hard  to 
cut,  despite  the  marginal  dens i f ica t ion . 

Specific  heat  data  were  measured  2-34  K  by  the  dynamic  pulse 
technique  in  the  adiabatic  calorimeter,  and  the  two  samples  were 
measured  simultaneously. 

Specific  heat  data,  C  vs.  T,  for  the  CdCr20^  disk  are  shr'wn 
in  Fig.  4-9  compared  to  the  specific  heat  of  CCN(9/1).  In  the 


spinel  disk,  a  small,  bfoad  maximum  in  C  occurs  at  5.8  K,  in  good 
agreement  with  the  x~data  in  Fig.  5-x  below.  However,  this 
maximum  is  smaller  and  much  broader  than  the  C-maximum  in 
CCN(9/1)  at  8.0  K  in  Fig.  2-1. 

The  corresponding  specific  heat  data  for  the  ZnCr^O^  spinel 
and  ZCN(9/1)  are  shown  in  Fig.  4-10.  Here  a  different  picture 
emerges,  as  the  spinel  disk  has  a  distinct  peak  in  C  at  9.7  K, 
compared  to  the  C-peak  in  ZCN(9/1)  at  10.5  K.  The  temperature  of 
the  C-maximum  for  the  ZnCt^O^  spinel  is  in  good  agreement  with 
the  x~data  in  Fig.  5-y  below. 

3 

The  behavior  of  C/T  often  provides  useful  insights,  and  in 
Figs.  4-11  and  4-12  are  shown  C/T^-data  for  the  CdCi^O^  and 
ZnCr20^  materials,  respectively.  In  both  of  these  plots  are 
shown  data  for  the  disks  measured  here  plus  separated  data  for 
the  spinel  phases  in  CCN(9/1)  and  ZCN(9/1)  from  Fig.  2-1.  In  all 
cases  in  Figs.  4-11  and  4-12,  there  is  a  rapid  rise  in  C/T^  with 
decreasing  temperature  at  the  lowest  temperatures,  and  this  has 
previously  been  attributed  to  a  Schottky  term  due  to  f ree  spins . 
Using  this  interpretation,  one  would  say  from  Fig.  4-11  that  the 
density  of  free  spins  in  the  CdCr20^  compacted  disk  is  ~  6-7 
times  larger  than  for  the  spinel  phase  in  CCN{9/1).  However, 
this  picture  does  not  agree  with  the  x“data  in  Fig.  5-x  below. 
From  Fig.  4-12  for  the  ZnCr20^  materials,  one  would  speculate 
that  the  density  of  free  spins  in  the  ZnCr20^  compacted  disk  is 
^  2  times  larger  than  for  the  spinel  phase  in  ZCN(9/1),  and  this 
interpretation  is  not  necessarily  out  of  line  with  the  x'^^^ta  in 
Fig.  5-y  below. 

It  is  interesting  to  observe  how  large  the  specific  heat  of 

the  CdCr20^  disk  really  is  at  the  lowest  temperatures.  We 

previously  determined  a  reliable  Debye  temperature  for  spinel 

CdCi’TO.  ,  0p,  =  420  K,  and  this  means  that  the  Debye  limit  is 
^  ^  ^  -1-4 

(C/T  )  =  6.55  erg  g  K  .  From  Fig.  4-11  for  the  CdCi-  0.  disk 

u  3  3 

at  the  lowest  temperatures,  C/T  =  3x10  --  i.e.,  the  Debye 

contribution  is  only  ~  0.2%. 

It  is  difficult  to  reconcile  the  C-data  reported  here  and 
the  x'data.  First,  given  the  very  slight  dens  i  f  icat  ion  of  these 


spinel  disks,  it  is  reasonable  to  consider  these  disks  and  the 
respective  powders  as  being  essentially  identical  (e.g.,  same 
surface  area).  Going  further,  the  grain  sizes  in  the  powder  and 
disk  are  most  likely  also  identical,  for  the  same  reason. 

Electron  microscope  studies  are  needed  here.  The  broad  C-peaks 
in  Figs.  4-9  and  4-10  may  be  due  to  short  correlation  lengths  in 
the  disks,  since  it  is  believed  that  the  spinel  grains  in 
CCN(9/1)  and  ZCN(9/1)  are  quite  large  {~  lOum)  and  have  long 
correlation  lengths,  leading  to  sharp  C-peaks. 

However,  the  x~data  in  Section  V  below  consistently  show 
that  the  density  of  free  spins  increases  as  the  surface  area 
decreases  (e.g.,  powders  have  few  free  spins).  On  the  other 
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hand,  the  very  large  value  of  C/T  for  the  CdCr20^  disk  at  the 
lowest  temperatures  in  Fig.  4-11  suggests  a  very  large  density  of 
free  spins. 

Magnetic  Field  Dependence.  The  magnetic  field  dependence  of 
the  specific  heats  of  the  CdCr204  and  ZnCr20^  disks  were  measured 
in  a  field  of  10  T  at  the  National  Magnet  Laboratory.  A  so- 
called  drift  technique  developed  by  CeramPhysics  was  used  in 
these  measurements  (Lawless  e^  al.,  1982). 

Magnetic-field-dependent  specific  heat  data  in  the  neighbor¬ 
hood  of  T.,  for  the  CdCr_0.  and  ZnCr~0.  disks  are  shown  in  Figs. 
4-13  and  4-14,  respectively,  at  H  =  10  T  and  H  =  0.  In  the 
paramagnetic  regions,  the  10  T  field  decreases  the  specific  heat 
of  CdCr20^  by  ~  35%,  Fig.  4-13,  and  this  is  thermodynamically 
consistent  with  the  magnetocaloric  phase  diagram  shown  in  Fig. 
4-24  below.  That  is,  based  on  simple  entropy  arguments,  magneti¬ 
zation  heating  requires  that  a  magnetic  field  suppresses  the 
specific  heat. 

However,  a  magnetic  field  increases  the  specific  heat  of 
ZnCr20^,  Fig.  4-14,  and  this  is  thermodynamically  i neons isten t 
with  the  magnetocaloric  data  in  Fig.  4-30  below. 

Our  interpretation  here  is  that  the  carbon-chip  thermometer 
on  the  ZnCr20^  sample  probably  sh i f ted  in  calibration  in  the  10  T 
field,  leading  to  the  inconsistent  results  in  Fig.  4-14.  Along 


this  line,  this  thermometer  was  observed  to  shift  in  calibration 
during  the  magnetocaloric  measurements  (see  i^elow).  Further 
support  for  this  interpretation  comes  from  the  fact  that  a  field 
suppresses  for  an  anti  ferromagnetic  transition,  as  in  Fig. 
4-13,  whereas  the  opposite  is  seen  in  Fig.  4-14. 

This  interpretation  would  seem  to  cast  some  doubt  on  the 
ZnCr20^  magnetocaloric  measurements  below;  however,  there  are 
major  differences  in  the  two  experiments.  First,  calibration 
points  were  taken  at  practically  all  magnetocaloric  set  points, 
and  there  was  internal  consistency  amonst  the  sets  of  magneto¬ 
caloric  data.  And  second,  magnetocaloric  measurements  involved 
AT's. 

Data  Analyses.  We  now  consider  phenomenological  analyses  of 
the  specific  heat  data  for  the  pure  spinels  shown  in  Figs.  4-9 
through  4-12  for  comparison  with  similar  analyses  of  the  data  for 
the  spinel  phases  in  CCN(9/1)  and  ZCN(9/1)  (Eckels  ^.  ,  1985). 
To  review,  the  spacing  of  quantized  energy  levels  gives  rise  to  a 
so-called  Schottky  specific-heat  term.  For  a  two-level  system 
where  5  is  the  energy  separation,  for  example,  the  high-tempera¬ 
ture  (T  >>  6 )  form  of  the  Schottky  term  is 

Csch  =  nRgQg^(go  +  g^)”^  i6/T)^  (4-2) 

where  gg  and  g^  are  the  degeneracies  of  the  two  levels,  and  n  is 

the  number  of  level  systems  per  formula  weight.  The  form  of 
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Eg. (4-2),  “  T  ,  is  correct  for  the  genera  1  Schottky  term, 

and  in  the  Debye  limit  for  the  lattice  contribution  (i.e., 

T^)  a  plot  of 

CT^  =  mT^  +  b  (4-3) 

is  a  straight  line  where  m  is  related  to  the  Debye  temperature 
(0j^)  and  b  /  0  is  the  coefficient  of  the  Schottky  term. 

A  subtlety  enters  the  m-coe f f ic ien t  in  Eq.(4-2),  as  follows: 
If  an  anti  ferromagnetic  transition  occurs,  then  at  temperatures 


below  a  T  spin-wave  contribution  to  the  specific  heat  is 

present,  and  in  Eq.(4-3)  m  =  where  is  related  to  the 

Debye  temperature  and  m^^  is  the  coefficient  of  the  spin-wave 

term.  For  a  ferromagnetic  transition,  the  spin-wave  contribution 
^3/2 

varies  as  T  . 

Thus,  these  Schottky  analyses  of  specific  heat  data  yield  a 

great  deal  of  inf onnation :  (1)  The  b-coef f icient  in  Eg. (4-3) 

2 

yields  data  for  n6  from  Eg. (4-2),  given  certain  assumptions 
regarding  and  g^^;  and  (2)  The  m-coeff  icient  in  Eg.  (4-3)  yields 
information  on  both  the  Debye  temperature  and  the  presence  of 
anti  ferromagnetic  spin  waves. 


tures  for  the  spinel  phases  in  the  spinel  +  columbite  densified 
ceramics  CCN(9/1)  and  ZCN(9/1). 

Examples  of  these  types  of  Schottky  analyses  have  previously 
been  performed  for  the  CdCr20^  and  ZnCr20^  spinel  phases  in  the 
densified  ceramics  CCN(9/1)  and  ZCN{9/1),  respectively,  and  these 
analyses  are  reproduced  in  Fig.  4-15  (Eckels  e^  al . , 1985 ) . 
Analyses  at  T  >  T^^,  Fig.  4-15(b),  are  not  complicated  by  spin- 
wave  terms,  whereas  at  T  <  T^,  Fig.  4-15(b),  this  complication 
may  enter.  Excellent  fits  to  Eg. (4-3)  are  obtained  in  Fig.  4-15 
with  b  /  0  (note  that  the  deviations  at  the  lower  temperatures 
are  due  to  the  T  >>  6  approximation),  and  we  note  in  particular 


that  the  excellent  T  <  fits  in  Fig.  4-15(3)  are  strong  evi¬ 
dence  for  an t i f e r romagne tic  transitions.  We  shall  retain  below 
to  the  details  of  these  fits. 

Returning  to  Figs.  4-11  and  4-12,  we  note  that  the  behavior 
of  C/T^  is  basically  the  same  at  T  >  T,,  for  the  spinels  whether 

-  tN 

in  the  compacted  disks  or  as  the  spinel  phases  in  the  densified 

spinel  +  columbite  ceramics.  Likewise,  at  T  <  T  there  is  a 

3  ^ 

marked  rise  in  C/T  with  decreasing  temperature  in  all  cases. 

Therefore,  by  analogy  with  the  data  in  Fig.  4-15,  Schottky 

analyses  at  T  <  and  T  >  were  carried  out  for  the  specific 

heat  data  measured  on  the  compacted  disks  of  CdCi^O^  and  ZnCr20^ 

and  these  data  are  shown  in  Figs.  4-16  and  4-17,  respectively. 

Once  again,  excellent  plots  according  to  Eg. (4-3)  are  obtained. 

The  results  of  leas t-sguares  fittings  of  the  data  for  all 

materials  are  summarized  in  Table  4-1.  In  all  these  analyses, 

the  assumption  g^  =  g^^  in  Eg.  (4-2)  has  been  made. 

Table  4-1 

Schottky  Analyses,  Eg. (4-3)* 


CdCr  .,0^ 

ZnCr ^0^ 

Quantity 

In  CCN(9/1) 

Disk 

In  ZCN(9/1) 

Disk 

m  (High-T) 

6.546 

6.546 

5.875 

6.765 

m  (Low-T) 

470.6 

2764 

136.1 

2182 

(High-T 

420 

420 

463 

442 

9p  (Low-T) 

101 

56 

162 

139 

n5 ^  (High-T) 

1679 

2168 

1264 

1377 

n6  ^  (  Low-T  ) 

0.109 

0.761 

0.135 

0.637 

n_C. (k/2Js)  ^ 

1 . 57x10"^ 

9.30x10”^ 

3.66x10"'^ 

6.11x10 

*Units  are  consistent  with  the  specific  heat  in  units  of 
erg  g  K  ;  note  also  that  g^  =  g^^  in  Eg.  (4-2). 

Also  included  in  Table  4-1  are  data  for  the  spin-wave  specific 
heat  term. 


=  n_  R  C_  (k/2Js)^T^ 


-6C- 


_ ^'t-17.  SuhiUtky  fits  to  the  '1'  >  specifie  heat  data  fof 

ZnCr.O  aceording  to  I'.q  .  ( i )  . 


computed  from  the  (iitfetences  in  the  m-coe  1 1  i  c  i  e  n  t  s  i  i  .  , 

m  =  m^  +  m  =  m^^  +  n  R  C  (-: /2J5)  1,  whore  n  is  t:.o  density  ot 

spin  waves  . 

A  rather  staggering  amount  of  experimental  data  analyses  art- 
summarized  in  Table  4-1,  an  i  these  results  fall  into  s--veral 
separate  categories,  as  follows: 

Debye  Temperatures.  The  agreement  between  the  9q's  for  the 
spinels  from  the  h i gh- tempo ra tu re  analyses  are  very  good  (exact 
in  the  case  of  CdCr^O^).  For  comparison,  the  Lindemann  equation 
relates  6^  to  other  properties. 

e,,  s  B(T  (4-5) 

D  m 

where  T^  is  the  melting  point,  M  is  the  molecular  weight,  V  is 
the  molecular  volume,  and  B  is  constant  (=  115).  Adopting  ISOO^C 
melting  points  for  both  spinels/  we  find  from  Eg.  (4-5)  that 

(CdCr20^)  =  414  K 
(  ZnCr20^  )  =  469  K 

in  remarkably  good  agreement  with  the  data  in  Table  4-1  [Note  in 
this  regard  that  changing  by  ±  200°C  in  Eg. (4-5)  changes  9^  by 

only  ±  5%J.  We  also  point  out  that  these  0p.'s  are  large  enough 
to  justify  the  T  approximation  for  the  lattice  specific  heats 
implicit  in  Fig.  4-17. 

The  Deybe  temperatures  in  Table  4-1  from  the  low-temperature 
Schottky  analyses  are  clearly  too  small,  and  this  is  firm 
evidence  for  anti  ferromagnetic  spin  waves. 

Schottky  Terms,  T  >  T^j.  The  Schottky  specific  heat  term  in 

the  paramagnetic  phase  (T  >  T^^ )  is  due  to  the  spins  that  partake 

in  the  transition  at  T.,,  and  it  seems  reasonable  to  assume  that  -5 

N 

is  the  same  for,  say,  ZnCr2b)4  in  both  the  compacted  disk  and  as 
the  spinel  phase  in  the  spinel  columbite  densified  c<?ramic. 


This  argument  can  be  generalized  further  by  observing  that 

whatever  the  form  of  the  actual  Schottky  term,  the  h igh-tempera- 
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ture  expansion  will  always  have  the  functional  form,  n/T  , 

where  n  is  the  density  of  contributing  spins.  Thus,  6  becomes  a 
constant  of  proport  iv^na  1  i  ty  which  is  assumed  the  same  for  the  two 
forms  of  ZnCi'  O  (and  CdCr  O .  )  .  Proceeding  along  this  line,  the 
ratio  of  the  n5  values  should  correspond  to  the  ratio  of  spin 
densities  in  the  two  forms,  and  from  Table  4-1  we  find: 

Table  4-2 


Spinel 


CdCr^O 
ZnCr„0 


Paramagnetic  Spin  Densities,  T  >  T^ 

n ( Disk )/n( spinel  +  columbite) 


There  is  satisfactory  agreement  between  the  contributing 
spin  densities  in  the  two  forms  of  ZnCr204  in  Table  4-2  (i.e.,  ~ 
t  5%),  and  along  this  line  we  note  the  similarities  in  the 
specific  heat  maxima  in  Fig.  4-10  where  the  two  curves  appear 
shifted  by  about  1-2  K.  On  the  other  hand,  the  magnetocaloric 

phase  diagram  at  10  T  in  Fig.  4-31  below  for  the  two  forms  of 
ZnCr20^  suggests  that  the  spin  density  in  the  compacted  disk  is  ~ 
2  times  larger  than  the  spinel  +  columbite  ceramic.  However,  the 
magnetocaloric  data  reflect  the  total  spin  density  whereas  the  T 
>  Tj^  Schottky  term  reflects  the  spins  partaking  in  the  transition 

For  CdCr_,0.,  the  Table  4-2  data  indicate  that  about  30%  more 
2  4 

spins  partake  in  the  transition  in  the  compacted  disk  than  in  the 
sp ine 1 : columbi te  ceramic. 


Schottky  Terms,  T  <  T^^.  Below  the  transition  it  is 
reasonable  to  assume  that  the  f ree  spins  give  rise  to  the 
Schottky  term.  Proceeding  as  above,  we  find  for  the  ratio  of 
free-spin  densities: 
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Table  4-3 

Free  Spin  Densities,  T  <  T^j 
_ n(Disk)/n  (Spinel  +  Columbite) 


CdCr-0 . 
2  4 


ZnCr^O^ 


6.98 

4.72 


Here,  we  find  a  wide  disparity  in  the  free-spin  densities 


between  the  two  forms  of  both  CdCr20^  and  ZnCr20^.  The  origin  of 


these  disparities  is  simply  the  very  large  differences  between 
the  specific  heats  of  the  compacted  disks  and  the  spinel  + 
columbite  ceramics  in  Figs.  4-9  and  4-10. 


Our  interpretation  of  these  T  <  Tj^  Schottky  terms  may  be 


incorrect,  for  two  complementary  reasons:  (1)  The  susceptibility 
data  clearly  show  that  the  spinel  +  columbite  ceramics  have  much 
larger  free  spin  densities  compared  to  the  powders  (see  Section 
V ) ;  and  (2)  The  magnetocaloric  phase  diagrams  below  support  these 
susceptibility  data. 


The  question  then  is  the  origin  of  these  T  <  T^^  Schottky 


terms,  since  in  the  simple  model  of  the  two-spin  systems  the 
correlated  spins  contribute  the  T^  spin-wave  term,  and  the 


uncorrelated  spins  contribute  the  Schottky  term.  Possibly  a 
third  spin  system  is  involved  as  suggested  by  some  of  the 
susceptibility  data  in  Section  V. 


Spin-Wave  Densities.  As  a  final  exercise,  we  consider  the 


densities  of  the  spin  waves,  n  ,  from  Table  4-1,  as  given  py  the 

3  ® 

ratios  of  n  C  {k/2Js) 

S  3 


,  and  we  recall  that  these  densities  come 
from  the  m-coe f f i c ien ts  in  Eg. (4-3)  and  are  independent  of  the 
Schottky  terms.  Proceeding  as  above,  we  find 


Table  4-4 

Spin  Wave  Densities 


Spinel 


(Disk) /n( spinel  +  columbite) 


CdCr-,0  . 
2  4 


-s 
5.92 


ZnCr  ^0 . 
- 2-4- 


16.7 


Here  also  we  find  a  wide  disparity  in  n^  from  the  two  forms 
of  these  spinels.  Recalling  that  n^  is  also  the  density  of 
an t i f e r romagne t ica 1 ly  correlated  spins  and  that  these  spins  give 
rise  to  magnetization  cool i ng  (see  below),  the  Table  4-4  data 
suggest  that  such  cooling  phenomena  should  be  much  more 
pronounced  in  the  compacted  disks  than  in  the  spinel  +  columbite 
ceramics  --  but  this  is  not  borne  out  in  the  magnetocaloric  phase 
diagrams  below.  Along  this  line,  however,  we  do  note  from  Table 
4-1  that  the  quantity  n  C  (k/2Js)^  is  larger  for  CdCr^O.  than 
for  ZnCr20^,  and  magnetization  cooling  is  always  seen  in  the 
C6Qr 2'^ ^  materials  but  not  in  the  ZnCr20^  materials,  as  we  shall 
see  shortly. 

The  data  and  analyses  reported  here  support  the  earlier 
findings  on  the  CCN(9/1)  and  ZCN(9/1)  ceramics  and  lead  to  a 
clear  picture  of  the  specific  heat  properties  of  these  spinels  at 
temperatures  well  removed  from  T^: 

1.  The  lattice  modes  in  CdCr-TO.  and  ZnCroO.  are  well 

2  4  2  4 

characterized  by  Debye  temperatures  of  420  and  about  450 
K,  respectively,  in  excellent  agreement  with  the 
predictions  of  the  Lindemann  equation. 

2.  The  excellent  low-temperature  Schottky  fits  leave  little 
doubt  regarding  the  presence  of  anti  ferromagnetic  spin 
waves  . 

3.  From  the  high- temperature  Schottky  fits,  there  is  good 
agreement  for  the  density  of  contributing  spins  in  the 
two  forms  of  CdCr20^  (and  ZnCr20^). 

Attempts  have  been  made  here  to  estimate  the  ratios  of  free 

spin  densities  below  T^^  and  of  antiferromagnetic  spin  waves 

between  the  two  forms  of  CdCr-,0.  and  of  ZnCr_  0 .  .  However,  these 

2  4  2  4 

estimates  depend  critically  on  the  assumptions  that  5,  g^,  g^,  J, 
etc.  are  the  same  in  the  two  forms,  whereas  in  fact  the  stress 
conditions  most  probably  are  not  the  same. 

C .  Magnetocaloric  Measurements 

Measurements  of  magnetocaloric  effects  under  adiabatic 


•-  n  n  — 


conditions  were  performed  over  broad  temperature  ranges  (2-20  K) 
and  magnetic  field  ranges  (<  15  T)  on  samples  of  CCN(9/1), 
ZCN(9/1),  CdCr20^,  and  ZnCr20^.  These  were  the  same  samples  used 
for  the  specific  heat  measurements  above. 

Magnetocaloric  measurements  are  a  sensitive  test  of  the 
contributing  spin  systems  for  the  s imple  reason  that  uncorrelated 
spins  give  rise  to  magnetization  h‘~‘at i nq ,  anti  ferromagnetically 
correlated  spins  give  rise  to  magnetization  cooling,  and 
ferromagnetically  correlated  spins  give  rise  to  hysteret ic 
effects . 

There  are  several  experimental  details  involved  in  these 
measurements  (e.g.,  thermal  time  constants,  addenda  corrections, 
magnetoresistive  effects  in  thermometry,  etc.).  These  details 
have  been  adequately  discussed  previously  (Eckels  ^  ^.  ,  1985) 
and  in  Interim  Reports  and  will  not  be  discussed  here. 

In  addition,  a  voluminous  amount  of  magnetocaloric  data  were 
measured  here,  as  reported  in  Interim  Reports,  and  rather  than 
present  all  these  data  in  this  Annual  Report,  we  shall  present 
selected  and  summary  data. 

Data  Measured  on  CCN(9/1)  and  CdCr20^»  Considering  the  Cd 
spinel  materials  first,  the  magnetocaloric  AT^  data  for  CCN(9/1) 
measured  above  and  below  T^  are  shown  in  Figs.  4-18  and  4-19, 
respectively.  These  are  the  reversible  AT-components ;  the 
hysteretic  components  while  small  (<  5%)  have  been  separated  out. 

In  a  narrow  temperature  range  around  T^^  (i.e.,  7.5  <  T  <  9), 
CCN(9/1)  displays  magnetization  cooling;  at  all  other  tempera¬ 
tures  CCN(9/1)  displays  magnetization  heating.  We  note  in 
particular  in  Fig,  4-18  that  a  AT  “  H  relatio  is  followed 
(inset),  and  at  the  lowest  temperature  AT/T  -  100%. 

Data  at  the  lowest  temperatures  for  CdCr204  are  shown  in 
Fig.  4-20,  and  here,  in  contrast  to  CCN(9/1),  a  significant 
hysteretic  component  is  found  which  indicates  a  presence  of 
ferromagnetically  correlated  spins.  In  addition,  the  data  in 
Fig.  4-20  show  a  dependence  on  the  ramp  ra te ;  that  is,  40  s  ramp 
rates  were  employed  except  Cor  the  4-s  rates  indicated  in  Fig. 


L 1  on 


4-20. 

An  unexpected  discovery  in  CdCr20^  is  shown  in  Fig.  4-21  -- 
namely,  a  field-induced  change  from  magnetization  heating  to 
magnetization  cooling.  The  hysteretic  components  seen  in  Fig. 
4-20  are  negligible  at  the  temperatures  in  Fig.  4-21,  and  these 
latter  data  are  not  sensitive  to  the  ramp  rate.  Note  in  this 
regard  that  the  CCN(9/1)  data  in  Fig.  4-18  and  4-19  do  not  show 
even  a  hint  of  this  crossover  behavior. 

A  crossover  magnetic  field  for  CdCr20^  can  be  seen  in  Fig. 
4-21  (i.e.,  that  field  where  AT^/AH  changes  sign)  ,  and  the 
dependence  of  this  crossover  field  on  temperature  is  shown  in 
Fig.  4-22.  As  seen,  the  minimum  crossover  field  occurs  at  about 
4.8  K. 

Magnetocaloric  data  for  CdCr20^  at  temperatures  where  only 
magnetization  heating  is  seen  Are  shown  in  Fig.  4-23,  and  here 
also  there  is  no  dependence  on  the  ramp  rate. 

Finally,  in  Fig.  4-24  is  shown  a  f ie Id- temperature  phase 
diagram  for  reversible  magnetocaloric  effects  in  CdCr20^,  and  in 
Fig.  4-25  are  shown  combined  phase  diagrams  for  CCN(9/1)  and 

CdCr tO . . 

2  4 

The  phase  diagrams  in  Fig.  4-25  are  very  similar  and  appear 
offset  in  temperature  by  the  difference  in  the  T^^'s  (indicated). 
However,  at  temperatures  well  below  T^^  the  phase  diagrams 
diverge,  indicating  a  much  larger  density  of  free  spins  in 
CCN(9/1)  compared  to  CdCr20^. 

Data  Measured  on  ZCN(9/1)  and  ZnCr^O^.  Following  the  same 
organization  of  the  data  as  for  the  Cd  spinel  materials  above, 
the  magnetocaloric  data  for  ZCN(9/1)  abvoe  and  below  T^^  are  shown 
in  Figs.  4-26  and  4-27,  respectively.  We  note  that  ZCN(9/1) 
displays  nonhysteretic  magnetization  heating  at  all  temperatures, 
and  at  the  lowest  temperature  in  Fig.  4-26  the  magnetocaloric 
effect  is  not  as  large  as  in  CCN(9/1),  Fig.  4-18. 

Complementary  data  measured  on  ZnCt^O^  are  shown  in  Figs. 
4-28  and  4-29,  respectively.  As  with  the  ZCN(9/1)  data  in  Figs. 
4-26  and  4-27,  nonhysteretic  magnetization  heating  occurs  at  all 
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temperatures,  independent  of  the  ramp  rate. 

The  temperature-magnetic  field  phase  diagram  for  ZnCr20^  is 
shown  in  Fig.  4-30,  and  is  indicated.  Combined  magnetocaloric 
phase  diagrams  for  ZCN(9/1)  and  ZnCr20^  are  shown  in  Fig.  4-31. 

The  combined  data  in  Fig.  4-31  for  the  Zn  spinels  are  quite 
dissimilar  compared  to  the  combined  data  for  the  Cd  spinels  in 
Fig.  4-25.  The  data  in  Fig.  4-31  do  not  scale  above  T^^  based  on 
the  difference  in  '  s ,  as  the  data  in  Fig.  4-25  do.  However,  as 
with  CCN(9/1),  the  ZCN(9/1)  material  has  a  much  larger  density  of 
free  spins  at  the  lowest  temperatures  compared  to  the  pure 
spinel. 

"Accidental"  Maqnetocaloric  Event.  During  the  course  of 
magnetocaloric  measurements  on  CCN(9/1)  and  ZCN(9/1),  an 
inadvertent  col  lapse  of  the  magnetic  field  occurred  during 
up-ramp,  and  the  measured  behavior  of  the  samples  during  this 
event  may  provide  important  clues  regarding  the  spin  systems  in 
these  materials. 

The  temperature- time  and  magnetic  field-time  records  for 
this  event  are  shown  in  Fig.  4-32. 

The  Fig.  4-32  data  reveal  two  interesting  features:  (1)  The 
sum  result  of  the  event  was  to  heat  the  samples  far  above  the 
reservoir  temperature,  CCN(9/1)  in  particular;  (2)  Following  the 
field  collapse,  an  instantaneous  cooling  occurred  followed  by  a 
heating  process  with  a  characteristic  time  ~  1-2  sec. 

These  results  might  be  understood  in  terms  of  the  the  model 
of  two  spin  systems:  The  uncorrelated  spins  instantaneously 
demagnetization-cool  with  a  very  rapid  spin-phonon  relaxation 
time;  and  the  correlated  spins  demagnet i zat ion-heat  with  a  long 
spin-phonon  relaxation  time.  The  explanation  then  is  that  on 
field  collapse  the  uncorrelated  spins  first  cooled  the  lattice 
phonons  to  a  very  low  temperature,  followed  by  the  correlated 
spins  slowly  heating  the  phonon  field  back  to  the  reservoir 
temperature  or  higher. 

The  question,  of  course,  is  what  happened  to  the  entropy  in 
the  uncorrelated  spin  system  since  clearly  when  the  field  is 


slowly  ramped  down  the  net  etlect  is  demagnetization  cooling  back 
to  the  reservoir  temperature.  The  answer  here  may  lie  in  the 
zero-field  splitting  of  the  uncorrelated  spins  which  ultimately 
limits  how  low  in  temperature  this  spin  system  can  cool  and 
therefore  limits  the  recoverable  (i.e.,  reversible)  entropy. 

Under  slow-ramp  conditions  the  spin  systems  are  presumably  in 
equilibrium  with  each  other  and  the  phonon  field  so  that  zero- 
field-splitting  effects  aren't  seen. 

In  any  case,  this  accidental  event  seems  to  lend  additional 
weight  to  the  notion  of  two  spin  systems  in  these  spinels.  Also, 
it's  interesting  to  observe  that  the  effect  of  the  correlated 
spins  in  ZCN(9/1)  in  Fig.  4-32  is  weaker  than  in  CCN(9/1),  and 
this  reflects  the  behavior  seen  above  (Figs.  4-18  and  4-26). 

D.  Magnetothermal  Data  Analyses 

In  the  above  sections,  a  large  amount  of  specific  heat  and 
magnetothermal  data  have  been  reported  on  materials  composed  of 
the  spinels  CdCr20^  and  ZnCr20^.  We  now  take  up  analyses  of 
these  data  based  on  the  TdS  equation  for  magnetic  insulators. 
Additionally,  we  will  give  estimates  of  the  densities  of  free 
spins  in  these  materials.  We  begin  wi*-h  the  TdS  equation  for 
magnetic  insulators, 

TdS  =  mC^dT  +  u  mT(3M/3T)^dH  (4-6) 

n  O  n 

where  m  is  the  sample  mass  and  C.,  is  the  specific  heat  at  field. 
Next,  guided  by  the  linear  M  -  H  characteristic  in  Fig.  2-4,  we 
wr  i  te 


M  =  xH, 


(4-7  ) 


and  substituting  in  Eq.(4-6)  for  adiabatic  conditions,  we  have 


(C^/T)dT  =  u^(d  )/dT)dH^. 


(4-8  ) 


The  specific  heat  'aeasuremen ts  on  these  material,  in  intense 


magnetic  fields  have  not  revealed  a  strong  field  dependence. 
Approximately  ,  we  have  upon  integration  that 


u^(dx/dT)  =  -2ASyH- 


(4-:/ ) 


At  the  lowest  temperatures. 


where  AS^  is  the  change  in  the  zero-field  entropy  on  adiabatical- 
ly  ramping  0  +  H.  Thus,  by  this  convention,  if  magne t i za t i on 
heating  occurs  fAS^  >  0),  then  dx/dT  <  0. 

It  will  prove  useful  to  review  the  elementary  dependence  of 
dx/dT  on  x(T),  and  this  is  shown  in  Fig.  4-33.  T^^  can  be  defined 

by  the  maximum  in  dx/dT,  as  shown.  At  the  lowest  temperatures, 
the  presence  of  f ree  spins 

causes  x  to  rise  rapidly  (dashed  L 

\  /  \ 

curve),  and  this  in  turn  will  i  y  \ 

cause  dx/dT  to  plunge  rapidly,  as  ^  \  /{ 

shown.  \ 


In  arriving  at  Eg. (4-9),  the 
most  serious  approximation  is  that 
of  a  1  inear  paramagnet,  Eg. (4-7), 
and  it  is  instructive  to  consider 
the  corrections  to  this  approxima¬ 
tion.  In  the  most  general  case, 
we  have  the  Landau-type  expansion. 


M  =  xH  +  C  ^  -t 


(4-10) 


Substituting  in  Eg. (4-8)  and  re¬ 
arranging,  we  have 


Figure  4-33.  Sketch  of 
the  variation  in  dx/dT 
due  to  the  idealized 
change  of  x  with  tempera¬ 
ture. 


UQ(dx/dT)  =  -2ASyH^  +  u^(|(dC/dT)H^  +  |(dt/dT)H‘^ 


(4-11 


We  note  two  features  from  Eg. (4-11):  First,  even  in  the  case  of 
a  nonlinear  paramagnet,  if  ^  and  t  are  tempera tu re- i ndependen t 
then  Eg. (4-11)  reduces  to  Eg. (4-9).  And  second,  under  isothermal 
conditions  (i.e.,  dx/dT  is  constant),  if  AS  /h"  is  not  constant 


then  the  importance  of  the  higher  terms  in  Eg. (4-11)  is  clearly 
indicated.  We  shall  see  below  that  these  considerations  become 
important  in  interpreting  the  phonomenological  results. 

Finally,  anticipating  that  AS^  tables  will  be  generated, 
another  interesting  analysis  of  the  data  becomes  possiole. 

Namely,  tht  maximum  entropy  available  in  a  free-spin  system  is 
simply 

AS  -  n  R  ^n  (2s  +  1 )  (4-12  ) 

m  p 

where  n  are  the  number  of  free  spins  per  formula  weight  (e.g. , 

P 

n^  =  2  in  tfiO  h  i  gh- tempera  ture ,  paramagnetic  phase  of  CdCr^O^) 
and  s  is  the  spin  value.  If  we  now  assume  that  at  some 
sufficiently  large  field  strength  all  of  these  free  spins  are 
aligned,  then  upon  demagnetization  the  temperature  change  of  the 
sample  is  due  to  the  randomization  of  these  spins.  Stated 
di  f ferently , 

AS  =  n  R  iln  ( 2s  +  1  )  (4-13 ) 

op' 

where  As^  is  the  entropy  change  of  the  sample,  as  above.  In  this 
fashion,  the  free-spin  density,  n^,  can  be  estimated  from  the 
very  high  field  magnetocaloric  data. 

The  first  step  in  these  analyses  was  to  develop  AS^-tables 
from  the  four  sets  of  zero-field  specific  heat  data  in  Figs.  4-9 
and  4-10.  This  was  done  by  numerically  integrating  the  experi¬ 
mental  Cq/T  data,  and  then  these  data  were  smoothed  for  greater 
accuracy.  Next,  a  code  was  written  for  converting  a  "*^2 

netocaloric  event  on  0  -*■  H  up  ramp  to  the  ratio  AS^/H^  according 
to  Eg. (4-9).  The  average  temperature  of  the  event  is  (T^  + 

and  only  the  reversible  AT ' s  were  considered.  In  this 
fashion,  two  types  of  phase  diagrams  were  generated  for  each  of 
the  four  samples:  (1)  dx/dT  T  at  constant  H;  and  (2) 

dx/dT  v_s.  H  at  constant  T.  In  the  former  type  of  diagram,  the 
data  are  limited  to  H  =  5  and  10  T  because  complete  data  sets 
were  unavailable  at  other  fiuid  levels,  in  particular  at  lowe r 


field  levels. 


CdCL'^O^  Disk.  The  magnetotheemal  phase  (iiagrams  fes  the 
compacted  CdCi  ^O^  disk  are  shown  in  Fig.  4-34.  The  cJx/dT  vs . 

T  diagram  in  Fig.  4-34(a)  at  5  and  10  T  s'-aw  that  f  i  e  1  d-dependen  t 
magnetization  coo  ling,  dx/dT  >  0,  is  observed  at  and  below  T^.,  as 
discussed  previously. 

We  note  here  that  for  an  ideal  Curie  magnet  x  “  4’  ^  /  so  in 
some  temperature  range  dx/dT  «  -T  ,  as  illustrated  in  the 
idealized  sketch  in  Fig.  4-33.  The  Fig.  4-34  (a)  data  rjo  not 
reflect  this  up  to  23  K ;  otherwise,  these  Fig.  4-34(a)  data 
follow  the  idealized  behavior  in  Fig.  4-33  and  indicate  that  a 
free-spin  contribution  exists  below  about  3  K. 

The  dx/dT  v^ .  H  magnetothermal  phase  diagram  in  Fig, 

4-34  (b)  is  shown  toi  temperatures  below  and  above  T.^;  at  both 
temperatures  magnetization  heating  is  observed.  For  both 
temperatures,  d;</dT  decreases  rapidly  as  H  decreases,  so  little 
quantitative  agreement  can  be  expected  with  the  low-field  x”Ci'3t‘-'3 
in  Section  V  below. 

The  dx/<J'l“data  at  3,47  K  in  Fig.  4-34(b)  involve  contribu¬ 
tions  from  both  correlated  and  uncorrelated  spins.  However,  at 
10.6  K  one  expects  all  the  spins  to  be  uncorrelated,  and  hero  the 
non--cons  ta  ncy  of  Kg.  (4-9)  is  clearly  seen.  The  higher-order 
terms  in  l.q.M-ll)  are  needed  to  explain  these  data,  but  this  was 
not  pursued  due  to  the  lack  of  data  at  lower  H-valuos  in  Fig. 
4-34(b).  r.’-dc*  in  this  regard  that  the  measurement  of  magneto¬ 

caloric  data  It  1  jwer  H-values  is  subject  to  the  experimental 
constraint  of  restjlving  the  associated  small  temperature  changes 
in  intense  m a g n e t i c  fields. 

Finally,  from  Fig.  4-34(b)  it  appears  that  saturation  occurs 
by  about  15  T. 

CCN ( 9 / 1 ) .  The  dx/dT  v^.  T  magnetothermal  phase  diagram 
for  CCN{9/1)  IS  shown  in  Fig.  4-35(a),  and  this  diagram  closely 
resembles  that  of  CdCr29^  in  Fig.  4-34(al  given  the  shift  in  the 
T^-values.  Note  also  that  the  scales  (i.e.,  magnitudes)  of  the 
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data  in  Figs.  4-34(a)  ^  4-34(d)  ai'e  commensurate. 

The  dx/dT  H  diagram  for  CCN(9/1)  is  shown  in  Fig. 

4-35(b)  at  two  temperatures  above  and  below  T^,  and  once  again 
the  non-constancy  of  Eq.(4-9)  is  clearly  evident. 

The  larger  density  of  free  spins  below  T^^  in  CCN(9/1) 
compared  to  the  CdCr^O^^  disk  is  clearly  seen  in  both  Figs. 

4-35(a)  and  4-35(b)  by  the  much  larger  absolute  values  of  dx/dT 
(see  be  low ,  also). 

ZnCr ^0^  Disk.  The  dx/dT  vs.  T  magne to  the rma 1  phase 
diagram  for  ZnCr.,0,  in  Fia.  4-36{a)  is  in  marked  contrast  to  the 
equivalent  phase  diagram  for  CdCr20^  in  Fig.  4-34(a)  and  also 
does  not  resemble  the  idealized  behavior  in  Fig.  4-33.  In  Fig. 
4-36(a)  dx/dT  goes  through  a  deep  minimum  at  T^  at  both  5  and  10 
T,  indicative  of  a  near-singularity  in  x  y_s.  T.  At  all  tempera¬ 
tures  in  Fig.  4-36(a),  dx/dT  <  0,  indicatin.g  that  x  always  in¬ 
creases  with  decreasing  temperature  at  these  field  strengths. 

The  dx/dT  H  magne tothermal  phase  diagram  for  the 

ZnCr20^  disk  is  shown  in  Fig,  4-36(b),  again  at  two  temperatures 
above  and  below  T.^,.  As  in  CdCr20^  [Fig.  4-34  (b)],  dx/dT  de¬ 
creases  very  rapidly  with  decreasing  H  indicative  of  the  impor¬ 
tance  of  the  h i ghe L-order  terms  in  Eg. (4-11).  Also,  the  appear¬ 
ance  of  saturation  is  seen  at  about  15  T. 

AnotheL'  striking  difference  between  CdCr20^  and  ZnCr20^  is 
seen  on  comparing  the  equivalent  phase  diagrams  in  Figs.  4-34(b) 
and  4-36(b).  Namely,  at  T  ~  2  T^^ ,  the  dx/dT-values  for  Zncr20^ 
are  substantially  larger  in  magnitude  compared  to  CdCr20_^. 

ZCN ( 9/1 ) .  Ihe  dx/dT  vs.  T  magnetothermal  phase  diagram 
for  ZCN(9/1)  in  Fig.  4-37(a)  closely  resembles  that  for  the 
ZnCr20^  disk  in  Fig.  4-36{a)  --  note  in  particular  the  deep  mini¬ 
mum  at  T,,  and  the  similarity  of  the  scales. 

The  dx/dT  .  H  phase  magnetothermal  diagram  for  ZCN(9/1) 
is  given  in  Fig.  4-37(b)  for  three  temperatures  —  above,  at,  and 
below  T  .  Again  the  similarity  to  the  equivalent  diagram  for 
ZnCr^O.  [Fig.  4-36 (b)l  is  observed,  and  we  note  the  huge 


ZCN(9/1) 


t  i>t  liermii  1  plKi.se  tlingr.ims  for  ZCN(9/1)  according,  to 
i(|.('t-9)  as  fiuiot-ioiis  of:  (a)  Toiii()era  t  u  re ,  and  (h) 


2/2 


(absolute)  values  for  d^/dT  in  ZCN(9/1)  at  T  ~  in  Fig. 
4-37(b) . 


Free  Spin  Densities.  As  a  final  exercise,  we  consider  the 
free-spin  density  n^  computed  from  Eq.(4-13).  The  magnetocaloric 
data  used  here  to  obtain  AS  are  taken  at  the  highest  field 
strength,  15  T,  and  the  y  ^  dy/dT  y_s.  H  phase  diagrams  suggest 
that  this  field  strength  is  very  near  the  saturation  level  at 
certain  temperatures.  The  free-ion  spin  value,  s  =  3/2,  was  used 
here.  These  computed  n^  vs.  T  data  for  the  four  samples  are 
shown  in  Fig.  4-38.  Only  one  point  each  for  the  spinel  + 
columbite  ceramic  is  shown  in  Fig.  4-38  (i.e.,  at  15  T). 


The  n^  data  in  Fig.  4-38  at  temperature  above  T^^  are  sur¬ 
prisingly  consistent  for  CdCr20^  and  ZnCr20^,  n^  =  0.04,  but  far 

below  the  expected  value  of  n  =2.  To  bring  these  data  into 
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line  with  n  =2  requires  s  ~  6x10 

P 

The  data  for  n^  at  temperatures  well  below  T^^  are  also 
surprisingly  consistent  for  all  four  samples  in  Fig.  4-38,  n^  = 
0.003-0.007.  The  data  for  CdCr20^  in  Fig.  4-38  do  not  extend 
below  7  K  because  it  will  be  recalled  that  there  is  a  field- 


enforced  crossover  in  the  sign  of  dx/dT  for  this  material  at  the 
lower  temperatures. 


We  point  out  that  Eq.(4-13)  is  not  strictly  valid  at  temper¬ 


atures  below  Tj^.  Namely,  in  the  simple  two-spin  model  the 

uncorrelated  free  spins  contribute  a  dx/dT  <  0  component,  the 

correlated  spins  contribute  a  dx/dT  >  0  component.  Consequently, 

n  below  T„  is  a  measure  of  the  net  uncorrelated  spins, 
p  N  -  ^ 


E.  Conclusions  from  Thermal  Data 

A  large  amount  of  experimental  data  have  been  measured  on 
the  same  compacted  disks  of  CdCt^O^  and  ZnCi^O^  --  zero-field 
specific  heat  data,  field-dependent  specific  heat  data,  and 
magnetocaloric  data.  The  magnetocaloric  data  have  been  reduced 
and  presented  in  the  form  of  phase  diagrams.  The  zero-field 
specific  heat  data  have  been  analyzed  in  terms  of  Schottky  and 
spinwave  models.  Finally,  all  the  magnetocaloric  data  have  been 


analyzed  according  to  the  TdS  equation,  making  use  of  the  speci¬ 
fic  heat  data.  In  what  follows  below,  the  results  of  all  these 
studies  are  summarized: 

1.  In  contrast  to  ZnCr20^,  CCN(9/1),  and  ZCN(9/1),  the  CdCr^O^ 
disk : 

(a)  Does  not  display  a  sharp  maximum  in  the  specific  heat; 

(b)  Displays  a  field-enforced  crossover  between  dx/dt  >  0 

and  d  x/dT  <  0  below  and 

(c)  Exhibits  a  substantial  hysteretic  component  below  in 
the  magnetocaloric  data. 

2.  Except  at  the  lowest  temperatures,  there  is  a  great  deal  of 
similarity  in  the  AT^  phase  diagrams  of  CdCr20^  and 
CCN(9/1),  Fig.  4-25.  The  corresponding  diagrams  for  the 
zinc  spinels  are  not  similar.  Fig.  4-31. 

3.  Both  CdCr20^  and  CCN(9/1)  exhibit  magnetization  cooling  in  a 

limited  temperature  range  at  and  below  Fig  4-25;  in  con¬ 

trast,  ZnCr20^  and  ZCN(9/1)  exhibit  magnetization  heating  at 
all  temperatures.  Fig.  4-31. 

4.  Both  CdCr20^  and  ZnCr20^  show  strong,  nonlinear  departure 
from  the  linear  relation  Eq.(4-9),  as  seen  in  Figs.  4-34(b) 
and  4-36(b).  This  was  unexpected,  given  the  linear  M  H 
data  in  Fig.  2-4,  and  frustrates  meaningful  comparisons  with 
the  low-field  susceptibility  data  in  Section  V  below. 

5.  Specific  heat  data  for  CdCr^O.  and  CCN(9/1)  are  both  consis- 

3  ^  ^  . 

tent  with  a  T  anti  ferromagnetic  spin-wave  contribution  at 
the  lowest  temperatures  and  a  =  420  K  acoustic-mode  con¬ 
tribution  . 

6.  Specific  heat  data  for  ZnCr20^  and  ZCN(9/1)  are  both  consis¬ 
tent  with  a  antiferromagnetic  spin-wave  contribution  at 
the  lowest  temperatures  and  a  0^^  ~  450  K  acoustic  mode  con¬ 
tribution  . 

7.  For  both  the  cadmium  and  zinc  spinels,  the  Debye  tempera¬ 
tures  are  in  excellent  agreement  with  predictions  from  the 
Lindemann  equation. 

8.  From  the  T  >  Schottky  fits  to  the  specific  heat  data  of 
all  four  samples,  there  is  satisfactory  agreement  for  the 


number  of  spins  partaking  in  the  transition  at  between 
CdCt^O^  and  CCN(9/1)  and  between  ZnCr20^  and  ZCN(9/1). 

From  the  T  <  Schottky  fits  to  the  specific  heat  data,  the 
densities  of  free  spins  in  the  disks  of  CdCr20^  and  ZnCt^O^ 
are  ~  5-7  times  larger  than  in  the  densified  ceramics 
CCN(9/1)  and  ZCN(9/1).  This  is  in  semi-quantitative 
agreement  with  the  magnetocaloric  data  here  based  on  entropy 
arguments.  Fig.  4-38.  On  the  other  hand,  the  low-field  dc 
susceptibility  data  below  suggest  that  the  ceramics  CCN(9/1) 
and  ZCN(9/1)  have  the  larger  densities  of  free  spins. 

From  the  spin-wave  terms  in  the  specific  heats  below  the 

densities  of  the  antiferromagnetically  correlated  spins  are 
much  larger  in  the  disks  of  CdCr20^  and  ZnCr20^  than  in  the 
ceramics  of  CCN(9/1)  and  ZCN(9/1),  respectively.  Nonethe¬ 
less,  the  spin-wave  terms  do  indicate  that  magnetization 
cooling  should  be  more  pronounced  in  the  cadmium  spinel  than 
in  the  zinc  spinel,  in  agreement  with  the  magnetocaloric 
data  [We  draw  attention  to  the  assumption  made  in  analyzing 
the  specific  heat  that  the  zero-field  splitting,  etc.  is  the 
same  in  both  forms  of  CdCr20^  (and  ZnCr20^),  whereas  in  fact 
different  stress  conditions  exist  in  the  two  forms). 

In  the  paramagnetic  region,  T  >  for  both  CdCr20^  and 

ZnCr20^  and  for  magnetic  field  strengths  at  or  very  near 
saturation  (15  T),  the  densities  of  free  spins  estimated 
from  entropy  arguments  are  ~  0.04  —  i.e.,  ~  50  times 
smaller  than  anticipated.  Fig.  4-38. 

In  all  the  thermodynamic  analyses  here,  the  field  dependence 
of  the  specific  heat  has  been  ignored  (i.e.,  C„  =  C  ).  This 

rl  O 

is  not  a  serious  assumption  because  C„/C  <  35%  near  T.,  at 

^  H  O  ~  N 

10  T,  Fig.  4-13. 


V.  MAGNETIC  SUSCEPTIBILITY  AND  RESONANCE  STUDIES 

Magnetic  Susceptibility  (x)  and  EPR  Studies 

This  section  summarizes  the  magnetic  susceptibility  ( x)  and  EPR 
measurements  on  several  CdCr204-  and  2nCr204-based  spinels,  made  (by  the  WVU 
group)  with  a  view  of  providing  new  clues  as  to  why  these  spinels  exhibit 
exceptionally  high  heat  capacities. 


V . 1  Experimental  Methodology,  Sample  Preparation  and  Sample  Notation 

The  X  measurements  were  made  at  the  National  Magnet  Laboratory  of  MIT, 
using  a  SHE  squid  magnetometer.  The  EPR  measurements,  on  the  other  hand,  were 
made  at  West  Virginia  University  with  an  IBM-Bruker  Model  ER-200D  spectrometer 
which  uses  an  ASPECT  2000  microcomputer  for  data  acquisition  and  analysis. 

The  description  of  the  samples,  investigated  in  this  report  is  given  in 


Table  5.1  al'.'ough  measurements  were  made  on  many  others  as  well. 


Table  5.1  Investigated  Samples 


Sample  Name 

Notation 

Source 

CdCr204  (powder) 

CdCr204  (P) 

Penn  State  (PS)  Univ. 

CdCr204  (powder  disk)^ 

CdCr204  (PD) 

PS  Univ. 

CdCr204  (densified  dlsk)^ 

CdCr204  (DDCP) 

Ceram  Physics  (CP) 

CdCr204  (densified  disk)*^ 

CdCr204  (DPPS) 

PS  Univ. 

ZnCr204  (powder) 

ZnCr204  (P) 

PS  Univ. 

ZnCr204  (powder  disk)^ 

ZnCr204  (PD) 

PS  Univ. 

ZnCr204  (densified  dlsk)^ 

ZnCr204  (DDPS) 

PS  Univ. 

ZnCr204  (densified  disk)^ 

ZnCr204  (DDCP) 

Ceram  Physics 

a.  the  powder  compressed  into  a  disk  at  1290“C/1  hr  for  CdCr204  (p  =  4.4  g/cc) 

and  at  1650“C/1  hr  for  ZnCr204  (p  =  3.8  g/cc)  (Section  III). 

b.  contains  columbite 

c. -  CdCr204  +  5%  Nb205  heated  at  1320®C/1  hr  (p  =»  5.3  g/cc)  (Section  III). 

d.  ZnCr204  +  5%  Nb205  heated  at  1350°C/1  hr  (p  «  4.87  g/cc)  (Section  III). 


V.2.  Experimental  Results  on  Susceptibility  (x) 


V.2.a. 


vs .  T 


Figure  5.1  shows  the  temperature  dependence  of  x  the  four  samples 

of  CdCr204,  namely  powder  (P) ,  powder  disk  (PD)  and  the  densified  disks  (DDCP 
and  DDPS).  The  high  temperature  (T  >  10  K)  behavior  of  P  and  PD  is  very  nearly 
the  same,  as  is  also  the  case  for  DDCP  and  DDPS.  However,  the  low  temperature 
behavior  (shown  in  detail  in  the  inset),  reflecting  the  cooperative,  magnetic 
ordering,  is  different  suggesting  that  sample  preparation  conditions  do  affect 
the  microscopic  spin  arrangement  in  every  sample. 


TCK] 


Figure  5.1  Plots  of  susceptibility  (x)  vs.  temperature  (T)  for  the  four  samples 
of  CdCr204.  The  inset  shows  the  low-temperature  x  vs.  T  behavior. 

Figure  5.2  shows  the  x  vs.  T  behavior  for  the  various  samples  of  ZnCr204. 

It  is  clear  that  both  the  T  <  and  the  T  >  Tjj  parts  of  the  x  vs.  T  plots  are 
different  for  different  samples.  All  the  samples  except  PD  show  a  fairly  well 
defined  peak  at  or  Immediately  above  T»j.  The  PD  sample  exhibits  a  very  broad 
peak,  quite  untypical  of  an  antlferromagnetlc-paramagnetlc  phase  transition. 
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Figure  5.3.  x  versus  T  behavior  for  CdCr20it  (PD)  in  magnetic  fields,  of  0.2 
and  20  kG.  The  behavior  of  CdCr20it  (P)  at  0.2  kG  is  also  shown, 


The  changes  in  the  x  versus  T  behavior  of  the  ZnCr204  pvwder  are  quite 
dramatic  as  shown  by  the  results  in  Figure  5.4  at  the  fields  0.2  and  10  kG.  The 
features  which  show  up  markedly  at  0.2  kG  get  suppressed  considerably  at  10  kG. 


Figure  5.4.  x  versus  T  behavior  for  ZnCr204  (P)  in  two  fields  of  0.2  and  10  kG 


V.2.C  X  versus  T  in  zero  field  cooling 

Since  the  x  vs.  T  behavior  was  fo’ind  to  depend  on  the  field,  and  in  all 
the  measurements  the  samples  were  field-cooled,  it  was  decided  to  cool  the 
samples  in  zero  field  for  x  measurements.  Results  of  this  experiment  (Figure 
5.5)  on  ZnCrjOv  (PD)  from  6  K  to  20  K,  exhibit  the  same  behavior  as  that  of 
Figure  5.2  except  that  the  x  values  for  the  zero-field-cooled  samples  are  at 
least  25Z  larger  in  magnitude.  The  first  derivative  of  x  w.r.t.  T,  which  gives 
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the  Neel  temperature,  is  however  at  the  same  temperature,  9.4  ±0.1  K  (within 
experimental  accuracy),  as  in  the  field-cooled  case. 


Figure  5.5«.  x  versus  T  behavior  for  ZnCr204  (PD)  when  the  sample  was  cooled  in 
zero  field.  The  curve  peaking  at  9.4  ±  0.1  K  shows  the  temperature 
dependence  of  dx/dT. 


V.2.d  Magnetic  field  (H)  dependence  of  magnetization  (M) 

To  find  if  these  systems  exhibit  any  spin-flop  transition,  the 


magnetization  (M)  was  studied  as  a  function  .f  H,  up  to  45  kG.  Figure  5.6  shows 
the  M  vs.  H  plot  for  a  PD  sample  of  CdCr204,  the  same  sample  as  for  Figure  5.3. 
For  this  sample,  a  linear  dependence  of  M  on  H  is  seen  at  10  K  (T  >  Tjj)  with  no 
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measurable  intercept  on  the  H  =  0  axis.  At  T  =  1.91  K,  however,  a  small  but 
finite  intercept  was  detected.  This  small  value  of  the  Intercept  was  found  to 
be  sample  dependent,  since  another  sample  examined  did  not  yield  a  measurable 


Intercept,  as  may  be  seen  from  Figure  5.7. 
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Figure  5.6.  Typical  M  vs.  H  behavior  for  a  powdered  disk  (PD)  .sample  of  CdCr20i 


Similar  results  were  obtained  for  a  powdered  disk  (PD)  sample  of  ZnCr204. 
The  results  are  shown  below,  Figure  5.7. 
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Figure  5.7.  M  vs.  H  behavior  of  a  powdered  disk  (PD)  sample  of  ZnCr204. 


V.3  EXPERIMENTAL  EPR  RESULTS 


As  the  complement  to  the  magnetic  susceptibility  data,  EPR  measurements 
were  made  above  liquid  nitrogen  temperatures  for  all  the  samples  listed  in  Table 
5.1.  For  some  samples,  the  measurements  were  made  both  in  the  X-band  (~9.5  x 
10^  Hz)  and  Q-band  (33.5  x  10^  Hz)  regions  to  measure  the  spin  dynamics. 

The  lineshape  analysis  of  all  the  investigated  samples  shows  the  lines 
nearer  to  a  Lorentzian  than  a  Gaussian,  over  the  temperature  range  80-400  K. 

The  lines  are  mainly  exchange  narrowed. 

The  EPR  dependence  of  the  peak-to-peak  linewidth  ( ABpp)  for  the  powder 
samples  of  CdCr204  and  ZnCr204  is  shown  in  Fig.  5.8.  This  plot  corresponds  to 
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the  Q-band  microwave  frequencies  and  the  values  for  ABpp  agree  with  those 
obtained  using  X-band  frequencies.  The  data  for  other  samples  (PD,  DDPS,  DDCP) 
is  plotted  in  Fig.  5.9.  The  linewidth  increases  smoothly  with  decrease  in 
temperature  for  all  the  samples  except  CdCr204  (DDCP)  which  showed  a 
sorae-sort-of  discontinuity  around  305  K.  This  was  shown  to  be  due  to  the 
presence  of  some  Cr203  which  has  an  antiferromagnetic  transition  around  308  K. 


Figure  5.8  Dependence  of  the  EPR  peak-to-peak  linewidth  (ABpp) 
temperature  (T)  for  the  spinel  powders. 


as  a  function  of 


Figure  5.9  Dependence  of  the  EPR  linewldth  ABpp  as  a  function  of  temperature 
(T)  for  the  CdCr204  and  ZnCr20i4  samples. 

V.4  ANALYSIS  OF  THE  x-DATA 
V.4.a  Neel  Temperature  (Tfj) 

The  Neel  temperatures  were  obtained  to  within  ±0.1  K,  as  the  temperatures 
corresponding  to  peaks  In  the  (dx/dT)  vs.  T  plots,  using  the  earlier  discussed  x 
vs.  T  plots  as  the  base  data.  Figures  5.10  and  5.11  show  the  (dx/dT)  vs.  T 
plots  for  CdCr204  and  ZnCr204  respectively. 


ZnCr204 


Figure  5.13.  Schematic  display  of  the  Neel  temperatures  for  ZnCr204  samples. 


V.4.b  Susceptibility  behavior  for  T  >  % 


Since  CdCr204  and  ZnCr204  are  known  to  undergo  antif erroraagnetic 
orderings  at  Tjq,  the  temperature  dependence  of  x  T  >  Tfj  is  expected  to  obey 
the  Curie-Welss  law 
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T  >  Tn 
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where  C^f  =  Nu^/Sk  with  N  defined  as  the  number  of  spins,  u  the  effective 
magnetic  moment  of  the  spins  and  k  the  Boltzmann  constant.  A  plot  of  1/x  versus 
T  would  thus  be  a  straight  line  with  its  slope  equal  to  l/Cgf  and  the  Intercept 
equal  to  -9/c.  Experimental  results  of  1/x  versus  T  for  the  CdCr204  samples  are 
shown  in  Figure  5.14.  Linear  behavior  is  Indeed  observed  but  only  for  T  >  5  T^. 
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This  result  clearly  implies  the  existence  of  short-range  order  between  the  spins 
for  Tfj  <  T  <  5  Tfj.  The  values  of  9  and  C^f  thus  obtained  are  listed  in  Table 

5.2. 


Table  5.2.  Values  of  the  parameters  C^f  and  9  obtained  by  fitting  the  T  »  T«^ 
data  of  CdCr204  samples  to  the  expression  x  =  Caf/T-9 


powder 

powder  (*] 

powder  disk  (PD) 

DDPS 

DDCP 


Cgf  [emu'K'g'M 

0  [K] 

1.2  X  10-2 

-83 

- 

-83 

1.14  X  10-2 

-84 

1.10  X  10-2 

-80 

0.81  X  10-2 

-60 

Figure  5.15,  l/x  versus  T  for  ZnCr20H  samples  at  T  > 


The  data  for  the  various  samples  of  ZrCr204  were  also  plotted  as  (l/x) 
versus  T  (Figure  5.15)  and  yielded  the  parameters  Cgf  and  9  as  given  in  Table 
5.3.  An  overwhelming  presence  of  short-range  spin  order  can  clearly  be  seen 
for  all  the  samples. 


Table  5.3.  Values  of  the  parameters  Cgf  and  0  obtained  by  fitting  Che  high 

temperature  data  of  ZnCr204  samples  to  the  expression  x  =  Cgf/T-B 


— - - - - — 

Caf  (emu*K’g“M 

0  [K] 

powder 

1.17  X  10-2 

-331 

powder  [*] 

- 

-392 

powder  disk  (PD) 

1.30  X  10-2 

-396 

DDPS 

• 

X 

o 

1 

-407 

DDCP 

1.26  X  10-2 

-389 

(*]  Baltzer  et  al .  ,  Phys .  Rev.  151  ,  367  (  1966). 
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c  Susceptibility  behavior  below  T(^ 


One  region  of  particular  interest  below  Tjij  is  where  x  increases  with  a 
decrease  in  temperature,  the  Curie-like  behavior.  This  is  believed  to  be  due  to 
the  presence  of  "free"  paramagnetic  Curie-like  spins.  Assuming  that  the 
contribution  to  x  from  the  ant  if erroraagnetlcally-ordered  spins  is  Independent  of 
temperature  in  the  range  '  ^  T  <  0.5  Tjj,  the  experimental  data  can  be  fitted  to 
the  following  expression 

X  =  Xg  +  Cpara/T  T  <  Tn  (V.2) 

where  Xg  =  X  (T  -*•  0)  and  Cpa^a  =  Nu2/3k.  A  plot  of  x  versus  1/T  would  thus 
yield  Cpara  hence  the  number  of  paramagnetic  spins.  The  plots  of  x  versus 
1/T  are  given  in  Figures  5.16  and  5.17  for  CdCr204  and  ZnCr204  samples, 
respectively.  The  values  of  Xg  and  Cpara*  obtained  from  the  slopes  of  these 
plots,  are  collected  in  Tables  5.4  and  5.5. 


DDPS 


7  0.4. 

cr) 
a  . 

E 


^  DDCP 


ZnCr^O^ 


1/T  £K"’] 


Figure  5.17.  x  versus  1/T  for  ZnCr20i,  samples. 


Table  5.4.  Values  of  the  parameters  Xp  and  Cpara  obtained  by  fitting  the  low 

temperature  data  of  CdCr204  samples  to  the  expression  X  =  Xq  +  Cpara 


CdCr 2O4 

Xq  (erau’g'^) 

1  .  -  * 

Cpara  (o'’i'u’K*g 

powder 

1 

1.18  X  lO-** 

1.96  X  10"^ 

PD 

1.28  X  lO"** 

2.14  X  10"® 

DDPS 

1.21  X  10~'* 

2.86  X  10~5 

DDCP 

0.92  X  10~‘* 

9.60  X  10~^ 

Table  5.5.  Values  of  the  parameters  Xp  and  Cpaj-a  obtained  by  fitting  the  low 
temperature  data  of  ZnCr204  samples  to  the  expression  x  =  Xp 


ZnCr204 

Xp  (emu'g'O 

Cpara  (orau*K"g 

powder 

0.38  X  IQ-"* 

1.43  X  10"® 

PD 

0.24  X  10~** 

2.50  X  10"^ 

DDPS 

0.32  X  lO'** 

2.33  X  10"5 

DDCP 

0.22  X  lO-** 

1.92  X  10"5 

V.4.d  "Free"  spin  concentration 


An  important  parameter  <Iodnced  from  the  above  analysis  is  the  fraction 
of  the  "Curie-like"  or  the  "free"  spins,  as  given  by  the  ratio  Cp^  +  C^f . 

Tables  5.6  and  5.7  list  the  fraction  of  the  "free"  spins  in  the  various  samples  of 
CdCr204  and  ZnCr204  respectively.  It  is  clear  that  the  densifled  disks  made  by 
Ceram  Physics  hare  the  maximum  number  of  "free"  spins. 

Table  5.6.  Fraction  of  the  "free"  spins  in  various  samples  of  CdCr204 


^para^*^para  ^af 

PD 

1.9  X  lO-** 

DDPS 

2.6  X  10"3 

DPCP 

1.2  X  10-2 

_ 1 

Table  5.7.  Fraction  of  the  "free"  spins  in  various  samples  of  ZnCr204 


C^0O(V(Cp^+^  Cgf) 

PD 

1.9  X  lO-"* 

DDPS 

1.6  X  10-3 

DDCP 

1.5  X  10-3 

V.5  ANALYSIS  OF  THE  EPR  DATA 


Due  to  experimental  difficulties  in  obtaining  reliable  data  for 
temperatures  below  ~80  K,  we  did  not  attempt  to  analyze  the  EPR  llnewldths 
quantitatively  using  rigorous  expressions.  This  must  await  a  little  longer  when 
we  expect  to  obtain  extensive  lineshape  data  at  liquid  helium  temperatures. 


Nevertheless,  It  is  fairly  clear  fron  the  magnitudes  of  the  peak-to-peak 
linewidths  (ABpp)  for  the  CdCr204  and  ZnCr204  •■amples  (Figures  5.8  and  5.9)  that 
the  exchange  narrowing  is  A-5  times  more  effective  in  ZnCr204  than  CdCr204 
samples.  This  is  essentially  what  is  expected  from  the  susceptibility  analysis 
above  Tf^,  since  the  values  of  the  paramagnetic  temperature  6  for  ZnCr204  and 
CdCr204  samples  are  ~390  and  '80  K  respectively.  The  value  of  J,  the  exchange 
constant,  which  is  directly  proportional  to  9  is  thus  about  5  times  larger  i" 
ZnCr204  samples.  The  EPR  results  thus  provide  a  direct  support  for  the 
essential  correctness  of  the  basic  model  for  the  temperature  dependence  of  x 
all  these  samples. 

V.6  DISCUSSION 

In  this  section,  we  make  a  critical  evaluation  of  the  implications  of  the 
X-data  with  regard  to  the  extensive  results  of  the  specific  heat  and 
raagnetocaloric  experiments.  We  focus  our  attention  on  two  samples,  PD  and  DDCP, 
for  which  a  complete  set  of  measurements  is  available. 

V.6. a  Magnetocaloric  vs.  Magnetic  Susceptibility  Results 

ATp,  the  reversible  part  of  the  raagnetocaloric  change  in  temperature, 
has  been  found  to  go  through  a  narrow  negative  range  (implying  magnetization 
cooling)  for  both  the  PD  and  DDCP  samples  of  CdCr204  (Section  IV).  Whilst  the 
magnitude  of  ATj.  in  this  range  (magnetization  cooling)  is  nearly  the  same  for 
the  PD  and  DDCP  samples,  ATj-  is  much  larger  for  DDCP  than  for  PD  below  Tj; 
(magnetization  heating).  Recall  that  magnetization  heating  results  from  the 


"free”  spins  whilst  rna^net Izat Ion  cooling  from  the  "ordered”  spins.  The 


presence  of  ~50  times  more  "free”  spins  in  DDCP  than  in  PD,  as  shown  by  the 
susceptibility  measurements  (Table  5.3),  is  in  line  with  much  larger  ATj-  for 
DDCP  than  PD.  As  for  the  magnetization  cooling  around  Tfj,  we  can  speculate  that 
the  "ordered”  spins  dominate  the  magnetization  heating  effects  of  the  "free” 
spins  and  give  rise  to  negative  ATj..  Since  the  amount  of  magnetization  cooling 
is  nearly  the  same  in  both  PD  and  DDCP,  one  would  expect  the  same  concentration 
of  the  "orderod"  spins  for  the  two  cases,  provided  the  exchange  coupling 
constant  (J)  among  the  ordered  spins  remains  the  same.  That  t'n.e  exchange 
coupling  constants  are  not  the  same  for  PD  and  DDCP  samples  of  CdCr204  can  be 
seen  from  the  different  values  of  9  (Table  5.2);  note  that  J  is  related  to  6  via 
the  expression  |j|  =  3kg9/2zS(S+l)  where  z  is  the  number  of  nearest  neighbor 
spins.  If  the  dependence  of  ATj.  on  J,  T  and  the  spin  concentration  of  the 
ordered  spins  were  known,  the  ratio  of  the  ordered  spins  at  Tj^  in  PD  and  DDCP 
can  be  calculated.  It  seems,  however,  that  ATj-  depends  rather  strongly  on  T 
since  the  raagnetocaloric  cooling  disappears  rapidly  as  we  go  from  T;»;  to  lower 
temperatures . 

In  the  case  of  ZnCr204,  ATj.  does  not  assume  negative  values  i.e.,  it  shows 
only  magnetization  heating.  This  Implies  that  here  the  free  spins  dominate  the 
ordered  spins.  Susceptibility  results  (Figure  5.2)  especially  for  PD,  show 
clearly  that  around  Tjq  significant  amount  of  short-range  order  (or  lack  of 
long-range  order)  exists  in  these  materials  (see  also  Figure  5.13).  This  could 
be  the  explanation  for  the  absence  of  magnetization  cooling  effects. 

In  a  recent  report  Lawless  and  Munson  (Section  IV)  have  derived,  using 
standard  thermodynamic  relationships,  (dy/dT)  plots  from  the  temperature 
dependence  of  AT^..  The  (dy/dT)  behavior  for  the  CdCr204  samples  as  obtained 


-113- 


from  the  ATj-  measurements  is  that  dx/dT  is  negative  at  T  <<  T>.j,  becomes 
positive  around  Tj^j,  becomes  negative  for  T  >  T.;,  and  that  dx/dT  is  highly  field 
dependent.  For  ZnCr204,  on  the  other  hand,  dy/dT  alway  r''"’.lns  negative  and 
shows  a  negative  peak  around  These  data  may  be  compared  with  the  dy/dT 

values  as  obtained  from  direct  measurements  of  x  shown  in  Figures  5.10  and 

5.11  for  all  the  four  samples  (namely  P,  PD,  DDCP  and  DDPS) .  As  nay  be  noted 
from  Figures  5.10  and  5.11,  all  the  samples  of  ZnCr204  do  show  a  behavior 
similar  to  that  for  CdCr204  samples  In  contrast  to  the  results  from  the 
magnetocaloric  measurements.  We  note,  however,  that  the  susceptlblllt;.'  (  y) 
measurements  were  made  at  200  G  whilst  the  magnetocaloric  experiments  were  done 
at  much  higher  fields,  which  makes  It  difficult  to  make  more  definitive  comments 
since  the  susceptibility  behavior  has  been  found  to  change  significantly  for  the 
ZnCr204  powder  by  going  from  0.2  kG  to  10  kG  (Figure  5.4). 


V.6.b  Specific  Heat  vs.  x  Data 

We  now  compare  the  x  results  with  those  from  the  specific  heat  analysis 
(Section  IV).  As  shown  in  Figures  5.12  and  5.13,  the  Tjg  values  from  the 
specific  heat  curves  are  always  higher  than  those  obtained  from  the  x  data.  For 
a  simple  anttf erromagnetic ,  it  is  known  that  the  magnetic  specific  heat  is  given 
by  the  Fisher  relation  (Fisher  1962). 


9(XT) 

8T 


(V.3) 


where  A  is  expected  to  be  a  slowly  varying  function  of  T  around  Tjvj.  Figures 
5.18  and  5.19  show  the  temperature  dependence  of  for  CdCr204  and 

ZnCr204  (PD),  respectively.  l!cre  A  has  been  assumed  to  he  tvmperatvire 
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Figure  5.19. 


Comparison  of  the  specific  heat  Cp  (continuous  curve)  with  the 
magnetic  contribution  Cm  (open  squares)  for  ZnCr20i+;  the  parameter 
A  was  set  to  10.  Circles  denote  the  same  data  (C^)  with  a 
different  A  and  after  renormalizing  Cp  and  C^  at  (T/Tf^)  =  0.8. 


Inversely  as  the  square  of  the  temperature  and  is  small  around  Tm  and  (iii)  the 


Fisher  relation  remains  valid  for  T/T^  =  1.0  ±  0.2  for  a  simple  ant  if erromagne t . 


We  further  assume  that  at  (T/Tfj)  =  0.8,  Cp  =  0^.  Choosing  A  =  17.84  Instead  of 


the  arbitrary  number  20  and  scaling  the  0^  data  such  that  it  peaks  at  (T/T»j)  =  1 


yields  data  points  shown  by  closed  circles.  Similarly  for  the  ZnCr204  PD, 


choosing  A  =  5.47  Instead  of  10  and  scaling  the  data  the  closed  circles  are 


obtained.  An  observed  significant  disagreement  for  the  ZnCr204  sample  than  for 


the  CdCr204  sample,  when  compared  with  the  Cp  data,  ind'.-ated  tliat  the  spin 
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ordering  in  ZnCr20i|  samples  are  more  complicated  than  for  the  CdCr20L,  samples. 
Recalling  that  the  (1/x)  versus  T  plots  of  Figure  5.15  indicated  persistence  of 
short-range  order  up  to  T  “  5  T|>j,  we  deduce  that  this  postulated  complicated 
spin  ordering  might  be  the  explanation  for  the  behaviors  found  in  Figures  5.15 
and  5.19.  On  the  other  hand,  the  CdCr204  PD  seems  to  behave  in  a  much  more 
expected  way  as  also  reflected  by  the  l/y  versus  T  plots  of  Figure  5.14. 

Lawless  (Section  IV)  has  analyzed  the  specific  heat  data  by  taking  into 
account  the  Debye  (D),  spin-wave  (SW)  and  Schottky  contributions.  That  is, 

Cp  =  (mg^^  +  mp)T3  +  b/T^  (V.4) 

or 

CpT2  =  (mgy  +  mj^)T5  +  b  (V.5) 

Thus  a  plot  of  CpT^  versus  is  expected  to  be  a  straight  line,  and  is  indeed 
found  to  be  so.  Let  us  look  at  the  implications  of  this  analysis  in  the  light 
of  our  susceptibility  results. 

A  Schottky  type  analysis  of  the  Cp  data  for  T  >  Tjq  (l.e.  a  linear  plot  of 
CpT^  versus  T^)  yields  b  which  is  proportional  to  n,  the  density  of  the 
paramagnetic  spins.  (Note  that  the  ordered  spin-;  behave  like  free  paramagnetic 
spins  above  T  "  Tj^)  .  The  ratio  of  b's  for  the  powder  disk  (PD)  and  the  DDCP  or 
equivalently  the  ratio  of  n  (PD)  and  n  (DDCP)  is  found  to  be  1.29  for 
CdCr204  and  1.09  for  ZnCr204  samples.  These  values  are  in  excellent  agreement 
with  our  susceptibility  data  which  independently  yield  the  values  1.41  for 
CdCr204  and  1.03  for  ZnCr204. 

Schottky  analysis  of  Cp  for  T  <  Tfj  does  not,  however,  give  results 
consistent  with  the  susceptibility  results.  In  the  model  comprising  of  "free" 
and  perfectly  "ordered"  spins,  it  is  reasonable  to  assume  that  the  free  spins 


.  •«.*  V*  •s." 


give  rise  to  the  Schottky  tern.  Witli  this  assunption  it  is  found  that 
n( PD) /n( DDCP)  =  6.98  for  CdCr2n4  and  4.72  for  2nCr204,  Mote  that,  this  ratio 
involves  the  numbers  of  "free"  spins  only,  l.e.,  it  exclvides  the  "ordered"  spins 
specifically.  From  our  susceptibility  analysis  we  find  n( PD) /n( DDCP)  =  .022  for 
CdCr204  and  0.130  for  ZnCr204.  The  value  from  the  Schottky  analysis  is  “320 
times  larger  for  CdCr204  and  “35  times  larger  for  ZnCr204.  As  shown  bv  our 
susceptibility  plots  and  also  suggested  by  Lawless  (Section  IV),  there  may  be  a 
large  number  of  spins  showing  only  short-range  order  Instead  of  a  long-range 
order.  The  contribution  of  these  spins  at  T  <  T^vj  might  be  one  of  the  reasons 
for  the  above  discrepancy.  To  this  end,  we  wish  to  note  a  few  points: 


(1)  Low  temperature  Schottky  analysis  yields  linear  plots  only  above  T  =  3.3  K 
for  ZnCr204  PD,  T  *>  2.2  K  for  ZnCr204  DDCP,  T  “  1.7  K  for  CdCr  204  DDCP  and  T  “  1 
K  for  CdCr204  PD.  Tlie  coefficient  b  which  depenls  on  n  and  6^  (<5  =  the 
two-level  spacing)  must  therefore  Involve  fairly  large  6  to  account  for  the 
deviations  below  the  temperatures  mentioned  above. 


(li)  C/T^  versus  T  plots  at  low  temperatures  yield  a  similar  behavior  for  the 
PD  and  the  DDCP  samples  whilst  the  number  of  "free"  spins  in  the  DDCP  samples  is 
much  higher  than  in  the  PD  samples. 


(ill)  The  Debye  contribution  is  negligibly  small  for  temperatures  below  Tnj. 
This  can  be  seen  from  the  ratio  of  the  spin-wave  densities  as  obtained 
rigorously  by  Lawless  (Section  IV)  and  those  obtained  by  just  dividing  the 
^™SW  ^  ”0^  coefficients  as  shown  in  Table  5.8. 
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Table  5.n.  Spln-w.ive  Densities 


Spinel 


CdCr  2O14 
ZnCr 20^ 


ncu  (PD)/n<.u  (DDCP) 


(mj;y  +  '^D^PD'^^^SW  ^D^DDCF 


5.«7 
I  6.00 


V.6.C  Shift  in  Tv 


As  shown  in  Flj^ures  3.12  and  5.13,  there  Is  .1  clear  shif"t  in  the  valne 
of  Tfj  as  the  powder  is  subject  to  different  heat  and  stress  t  reatnents .  To 
understand  this  effect  qualitatively  we  need  to  diy^ress  the  followin>;. 

Anderson's  theoretical  considerations  (Anderson,  1956)  drew  attention  to 
the  role  played  by  the  next-nearec t  neighbor  interactions  In  the 
antiferromagnetic  ordering  of  the  B  ions.  The  B  lattice  is  built  up  from  units 
of  four  nearest  neighbors  situated  on  the  corners  of  a  tetrahedron,  and  such  a 
tetrahedron  shares  one  B  site  with  an  adjacent  B  tetrahedron.  For  only  the 
nearest-neighbor  interactions  B-X-B,  the  ground  state  Is  determined  bv  the 
condition  that  the  spins  in  each  tetrahedron  cancel.  That  Is,  negative  nearest 
neighbor  in  a  completely  normal  spinel  does  not  load  to  ant Iferronagnet Ic 
long-range  ordering  since  the  ground  state  Is  largely  degenerate  duo  to  the 
frustration  of  the  lattice.  The  interaction;  between  the  R  sites  in  different 
tetrahedra  (B-X-X-B)  can,  however,  give  long  range  ordering  and  the  N’eel 
temperature  will  depend  on  these  interactions  in  a  complicated  way.  Denoting  h 
J  and  K  as  the  strengths  of  the  nearest  and  next-nearest  Interactions,  the 
magnetic  propert'es  can  bo  accounted  for  approximately  by  the  Heisenberg  spin 
haralltonlan  of  the  form 


■  -2.: 


2J  I  Hi  -  Si  -  2K  S  Si-S;^  -  g;:H^  E  Si^ 


(  \’  .  6  1 


The  number  of  nearest  neighbor;  In  a  perfect  spinel  in  6  whilst  that  of  next 
nearest  neighbors  is  36;  J  anc*  K  in  the  above  expression  denote  the  average 
strengths.  One  can  calculate  Tjj  with  this  hamiltonlan  and  study  its  dependence 
on  the  parameter  (K/J).  It  is  to  be  noted  that  such  a  philosophy  has  already 
been  applied  to  spinels  which  exhibit  ferromagnetic  behavior  (Baltxer  et  al. 
1966).  For  ferromagnetic  orderings,  it  was  found  that  slight  changes  in  the 
parameter  (K/J)  can  shift  the  Curie  temperature  by  a  few  d.egrees.  It  should  be 
noted  that,  according  to  the  Ohio  State  University  group,  many  features  of  the 
X-  and  Cp-data  on  these  spinels  can  be  accounted  for  by  distorting  the 
tetrahedra  of  tlie  spinels  to  remove  frustration  (Section  IV). 

In  short,  we  believe  that  the  value  of  T[vj  is  a  sensitive  indicator  of  the 
changes  in  the  next-nearest  neighbor  interactions. 

V.6.d  Field  dependence  of  x 

The  magnetic  field  dependence  of  the  specific  heat  and  the  nagnetocalnric 
heat  changes  seem  to  have  been  established,  although  it  has  not  been  pursued 
systematically.  Similarly,  the  field  depen, !ence  of  x»  although  noted,  has  not 
been  studied  systematically.  It  is  puzzling  to  find  that  the  x  versus  T 
behavior  depends  on  the  field.  For  CdCr204  PD,  the  magnetic  field,  raises  tb.o 
Neel  temperature,  contrary  to  expectation.  We  do  not  have  any  explanation  for 
these  features  right  now  and  would  like  to  investigate  this  aspect  to  bring 


together  the  susceptibility,  specific  he.it  and  magnetocalor ic  results, 
the  possibilities  is  that  these  compounds  exhibit  low-dimensional 


One  of 


antiferromagnetic  behavior.  However  the  low  temperature  (T  <  Tj^.)  neutron  or 
X-ray  diffraction  studies  and  single  crystal  samples  are  needed  to  confirm  t 
conjecture  which  provides  a  new  viewpoint  for  the  unlriuo  properties  of  tliese 
spinels . 
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VI.  THEORETICAL  STUDIES  OF  THE  SPINEL  LATTICE  :  GINZBURG-LANDAU  MODEL  AND 


MONTE  CARLO  OCMPUTER  STMULAlV.'M.S 

A.  Introduction 

This  section  is  devoted  to  the  theoretical  research  on  the  spinel  materi¬ 
als  under  investigation  in  this  contr.oct.  In  particular,  calculations  of  the 
specific  heat,  magnetic  properties,  and  dielictric  constant  v/ere  carried  out 
using  a  variety  of  complementary  calculational  approaches,  including  a  phe¬ 
nomenological  Gin2burg-Landau  model,  as  well  as  first  principles  Monte  Carlo 
computer  simulations. 

In  this  section  we  present  first  an  overall  picture  for  the  spin  ordering 
in  the  spinel  lattices,  then  discuss  the  results  of  the  Ginzburg-Landau 
calculation  of  the  magnetic  ordering  and  its  effect  on  the  dielectric  con¬ 
stant,  which  enables  us  to  explain  the  striking  dielectric  anomaly  reported  in 
section  IV.  Next  we  summarize  the  results  of  our  extensive  Monte  Carlo 
calculations  on  the  undistorted  cubic  phase  of  the  spinels,  and  then  present 
our  calculations  on  the  tetragonal  or  distorted  form  of  the  spinel  in  which 
the  frustration  has  been  reduced.  Finally  we  describe  how  low  temperature 
paramagnetic  tails  may  arise  naturally  from  the  composite  structure  of  the 
CdCr204  and  ZnCr204  materials.  In  each  part  we  will  discuss  the  connection 
with  the  relevant  experiments. 

B.  Spin  Ordering  ^  the  Spinel  Materials 

Spinels  are  ferrite  materials  with  the  generic  formula  AB2R4  where  R  is  a 
group  VI  element  (0,  S,  Se,  or  Te),  while  A  and  B  are  metal  cations  as  shown 
in  Fig.  6-1.  The  particular  spinels  CdCr204  and  ZnCr__04  are  high  specific 
heat  spinels  of  present  interos  ,  in  which  the  P  site  elements  Cr  form  a  3-D 
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Flfi.  6-1.  Spinel  Lattice;  B-sites 
have  magnetic  Cr  ions 
coupled  by  O^"  ions. 


Fig.  6-2.  Magnetic  B-sites  in  lattice 
structure  of  spinel.  Note 
chains  in  successive  layers. 


frustrated  magnetic  system.  It  has  been  found  that  there  is  a  phase  transi¬ 
tion  into  an  antiferromagnetic-like  state  at  low  temperature. 

We  start  with  a  brief  review  of  the  salient  features  of  the  spin  lattice 
structure  of  the  B-site  frustrated  spinel  lattice  as  shown  in  Fig.  6-2.  The 
results  obtained  so  far  have  suggested  an  attractive  picture  for  the  ordering 
phenomena  in  the  B-site  spinels,  CdCr204  and  ZnCr204.  The  results  have  shown 


that  at  least  two  types  of  magnetic  correlations  are  present,  antif errcmag- 
netic  and  paramagnetic,  and  that  frustration  and  the  presence  of  strong 
spin-lattice  coupling  play  an  important  role  in  the  anomalously  large  specific 
heats  and  thermal  co.nductivities  (Patton). 

In  addition  the  transitions  in  these  materials  were  seen  to  have  a  pecu¬ 
liar  nature  in  which  the  spins  order  weakly  in  a  lattice  which  has  a  high 
degree  of  frustration.  This  has  great  importance  for  our  u.nderstand i.n,g  of 


these  systems,  since  it  means  that  large  numbers  of  spins  can  remain  unordered 
below  the  transition,  resulting  in  anomalously  large  specific  heats,  and 
furthermore,  that  distortions  of  the  lattice,  which  remove  the  frustration, 
can  couple  strongly  to  the  spins,  thus  leading  to  dielectric  anomalies  and 
large  thermal  conductivities  due  to  spin  energy  being  transported  through  the 
spin-phonon  interaction. 

More  recently,  the  experimental  work  done  under  this  current  proposal  has 
revealed  much  about  the  properties  of  the  B-siie  spinels,  CdCo204  and  ZnCo204. 
In  particular,  the  recent  work  includes  1)  successful  fabrication  and  struc¬ 
tural  characterization  of  various  ceramic  phases  including  powder,  compressed, 
and  densified  samples  of  CdCo204  and  ZnCo204,  reported  in  section  III,  2) 
experimental  work  on  the  magnetization  and  susceptibility  which  reveal  an 
antiferromagnetic  transition  at  Tjj  z  5-llK,  section  v,  3)  experimental 
research  on  the  specific  heat  and  magnetocaloric  properties,  section  IV,  and 
4)  experimental  work  on  the  dielectric  anomaly  at  Tj^,  which  reveals  a  novel 
coupling  of  the  spin  and  lattice,  section  IV. 

The  main  experimental  results  we  will  focus  on  in  connection  with  the 
theoretical  calculations  are  the  following: 

a.  The  fabrication  experiments  (III)  give  clear  evidence  that  the  spinel 
powder  of  nominal  grain  size  10  microns  combines  with  additional 
mineralizers  to  produce  a  material  having  extra  paramagnetic  spins  at 
low  temperatures,  among  other  properties. 

b.  The  susceptibility  measurements  (V)  on  the  pure  spinel  powder  reveal 
little  low  temperature  param..  'netic  Curie  tail,  while  the  grains 
reacted  with  lO'/-  columbite  or  5'/.  Nb205  reveal  a  low  temperature 
paramagnetic  tail. 

c.  The  magnetocaloric  experiments  (IV)  reveal  no  trace  of  hy  ster ‘=■5  is, 
implying  the  phenomena  involved  in  the  low  temperature  transitions  in 
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the  Zn  and  Cd  spinels  are  second  order  in  nature.  The  pattern  of 
adiabatic  demagnetization  cooling  and  heating  at  different  temperatures 
is  indicative  of  antiferromagnetic  ordering  in  the  Cd  samples,  while 
all  samples  show  evidence  of  low  temperature  antiferromagnetic  spin 
waves. 

d.  Measurements  (IV)  have  discovered  a  striking  anomaly  in  the  dielectric 
constant  at  the  antiferromagnetic  transition,  which  clearly  shows  that 
the  ordering  of  t!'.e  frustrated  spin  system  has  a  large  effect  on  the 
structure  of  the  lattice.  The  dielectric  constant  in  large  magnetic 
fields  (Fig.  4-8)  shows  a  clear  shift  in  the  transition  for  one  sample. 

e.  Another  sample,  the  CdCo204  disk,  did  not  have  a  sharp  specific  heat 
maximum  and  had  a  strong  magnetic  field  dependence  to  8X/9T  below  Tjj. 

In  connection  with  results  a.  and  b.  above,  in  section  F  of  our  report  we 
show  that  by  allowing  our  clusters  to  have  a  finite  surface  area  we  also 
develop  (theoretical)  low  temperature  paramagnetic  tails.  Our  calculations  in 
sections  D  and  E  show  that  the  transitions  have  no  hysteresis,  and  are  charac¬ 
terized  by  basically  antiferromagnetic  interactions  as  found  in  c.  We  calcu¬ 
late  explicitly  the  dielectric  anomaly  of  d.  in  section  C  and  show  that  in 
general  it  shifts  slightly  in  a  magnetic  field.  Finally,  our  calculations  in 
section  D  on  the  frustrated  spinel  material  show  that  the  transition  as  seen 
in  the  specific  heat  is  broad  and  relatively  featureless,  which  may  explain 
the  sample  of  e.,  while  our  calculations  in  section  E  show  that  relaxing  the 
frustration  via  a  tetragonal  distortion  that  often  occurs  in  the  spinel  system 
gives  rise  to  an  enormous  specific  heat  singularity  as  seen  in  three  of  the 
samples. 
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C.  Ginzburg-Landau  Model  and  Calculation  of  the  Magnetization  and  D i 1  e c t r i c 
Constant. 

We  start  with  a  review  of  our  phenomenological  Ginzburg-Landau  model  which 
enables  us  to  interpret  and  relate  the  many  experimental  results  on  specific 
heat,  magnetic  susceptibility  and  dielectric  constant,  and  w:i  ch  has  motivated 
much  of  our  other  calculations.  In  particular,  this  approach  enables  us  to 
explain  the  very  intriguing  dielectric  anomalies  at  the  magnetic  transition. 

1.  Dielectr  1C  Constant  of  Spinels 

In  general,  the  spinels  exhibit  extraordinarily  high  dielectric  constants 
e  at  and  above  room  temperatures,  but  have  rather  weak  temperature  dependence 
at  low  temperatures.  In  the  present  case,  experimental  results  (section  IV) 
have  revealed  a  large  peak  in  the  temperature  derivative  of  the  dielectric 
constant,  de/dT,  at  the  Neel  temperature  (see  Fig.  4-8),  which  suggests  a 
coupling  between  the  charges  and  the  lattice. 

We  give  first  the  simplest  picture  for  the  dielectric  constant  in  an 
insulator  and  explain  why  it  can  be  affected  by  a  magnetic  transition,  which 
IS  very  unusual.  The  dielectric  constant  arises  from  the  polarization  of  the 
medium  via  virtual  transitions  from  a  filled  electronic  state  to  an  empty  one. 
Since  the  energy  denominator  in  such  a  transition  is  of  the  order  of  the  gap, 
the  dielectric  constant  will  increase  as  the  gap  decreases,  becoming  infinite 
in  the  limit  that  the  gap  vanishes  and  a  metallic  state  is  reached.  A  simple 
representation  of  the  band  structure  of  the  spinel  is  shown  in  Fig.  6-3;  the 
filled  valence  band  is  separated  by  a  gap  at  the  zone  boundary  k  -  tir/aQ  from 
the  empty  conduction  band.  The  basic  point  is  that  any  perturbation  which 
affects  the  band  structure  will  have  an  effect  on  the  diel-jitric  constant.  The 
expression  for  the  dielectric  const.ant  in  an  insulator  may  be  written  CJ.man, 
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If  we  consider  the  static  limit  here,  ftu)  -*  0,  replace  the  energy  denominator 
with  an  average  gap  A,  and  utilize  the  oscillator  strength  sura  rule,  then  Eq. 
(6-1)  takes  the  simple  ^rrm 
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Where  uip  is  the  plasma  frequency.  We  may  use  the  result  (6-2)  to  understand 
qualitatively  what  happens  when  the  ordering  of  the  spins  couple  to  a  distor¬ 
tion  of  the  lattice.  Neutron  scattering  experiments  (Oles,  et  al.,  1976) 
show  a  variety  of  patterns  in  the  low  temperature  spin  ordering  of  the  B-site 
spinels:  however,  a  characteristic  feature  is  a  doubling  of  the  unit  cell  from 
ap  to  2ao.  If  the  spins  couple  to  the  lattice  then  a  periodic  potential  acts 
on  the  electrons  which  induces  bandgaps  at  multiples  of  Tr/2ap  as  shown  in  Fig. 


6-4. 


Since  the  valence  band  in  the  insulator  is  completely  filled  and  the 


conduction  band  is  empty,  the  gap  induced  at  ±Tr/2aQ  has  little  effect; 
however,  the  effect  at  the  original  zone  boundary  iir/ag  is  to  incriase  the  gap 
already  present.  This  clearly  lowers  the  electronic  energy  since  occupied 
states  are  lo-..-ered  in  energy,  while  unoccupied  conduction  states  are  raised  in 
energy.  Thus  it  is  clear  that  the  average  gap  A  is  increased.  We  may  write  A 
in  terms  of  a  constant  part  Aq  and  a  temperature-dependent  part  6(T) 

A(T)  -  A  +  (5  fT)  (6-3; 

''o' 

where  6(T',  due  to  its  couplir.g  to  the  magnetic  order  parameter,  will  have  a 
temperature  dependence  below  Tj,j  as  shown  in  Fig.  6-5,  Taking  the  derivative  of 
Eq.  (6-2)  then  gives  the  result  shown  in  Fiq.  6-6  for  3e/3T.  The  experi.ments 
on  the  dielectric  constant  therefore  represent  a  very  important  probe  of  the 
spin  and  structural  ordering  in  the  spinel  materials  and  have  been  pursued  on 
all  available  samples. 
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Fig.  6-5.  Temperature  dependence  of  Fig.  6-6.  Temperature  derivative 
bandgap  6  arising  from  of  dielectric  constant 

''oupling  to  .magnetic  vs.  temperature. 


ordering. 
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2.  Ginzburg-Landau  Model  of  Spinels 

In  order  to  derive  the  proper  Ginzburg-Landau  dercripticn  of  ihe  rpinel 
materials,  we  note  that  since  the  experimental  susceptibility  (section  V)  and 
adiabatic  demagnetization  (section  IV)  are  both  characteristic  of  an  antifer- 
romagnetic  system,  the  natural  long  range  order  parameter  of  the  system  is  the 
sublattice  magnetization.  Furtheririore,  with  both  magnetic  and  electric  fields 
present,  then  the  sublattice  magnetization,  the  total  magnetization,  and  the 
polarization  should  all  appear  as  relevant  order  parameters.  We  write  dow-  a 
Ginzburg-Landau  theory  with  a  free  energy  containing  couplings  of  these  order 
parameters  and  use  it  to  analyze  the  system. 

We  start  by  expressing  the  free  energy  functional  F  in  terms  of  the  order 
parameters  L,  M,  P,  which  are  the  difference  of  sublattice  magnetic  moments, 
the  total  magnetic  moment,  and  the  electric  polarization,  respectively: 

F  =  F3f[L]  +  Fpjn[M]  +  F^LP]  +  Fap[L,M]  + 

where 

Faf[L]  :  J  dr  [  a|Ll^  +  b|L|^  +  ]  (6-5) 

IS  the  energy  associated  with  the  sublattice  magnetization  and  near  the  Neel 
temperature  Tj,j,  .1  ;  ao(T-TN)/TN, 

is  the  magnetization  energy  including  that  associated  with  the  external 
magnetic  field  H  with  a  >  0, 

(  2 

dr  \|P|  -  E*~ 

IS  the  dielectric  polarization  energy  including  the  coupling  to  the  e:<terr:al 
electric  field  F,  and  finally  the  coupling  terms. 
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Eq.  (6-8)  gives  the  coupling  between  the  magnetic  order  parameters  L  and  M;  we 
have  G  >  0  since  the  antiferromagnetic  and  ferromagnetic  states  tend  to 
exclude  each  other.  Eq.  (6-9)  describes  the  interaction  between  the  lattice 
distortion  involving  charje  displacement  and  the  sublattice  magnetication  !.. 
Both  these  couplings  are  required  by  1)  the  behavior  of  the  magnetic  suscepti¬ 
bility,  and  2)  the  experimental  connection  between  the  dielectric  constant  and 
the  magnetic  phase  transition. 

The  Eqs.  (6-4)-(6-9)  have  a  rich  structure  even  within  the  mean  field 
theory  approximation  (MET)  in  which  fluctuaticns  are  neglected.  Since  we  are 
concerned  primarily  with  the  equilibrium  state  of  the  system  and  the  nature  of 
the  phase  transition,  this  gives  a  reasonable  first  approximation.  Ve  include 
fluctuation  effects  in  the  next  section  where  we  consider  the  experimental 
data.  In  the  MET  we  have  the  following  equations  resulting  from  minimication 
of  the  free  energy  E[L,M,P],  namely, 


6F  (5F  _  6F  _ 
6M  '  6L  '  6P 


(6-10) 


thus 

2cM  -  H  +  2G1L1^  -  0 
2aL  +  2b|L|^L  +  2G1M1^L  +  2riP|^L  =  0 
2XP  +  2r|L|^P  -  E  --  0 


(6-11 
(6-12  ) 
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In  the  case  G  =  E  -  0,  we  may  identify  the  unrenormalized  magnet;  -  and 
dielectric  susceptibilities,  Xom  Xoe- 
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P  -  Xoe^  ” 


Then  from  Eq.  (6-13)  we  obtain  the  polarization  as 


P  5  XgE  =  - 7  --  - - 

2X+2r|Lr  U2ryoell^l 


(6-16) 


while  from  Eq.  (6-11)  the  total  magnetic  moment  is 


M  =  XmH  ^  - 2  ^  ^ - 2 

2a+2G|L|  l  +  2GXomli^l 


(6-17  ) 


The  sublattice  magnetization  L  is  determined  by  Eq.  (6-12),  which  gives  L  -  0 

above  the  transition  (a  >  0),  and  the  solution 


a  +  blir  V  GCXmH}^  +  r[XeE)^  --  0 


(6-18) 


below  the  transition  (a  <  0).  Eq.  (6-18)  is  a  non-linear  expression  in  |L|^ 

due  to  the  dependence  of  Xm  ^nd  Xe  on  |L  1^  as  indicated  in  Eqs.  (6-16)  and 

(6-17).  Defining  the  dimensionless  quantities:  sublattice  magnetization  1, 
temperature  t,  electric  field  e,  magnetic  field  h,  and  coupling  constants  y 


and  g 
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we  may  then  rewrite  Eq.  (6-18)  for  as 

- g--,  +  1"  +  t  =  0  (6-25) 

(1+2y1  )  (l+2gl  )" 

Solving  Eq.  (6-25)  yields  the  sublattice  magnetization  as  a  function  of  the 
temperature  and  fields  L  -  L(T,H,E).  We  note  that  when  y  is  not  equal  to 
Eq.  (6-25)  IS  a  quintic  equation  for  1^,  namely, 

2  22  22  2  22  2  22 
(1  +t)(l  +  2yl  )  (l  +  2gl  )  +  ye  (lt2gl  )  +  gh  (l  +  2yl  )  --  0  (6-26) 

The  value  of  |L|2  obtained  from  the  solution  of  Eq.(6-26)  determines  the 
dielectric  constant  e  and  the  magnetic  susceptibility  as  follows, 


where  from  Eq.  (6-16), 


and 


e  :  l  +  4TTXe 


Xe 


Xf)o _ 

l+2r|L|^Xoe 


(6-27  ) 


(6-25) 


Xja  --  - -  (6-29) 

1+2G|L|  Xom 

From  the  solution  to  Eq.  (6-26)  it  is  clear  that  as  the  magnetic  field 
increases  at  constant  temperature,  the  value  of  L  decreases.  The  field  at 
which  L  =  0  defines  a  critical  magnetic  field  h^^  r  (-t-ye2)/g,  cr 

Ho2  =  ^  (6-30) 

^'<ora 

at  which  the  system  is  driven  into  a  paramagnetic  stale  ly  the  external 
magnetic  field,  i.e.  a  kind  of  spin-f’.-'p  '■ransition. 
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As  the  electric  field  goes  to  zero,  or  if  the  coupling  constant  y  is  very 
small,  the  second  term  of  Kq  .  (6-26)  can  be  dropped.  Then  F)q  .  16-26)  reduces 

to 

(l"+t)(l  +  2gl^)^  +  gh^  I  0  (2  .i: ; 

The  numerical  solution  of  Eq  .  (6-31)  '  shown  :•  Fig.  6-7  Fig,  6-3  I's  two 

temperatures  below  the  transition,  t  =  -l/j  and  tr-i/-d;  the  normalized  suMat- 
tice  magnetization  squared  IL/Lol^  is  plotted  versus  the  normalized  mja.g.netic 
field  squared  (H/Ho)^  for  various  values  of  the  ar.tif erromagnetic- 

f erromagnetic  coupling  constant  g  where  Lo  is  the  value  of  L  at  H:0  and  Hq  is 
given  in  Eq  ,  (6-30) . 

It  IS  interesting  to  note  that  for  larger  values  of  g  the  curves  show 
re-entrant  behavior  which  suggests  a  discontinuous  (first  order)  transition 
those  cases,  We  may  determine  the  critical  value  of  g ,  ,  at  which  re¬ 

entrance  appears  as  the  value  of  g  for  which  the  slope  of  !L1^  as  a  function 
of  first  becomes  infinite  at  Hq  .  Equivalently,  we  may  use  (6-31)  to  express 
h^  as  a  function  of  1^  and  require 

^  =0  (6-22) 
dl  IrO 

With  the  result  that 

-  4|  I  I  ^  '^■33  ) 

Thus  the  I'endency  toward  a  first  order  tr.ansition  in  a  field  increases  as  the 
temperature  is  lowered.  Eq  .  (6-33)  predicts  a  .’..inimum  value  of  of  1.-4  at 
zero  temperature.  On  the  other  hand,  near  the  transition  Tfj  ,  t  '>  0,  and 
second  order  transition  in  field  will  always  occur. 


Fig.  6-9.  Change  of  dielectric  constant  vs.  external  magnetic  field  at  tern 
perature  T  :  3/4  Tjj  with  critical  g  =  gc  V  =  1  2. 
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Fig.  6-10.  Change  of  dielectric  constant  vs.  external  marnetic  field  at  te; 

perature  T  i  3/4  Tm  with  y  -  2  and  g  greater  >  he  ''ritic.il 


Using  Eqs.  (6-27),  (6-28)  and  (6-31),  we  can  plot  AC/fl-co)  ■/?..  as 

shown  in  Figs.  6-9  through  6-12,  for  various  fi::''-d  valurrs  of  the  tercpora  t  ur  e 
and  the  coupling  constants  g  and  y.  The  same  re-entrant  behavior  is  found  for 
g  values  greater  than  g^.  We  also  note  that  at  and  '  ve  the  critical  field 
Hq  the  dielectric  constant  anomaly  is  removed  for  all  values  of  g  and  y  ,  i.e. 
Ae  r  0.  This  occurs  because  above  Hq  the  magnetic  field  has  suppressed  the 
sublattice  magnetization  (L  -  0),  and  thus  the  interaction  between  th-'  charge 
and  the  lattice,  ylPhjLI^,  has  effectively  been  removed. 

To  summarize  this  section,  we  have  found  a  novf'l  Oinzbur.g-Landau  free 
energy  which  describes  tne  unusual  coupling  betw«»en  ^he  dielectric  constant 
and  the  antiferromagnetic  transition  in  the  B-site  spinels,  CdCr204  and 
ZnCr204.  We  next  consider  a  detailed  comparison  of  the  theory  with  the  recent 
experiments  in  section  IV  on  the  magnetic  field  dependence  of  the  dielectric 
constant. 

3.  Comparison  with  Dielectric  Constant  Measurements:  Effects  of  Fluctuations 

In  order  to  make  detailed  comparisons  with  the  experiments  it  is  necessary 
to  include  the  effects  of  fluctuations  on  the  expression  for  the  dielectric 
constant.  This  is  clear  from  the  fact  that  the  dielectric  constant  starts  to 
decrease  above  the  actual  transition,  which  occurs  at  the  point  of  maximum 
slope  (see  Fig.  4-1). 

The  simplest  approximation  that  includes  fluctuations  is  the  Hartr'='e 
approximation.  Consider  the  free  energy  F 

F:aL"  +  bL'  -  (a«-  bL*" )  L*" . 

In  the  Hartree  appro.xi ma 1 1 t.he  L'*  term,  w.hich  represents  the  interactions. 
IS  evalu.ted  using  its  a-erage  value,  namely,  2L^<;L^>.  Eq.  (6-3t)  t)'.  :. 


becomes 
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35) 


Self-consistency  is  imposed  by  calculating  tbe  average  of  L-  by  using 
Eq.  (6-35) 


^L  >- 


(Tr 


-0F 


L" )/  i  rr 


e 


(6-3  6) 


Together  with  Eq.  (6-35)  this  gives 


a  =  a  +  b/a 


Eq.  (6.-37)  may  be  solv explicitly  for  a.  The  reqi:ireTn.^nt  that  the  limit 
a  t,  a  be  achieved  at  te.mr--'rat ure  far  above  T^j  so  that  L-  approaches  the 
mean-field  result  implies  that  the  correct  solution  of  Eq.  (6-3^'  is 


(6-38  ) 


The  sublattice  magnetization  is  thus 

2  2  1/2 

L  :  2/[a+(a  +4b)  '^  ]  (6-39) 

In  the  vicinity  of  the  transition  temperature  Tj^,  where  the  reduced 
temperature  t  a  (T-Tfjl/Tjvj  satisfies  |t|  «  1,  the  Ginzburg-Landau  coefficient  a 
IS  giv-’n  by 

a  -  a  t  (6-‘'*0) 

0 

Although  Eq.  (6-40)  is  accurate  near  Tj.j,  it  gives  a  non-vanishing  derivative 
for  at  T  t  0  which  violates  the  second  l.iw.  A  be'  t^^r  f.'rmula  whu'h 
correct!'/  interpolates  between  the  high  and  Ts'W  t emperaturi'  limits  is 


a  :  a  [  (Tm/T)2-i  1/2  i  a  t(Ut/21  -  a  r 
0  0  0 


thus  we  have  as  a  result,  in  zero  magnetic  field 


T  +  (t"+  b) 


Rewriting  the  coupling  constant  y  in  (6-23)  in  terms  of 


Y  5y  (e  -1)  -  (e  -1  )  •  [Fa  T  /4TTb]  (0--+3  i 

10  o  0  N 

and  using  Eqs.  (6-42)  and  (6-43)  with  Eqs.  (6-27)  and  (6-25)  mvos  o\;r  final 
expression  for  the  dielectric  constant 

.  2vi(€^-  :  )i2.Cq 

l  +  2yi 

where  in  general  Eq  in  (6-44)  depends  on  t'^mperature  also.  Frcm  fitting  the 
experimental  curve  for  CCN(9/1)  in  Fig.  4-3,  we  find  e,j  -  1?.C15  -  0.2  t, 
yj:  4.33*e”5,  and  b  =  0.000.  The  theoretical  expression  from  Eq.  (6-44)  is 
plotted  together  with  the  data  roints  fr'-'i  experiment  (from  Fig.  4-1)  in  Fig. 
6-13. 

When  the  magnetic  field  is  included,  provided  th:  effect  of  electric  field 
IS  negligible,  the  sublattice  magnetization  is  given  by 


g  (cH  ) 


Tf  (t^+0.06)^'^^ 


1  + 


^  g 


(6-h5) 


x+  (t^+0.06)^'^^ 


where  H  is  in  units  of  Tesla.  From  fitting  the  curve  (in  Fig.  4-5),  g  and  c 
are  determined  to  be  g  ^  3.37»e‘3  and  c  :  2.-7  T.  We  note  that  the  coupling 
constant  g  is  small  enough  that  when  the  magnetic  field  increases  from  O.'^l  T 


to  1 

5.16  T  the 

change  in  t  is 

about  O.l  out  of  20. 

The 

field 

dependence  of 

the 

dielectric 

constant  from 

Eqs.  (6-44)  and  (6-45) 

IS 

s  h  e  w  n 

in  Fig.  6'  - 1  e . 

The 

theory  f  o  t 

ttie  magnetic  field  'dependence  is  also 

1'" 

good  a 

r  0  m  w  r.  t  w  ^  *  I  'l 

the  experimerFs,  indicating  that  e  in''reases  as  the  field  in;:  •vises 


D.  Monte  Carlo  Calcul.: '  .  ons 

1.  I :::  I'  action 

The  results  obtained  in  previous  work  on  this  contract  as  well  as  the 
G 1  n 2 b u r g  -  L a n d  a '  1  theory  of  the  previous  s e c 1 1 o n  h. a v e  all  s u p a e s  t e d  a  p i c  t u r e 
fvr  the  ordering  transitions  in  the  B-site  spinels,  CdCr-,C4  and  ZnCr-O^, 
involving  two  types  of  magnetic  correlations,  an  tiferrornagnetic  and  paramag¬ 
netic.  Frustration  and  the  pr"»nce  of  strong  spin-lattice  coupling  appear  to 
play  an  important  role  in  the  anomalously  large  specific  hea^.s  and  thermal 
conductivities.  In  addition  the  transi*.  ions  in  .‘..•re  m.aterials  seem.ed  tr 
involve  weak  spin  ordering  in  a  lattice  which  '  as  a  !i.,jh  degree  of  frustra¬ 
tion.  This  means  that  large  numbers  of  spins  can  remain  unordered  below  the 
transition,  resulting  in  anomalously  large  specific  heats,  and  furthermore, 
that  distortions  of  the  lattice,  which  remove  the  frustrati  .^n,  can  couple 
strongly  to  the  spins,  thus  leading  to  dielectric  anomalies  and  large  thermal 
conductivities  due  to  spin  energy  being  transported  through  the  spin-phoncn 
interaction.  The  most  powerful  way  to  examine  these  effects  is  through  "'onte 
Carlo  simulations  of  the  actual  spin  lattice.  Before  treating  the  actual 
magnetic  ion  lattice  formed  by  the  B-sites  in  spinels  CdCr204  and  ZnCr204,  we 
describe  the  computer  techniques  and  obtain  calibrations  for  comparison  pur¬ 
poses  on  a  simple  cubic  lattice  I  sing  model  spin  system. 

2.  Monte  Carlo  S imulations 

The  I-Ionte  Carlo  method  uses  a  compu'i’er  to  model  the  thermal  behavior  of  a 
physical  system  by  simulating  the  random  microscopic  processes  which  give  rise 
to  the  thermal  equilibrium  state.  The  Monte  Carlo  raethC'J  and  its  applica t ion.' 
in  various  fields  have  teen  der-onited  in  detail  (P.nder,  1979, 198e).  Tlie 
basic  idea  is  to  in  luce  the  system  ti  g--'  'hn.ush  .i  nmsl  's  t  a;  psmts  in 


phase  space,  calculating  equilibrium  thermal  quantities  along  the  trajector 


Many  different  sampling  techniques  m-'v  be  used.  The  .-irnplest  and  most  ef; 
cient  one  is  called  "importance  sampling",  which  chores  random  states  accor 
ing  their  importance  in  the  problem,  e.g.,  according  to  a  Soltr.m.ann  dirtrih' 
tion  in  the  canonical  ensemble. 

Choosing  the  Is’.nc  model  as  an  ''"ample,  the  usual  ram.plin.i  rroc'lure 
the  following; 

I.  Start  with  an  initial  state  at  temperature  T,  and  calculate  itr 

magnetization  M  and  energy  U; 

II.  Select  spin  to  flip,  ca’culate  the  "transition  probability" 

e;-:h(-AE/k3T),  where  AE  is  the  chan.ge  in  energy  a.ssociated  with  the 
spin  flip; 

III.  Compare  the  probability  with  a  random  number  .x  (Ofxil).  If  the 
probability  is  larger  than  x,  flip  the  spin;  otherwise  do  not; 

IV.  Repeat  ii  and  iii  for  j  from  1  to  N,  where  N  is  the  total  number  of 
spins  in  the  lattice,  then  cal--  late  E  and  U  again; 

V.  Repeat  ii,  lii,  and  iv  many  times.  Calculate  the  averages  of  M  and  U, 

which  are  just  the  equilibrium  values. 

There  is  some  flexibility  in  the  selection  of  the  starting  point  in  th 
procedure.  The  initial  state  can  be  chosen  as; 

(a)  an  arbitrary  state; 

(b)  an  ordered  state; 

(c)  the  equilibrium  state  at  a  temperature  close  to  the  tempera'ur? 

considered. 

The  last  one  is  used  m.ost  often  since  with  it  leads  te  equilibrium,  faster. 
The  free  energy  F  and  entropy  C  may  be  obtained  by  integration  :f  t 
following  t  h  e  r  m  0  d  y  n  a  .m  1  c  s  r  e  1  a  1 1  o  n  s ; 


(dS/dTlH  "  (dU/dTjTVT 


■  •”  rvdT:H 


The  expressions  for  S  and  F  are  provided  in  terras  of  integrations  either  f: 
0  to  T  or  fr.ira  :xi  to  T.  They  are 


VT 


S(T.  H)  -  S(0.  H;  +  f  (  TJ/dT  ’  )udT  ’ , 

Jo 

‘  -  t-t-iA-l  -  (dU/dT'  ;v;dT’/hoT’ 

kpT  Kb'I’  -■'■3  '-0 


or  eq’Jivalently 


FiT,  S(':o,  if) 


•f: 


i  />T 

U  dd/k^T) 
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f b-50  ) 


%1/kT 


SiT,  H)  :  S(0D,  H)  +  UiT,  H)/T  -  kg  f  U  d  f :  'koT) 

J  0 


(6-51  ) 


[IS.  an.l  S6-5’)  are  often  more  convenient  to  use  because  it  is  usually 

.nown  tor  an  s-sp.n  Ism g  mclel  that 


H):  I-Tg  ln'2s+l) 


(6-52  ) 


hen  F,T,;i;  ar.f  d.T.H,  can  be  easily  ccmputed  from  U(T,H), 

Iri  ;ur  oalcul.ition,  tlie  inost  r-:'!  :vant  experimental  quantities  are  the 
lerii.c  heat  ’  ant  the  ra  ,::netic  susceptil  ility  X'  k/e  can  calculate  C  directly 


C  (dlJ/dTlH 


'  6-53  ' 


-  •■■■•')  /  kgT-  (6-5-+) 

from  the  same  calculation  that  gives  oU  ■ 
E'l.  {'■■-5-V)  u.suaHy  gives  a  more  reliable  result 


The  susceptibr.; '-•/  can  be  defined  in  two  different  The  difference 

susceptibility, 

X  ;  M/Tf  (&-:!;} 

or  the  differential  susceptibility, 

X  ::  (WdH)'p  (6-5t  I 

The  latter  is  equivalent  to 

X  -  ug(<M^>  -  <M>^ )/kgT  (0-;^) 

which  IS  often  more  useful  than  (6-56)  in  practical  appliiiations. 

In  the  usual  procedure  described  above  for  Monte  Carlo  campling,  a  com¬ 
plete  cycle  '■  steps  (*i),  (iii)  and  :iv)  is  called  one  Monte  Carlo  step  per 
spin.  The  nu'r.ber  of  MC  steps  per  s:  ..  needed  is  dM  '-rmined  by  how  fast  the 
results  converge.  Good  convergenc-.-  r  .  luires  that  the  difference  between  the 
average  quantity  before  and  after  one  step  is  much  smaller  than  the  quantity 
.tself.  Usually  the  number  of  MC  steps  needed  ranges  from  10^  to  10^  for 
typical  spin  systems. 

In  a  Monte  Carlo  calculation,  the  systems  used  are  obviously  finite,  with 
the  number  of  spins  N  usually  ranging  from  10^  to  10^.  Since  a  phase 
transition  can  only  occur  in  an  infinite  system,  we  can  not,  strictly  speak¬ 
ing,  obtain  any  singular  result  from  our  calculation.  With  periodic  boundary 
conditions,  when  the  correlation  length  5  gets  larger  than  the  linear  dimen¬ 
sion  L  of  the  system  near  the  phase  transition,  the  fluctuations  become  "over 
correlated"  and  +he  calculated  results  are  rounded  and  smeared.  These  are 
called  finite-si2e  effects,  and  may  in  fact  be  used  to  determine  the  range 
correlations.  For  this  reason,  finite-size  scaling  theory  is  an  very  useful 
tool  to  investigate  properties  near  a  critical  point  (Binder,  1934). 

An  important  fact  should  be  stre.ssed  here,  r. 'imely,  the  '.fference  in  "Tie 


behavior  of  the  order  parameter  obtained  from  a  Monte  Carlo  calculation  near  a 
first  order  phase  transition  from  that  obtained  near  a  hisher  order  phase 
transition.  At  a  first  order  transition,  due  to  the  existence  of  metastable 
states,  the  order  parameter  shows  hysteretic  behavior,  while  at  a  higher  order 
transition  it  does  not.  As  shown  in  Fig.  6-15,  this  difference  is  '.ften  used 
to  distinguish  the  first  order  transition  from  higher  ord::  ones  (  Mouritsen, 
1984). 


a)  First  order  transition, 
dashed  lines  indicate 
metastable  states. 


(b)  Continuous  transition. 


Fig.  6-15.  Order  parameter  4>  vs.  temperature  T 


Compared  to  other  methods,  Monte  Carlo  simulations  have  two  obvious 
advantages:  first,  the  method  may  be  used  with  any  system  with  a  well  defi.ned 
hamiltonian,  and  second,  it  is  a  non-perturbative  method  which  relies  on  no 
approximations  of  the  underlying  fundamental  physical  laws,  an  approach  which 
may  bring  about  new  and  unexpected  discoveries. 


3.  Calculations  on  a  simp.'o  cubic  lattice 


To  make  sure  our  programs  are  correc\,  we  first  worked  on  a  simple  cubic 
lattice  as  a  trial  to  calibrate  our  tec:,  niques,  since  the  results  can  be 
compared  to  previous  work. 

The  lattice  we  used  in  the  present  case  v.'as  a  10x10x10  ferromagnetic  Ising 
(s  ;  1/2)  spin  array  with  periodic  boundar'-'  ;•< '.di  i  ions,  althoi'.gh  lattice  sices 
ranging  from  5x5x5  to  37x37x37  were  ''.so  sV'.died.  The  calculations  were  dene 
over  a  range  of  temperature  T  -  0  to  T  -  2.01V  (A.0T(-  for  U  and  M),  where 
T(,  =  4.51J  and  J  is  the  coupling  constant  b  ‘ween  nearest  neighbor  spins.  For 
these  calculations,  we  utilized  a  modest  ICO  MC  steps  per  spin  at  each  temper¬ 
ature,  which  proved  to  be  adequate  for  appr.-'ach  to  equilibrium  in  this  case. 

First,  the  energy  and  magnetization  were  evaluated  directly  from  the 
average  over  all  the  configurations  obtain-  1  in  the  random  procedure  at  each 
temperature;  the  results  are  shown  in  F-  u  6-15  and  Fig.  6-17,  respecti'.-ely. 


Fig,  6-16.  Energy  vs.  temper.itur e 
for  10x10x10  cubic  f er rom. icnet  1C 
Ising  spin  lattice. 


As  expected,  the  energy  is  equal  to  -30GGJ  fj  times  the  total  numher  of 
nearest  neighbor  bonds)  at  low  temperature,  and  approaches  cero  at  tempera¬ 
tures  much  above  the  transition,  indicating  that  the  spins  are  highly- 
decoupled  from  each  other.  Similarly,  the  .magnetisation  equals  the  total 
number  of  spins  at  low  te.mperature,  s.howing  that  all  the  spins  are  ordered  in 
one  direction,  but  is  nearly  sero  at  high  temperature,  corresponding  to  a 
completely  disordered  state. 

From  U  and  M  at  different  temperatures,  we  computed  the  free  energy  and 
entropy  from  Eq.  (6-50)  and  (6-51).  They  are  shown  in  Fig.  6-18  and  Fig. 
6-19.  We  find  a  reasonable  behavior  for  S,  as  well  as  F,  from  the  results. 
At  lo'iv  temperature,  S  equals  k2ln2  since  the  ground  state  is  2-fold  degener¬ 
ate.  At  high  temperature,  S  equals  Nkgln2  where  N  -  ICiOO  is  the  total  number 
of  spins. 


Fig.  6-13.  Free  energy  vs.  temper-  Ihg.  6-1).  uiitr.p-/  vs.  tempc-ru  turo 

ti  t  ij  r  0  for  1 0  A 1 C'  /.  I  ■  j  c  ’. !  b  i  c  f  o  r  r  o  rr^  r  -  i  s' r  1  i  '  a  I C'  a  ] ^  i  f  ‘  r  r  c  rr.  r,  r.  c  t . 
netic  I.sing  spin  '.  .ttice.  Ir.ing  .-^pin  l.itM-'f*. 


To  resolve  the  phase  transit if-ii 
at  T-,  the  specifi:  heat  was  calcu¬ 
lated  d'-ectly  from  the  derivativ.'  of 
energy  with  respect  to  temperature. 
The  result  is  plotted  in  Fig.  6-20, 
which  indicates  that  C  ex.hibits  a 
peak  anomaly  at  T^.  The  divergence 
of  C  IS  rounded  by  finite  size 
effects,  as  expected.  A  small  shift 
of  the  peak  to  lower  temperature  is 
also  seen  in  the  results. 

To  investigate  the  order  of  tran¬ 
sition,  we  calculated  U  and  M  both 


upon  heating  up  the  lattice  frcm  low 

temperature  and  by  cooling  it  down  ^.OQ.  Specific  heat  v; 

from  high  temperature.  The  results  perature  for  10x10x10  cubic 
turned  out  to  be  identical,  well  magnetic  Ising  spin  latti'^e 

within  the  numerical  uncertainty,  for  ulated  using  Eq.  6-5"' 

both  directions,  thus  no  hysteresis 
was  observed.  This  indicates  the 

transition  should  be  continuous  instead  of  first  order,  a  conclusion 
known,  of  course,  for  the  Ising  model  in  three  dimensions. 

All  the  above  results  are  consistent  ’.''ith  previously  published  w:r 
literature  Therefore  we  can  conclude  from  cur  trial  calcula 


analyses  that  the  Monte  Carlo  meth  1  and  our  program  .a- e  well  calirr 
further  studies  on  the  magnetic  spinels,  wliose  proper'uei-  can  re  r  1  ■: 


analogy  with  the  results  reported  here  for  'he  Simr’e 


E.  Monte  Car  1  o  Calculations  on  Frustrat '^■1  Spinel s 

1.  Spinel  Magnetic  Lattice 

As  discussed  in  previous  secticr^s,  the  r.agnetic  ions  on  the  3-sites  in 
spinels  CdCr204  and  ZnCr204  form  a  frustrated  spin  system.  The  st'^ucture  is 
shown  in  Figs.  6-2  and  6-21.  Since  v.'e  tro  focusing  here  on  just  tlie  ma.cnetic 
properties,  only  the  magnetic  atoms  are  draivn.  We  can  see  ir Fig.  6-^1  that 
the  whole  lattice  is  tuilt  up  of  inter c'^nnected  tetr.ihedral  units.  Each  s"  .n 
IS  part  of  two  adj.icent  units  which  means  that  for  each  tetrahedral  unit  there 
are  two  spins,  with  eight  units  in  a  pru.Mtive  cell  of  the  crystal. 

We  consider  the  basic  loops  composing  the  lattice  triangular  pladuettes 
which  form  the  faces  of  the  tetrahedral  units.  Around  each  of  these  loops, 
there  is  an  odd  number  of  ant  if errorna jnetic  bonds  which  has  the  consequence 
that  spin  alternation  around  the  loop  can  not  be  perfectly  satisfied,  as  s.h:7,n 
in  Fig.  6-22.  Considering  the  Ising  ca.se,  at  least  one  spin  in  the  loon  c  ui 
be  flipped  without  changing  the  energy.  This  is  called  "frustration"  in  spin 
systems  (Toulouse,  1977  ).  The  spin  which  can  be  flipped  freely  corresponds  t.  - 


a  (paramagnetic)  defect,  wnich  represer'-  an  ;  r.  f  :  t'r.d  in  ;:r:':r.d 

state  configuration.  Disorder  eKis*-,-  •.  ■  '-'"e  spins  c^n  f"? 

chosen  up  or  dovm  randomly.  Similar  cor  ■  ".r.  i:,:.  ha'/c  Leon  applied  to  the 

spin  glass  phase,  which  has  a  number  "i  onus-:-;',  proper t  le.s .  Thus  we  sh-uld 
not  be  surprised  to  find  some  peculiar  '  ".-ults  arisrnr  also  fr-m  the  case  if 
frustration  in  spinels. 

We  now  do.'^cribe  the  application  of  ‘he  Monte  Carlo  7/' •■■'dime  revr  iwed  in 
the  previous  section  to  CdCroO/^  and  ZnCr204.  As  mention'^d  before,  we  v/ill  use 
Monte  Carlo  simulations  to  investigate  ‘heir  magnetic  properties.  The  lattice 
formed  by  magnetic  atoms,  shown  in  Fig.  C-C  or  6-21,  is  frustrited  due  to 
triangular  plaquettes  of  the  corner  sh  ring  tetrahedra.  We  can  also  think  of 
the  spin  lattice  as  layers  of  spins  with  parallel  spin  chains  in  each  layer. 
The  Hamiltonian  we  will  use  is  the  Isir.;?  m.odel  with  nearest  neighbor  antifer¬ 
romagnetic  interactions; 

H  r  j  S  (5-58) 

<  1.  J  > 

where  J  >  0,  <i,j>  denotes  nearest  neighbors,  and  the  spins  and  Sj  take  on 
the  values  il.  We  n 'te  that  the  Ising  rodel  appears  *^  0  be  appropriate  to  the 
CdCr204  and  ZnCr204  systems  (rather  ;  an  the  isotropic  interaction  of  the 
Heisenberg  model,  for  example)  based  on  what  is  known  of  the  nature  of  the 
spin  ordering  (Oles  et  al.,  1976). 

The  next  section  is  devoted  to  the  ground  state  properties  obtained  from 
the  Monte  Carlo  calculations.  We  note  ir.  advance  that  although  .oome  of  these 
results  gave  encouraging  agreement  with  thu-  experiments  of  sections  IV  and  V, 
they  did  not  provide  evidence  of  a  char-'  transition.  This  suggest'^1  that  tr.'=‘ 
actual  spinel  compo.site  materials  are  able  to  reduce  the  frustration  in  rjp.e 
way,  which  will  be  discussed  in  the  ne”t  pa-t,  rect’ooi  F. 


Sirice  the  system  is  frustrated,  tlie  bonds  can  net  be  all  sa'is:i-d  -i‘ 
once.  However,  it  is  clear  that  "he  lowest  eneriiy  for  a  .-ir.dle  tetrahe  .I.’- :  r 
can  be  achieved  by  setting  the  four  spins  on  it  as  tv,'o  up  .and  tv/o  do-wn  un  ru;^ 
different  ways).  Tlie  ground  state  f-.r  the  whole  lattice  is  form’-l  by  pu’*i:'. 
many  such  tetrahed:  a  together.  The  lo-west  energy  for  .1  single  tetrahe  .ls'.  n  1: 


-2J.  If  there  are  a  total  of  !T  spins,  since  eic)'.  tetrahed 


4»(l/2)  =  2  spins,  there  will  be  11/2  t  r  fsahedra.  There!. re  the  .ground  ■*  i*u 
energy  Eq  is  equal  to 


•  (-2J)  -  -NJ 


The  ground  state  so  constructed  is  non-periodic  and  v“ry  dnsorderc-d .  If 
one  looks  at  a  single  chain  of  spins  running  through  the  spinel  lattice,  it 
can  have  any  pattern.  Obviously  the  ground  state  degeneracy  is  very  large.  A 
lower  bound  can  be  obtained  by  setting  all  chains  (in  x-y  plane)  in  a.nt.- 
parallel  configurations,  or  in  other  words,  satisfying  all  the  bonds  in  the 
x-y  plane,  then  counting  the  two  fold  degeneracy  of  each  chain.  The  number  of 
chains  are  (N/2)2/3_  and  each  chain  has  2  states,  yielding  a  lower  bound  of 

The  upper  bound  arises  as  follow-,  Divide  all  tetrahedra  i.ito  two  set.s 
such  that  the  tetrahedra  in  one  set  only  connect  to  tetrahedra  in  the  ot.her 
set.  Then  the  upper  bound  for  the  ground  state  degeneracy  can  be  obtained  bv 
setting  the  spins  on  each  tetrahedron  in  one  set  (which  ai'e  mutually  discon¬ 
nected'  ac  two  up  and  two  down;  the  number  of  possibili  ties  comes  out  to  be 
(jM/4,  an  u;  :  er  bound  since  the  spin  configurations  on  the  tetrahedra 

in  the  other  set  may  or  may  not  be  appropriate  for  the  ground  state.  If  b 


denotes  the  ground  state  de.gcnoracy,  we  have 


The^e  rounds  on  .Uun.-neru  :y  yi-.-iU  ..e  fellow  •'<; 


entropy  per  spin  s 


:k3.1:u:  :/M 


0.437  K, 


A  further  aFrro:<imate,  but  useful  ‘:'Sti.'nite  for  D  :.'i 


20X  (Anderson,  1950): 


N/2 

D  71  (3/2) 


which  im'lies  a  f:round  state  entropy  "'er  spin  of 


Sq  --  kg‘2-ln(3/2  )  :  0.203 


In  all  the  cases  considered  the  zero  temperature  entropy  has  .i  finite  valt 
a  result  of  the  frustration  of  the  spin  ordering.  As  we  dis''uss  in  secti 
below,  a  lattice  transformation  may  reduce  or  eli"iinate  the  frustration 
allow  the  entropy  to  go  to  zero  at  zero  t -"perature. 


3.  Monte  Carlo  C  alculations  of  tl:e  Thermal  and  M  ar;  net  ic  y  r  o  pe  rt  le  s  of 
Frustrated  Spinel  Lattice 

Using  Monte  Carlo  simulations,  the  thermal  and  mapn'=*ti'.'  pr'perties  of 
system  have  1  een  calculated.  Most  of  our  calculation.'  ■■'■'“re  ■•(''n"  on  a  l.’i 
of  5x5x5  unit  cells  containing  5x5xhxl&  -  2000  spin.'-,  wi^.h  ,terati;n  ■ 
ranging  from  2000  to  20000  M'l’  .steps. 

Direct  aver'iging  over  many  random  configurations  yiel  is  ■h'^  magnet:':. 


M(  D  and  ener.'y  U(T)  as  functions  of  temperatur M  appear 


r  to  gi 


zero  macr.olic  fiji';  ever  th 


.exper.'.  t  urc 


with  the  ant.‘''''rr  •  gn'.-'  ■  itur 


•- fv  .V  .V 


the  grovind  stat-'  er.'r  r!’  .  '.ip.  be  irr.’:;'- . :a ' ely  ce.l 'iile  te  J,  v/rj:  :.', 

"q~  . 3 

Our  Monte  Carlo  result  (t-t-o  is  in  goi  ;  agreement  v/ith  +h"  ontrepy  cbta 
from  the  estim. a''  e  :  i  ’■.he  ground  star-  degeneraoy,  E'l.  (6-63). 

One  cf  th.e  meat  ■us'f'.l  and  in  t  rr‘--s  ting  pr:  per  *  ic  .s,  the  specific  hea 
can  be  obtained  by  differentiating  t'(T  ;  the  •  ■'.'ult  is  shovn  n  Fig.  6-23 
low  fcro.id  maximu.m  in  C  is  found  at  T  h  0.6 J  with  Cr^^.j;-p:  C.28,  vihore  c  a  C/’j 
the  specific  heat  per  .spin.  normal  cinfr’i.strated)  .spin  ordering  transit 
typically  involve  maximum  specific  he->yr  an  order  :-f  magnitude  larger  ’ 
this.  For  such  a  low  a  pea.K,  it  appear.^  unlil'oly  tl'.at  a  sharp  ’■■-.ase  tra 
tion  IS  actually  present.  A  low  specific  heat  peak  could  also  occur  if 
transit  i  n  v.-  first  order,  but  near  a  second  order  critical  print  w 


Fig.  6-23.  Specific  heat  vs.  temp-  Fig.  b-24.  Energy  vs.  temperatur 
erature  for  txbxbxlj  spin  frus-  for  heat.ng  and  cooling  runs  at 

trated  spinel  lattu'e  at  H  -  C'.  H  r  n.  h  v ^s '  s  ebs-.u-ve.;. 


gives  weak  fluctuations.  However,  i';  that  care  a  l.;r, 
hysteresis  woulh  be  expeoteh.  In  Hi.j-  rl'!*:s  :  "  T 

and  cooling  r’lns  are  shcv/n.  No  hyst  -resis  i.c  observ'h 
transition,  if  there  is  one,  be  second  or  !er  or  hic'^n^r. 
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T 6  t, c* rn p 0 r  i  ^  li r o*  t.  ^ P o ci f  i  -  I'l  ’  * "t.  nui ^  <•  h ^ -  ^  j-i  j;  l f  o  r.  ' ^ ^  “  ■(■/ 1 

lattice  size  for  a  conventional  ic-iv, r  ;  ic  Iciii;:  transit-  -.  If  '. •_ 

the  temperature  of  the  maximum  of  the  .1  .-render. ce  of  t-----'-  di:: 
lattice  recces  was  calculated  and  the  results  plotted  in  "ih-  h-2f.  Su 
ingly,  we  find  good  agreement  with  finite-size  scaling  th-eory  -'r  large  L 
the  3-d  Ising  coherer:^  length  exponent  v  =  fi.f.3  i.<-  u.sed,  v  ■•sti’ig  a 
nuous  phase  transition.  It  is  interest!:  -  to  rote  that  'ih-, ...  '  '-ir-.-ases  v 
unlike  the  case  of  the  simple  Ising  magnet  where  it  decr-aases  WLth  L 
speculate  that  this  notable  difference  may  arise  fr'-'m  the  imur*  r  .ited  naf 
the  spinel  spin  system. 

.tnother  interesting  property,  the  magnetic  susceptibilit  y  x  (in 
field)  "  a  function  of  temperature,  is  calculated  fro.r,  the  f  luc  t-.:a  lions 


The  results  are  plotted  in  Fig.  6-27  and  shoiv  a  pat •. ..-rn  similar  to  that 


•-.'pc  L.  J'.  i  V 1 0  r  rince,  t  r  a ‘.e-l  below, 

alue  at.  whitb  j  larpe  portion  :;'  rpins  are 
sbj'wn  in  an  '■■npand  >  J  I'cate  in  F’l;;.  0-30, 

n  r  by  Mae  field  (this  is  more  clearl 

V  e  r  s  u  s  H  at  O'  ;  C . '  3  o  1  it  is  s ;  1 1  n  in  7 :  b . 

at  '3  ;  is  {jiven  by  t.he  dashed  li: 
'  ■■r.es  ip.ore  r 'ended,  o”  ' r. i  via  1 1 1'  '.'.'ir  e  its  ; 
)  structure  in  M  is  the  folic v: : n p, :  a  fie  1  d 
to  flip  a  spin  (in  the  ahsonce  of  thermal 


4.  Discussion 


From  the  specific  heat  calculci*  i'  '  ■■■'•?  do  net  fint  evider.'co  :;r  o 

sharp  transition.  The  specific  heat  i:  lead  shows  a  lov,'  m  a i  r.'.  u  m  and  a  s’,  c.v 
decrease  above  Tf,j.  The  absence  of  a  t.ansition  is  consistent  v/ith  the  early 
conclusion  of  Andor ■  on;  no  long  range  c'ler  can  ezist  in  the  3-site  spinels 
only  antif errorragnetic  coupling  betv/cen  nearest  neighbors  is  invvlvod  ',.\:ide:  - 
son,  1956).  The  large  value  of  the  specific  heat  above  the  peak,  on  the  cuh:-: 
hand,  implies  considerable  local  order,  a  result  seen  in  co.me  of  t.'ie  ther.m^’ 
and  magnetic  measurements  is  sections  IV  and  V. 

To  try  to  under st.and  v.'hat  kind  of  ''rder  is  d''velopin.e  ’'o  note  that  t.herc 
is  one  quantity  which  is  the  same  t  hro-' -’h'-'-t  the  ■'•.••'stal  at  h' •  -mrenature 
namely,  each  tetrahedron  has  2  up  and  2  down  spins.  That  .v. a n s ,  fr'^'m  thii 
fact,  we  may  extract  a  kind  of  long  r.-r.,te  ordering.  Assign  a  .number  to  eac." 
tetrahedron  as  follows;  if  the  i^*^  tetrahedral  has  2  up  and  2  down  sp.nr 
associate  ;  1  with  it;  otherwise  ’et  fj  =  -3/5.  Eefine  a  "t.itrah'' 1’ 

order  parameter  F  as 


a 


Where  we  have  normalized  the  order  parame*  -r  by  the  number  of  tetrahedra,  '.’/2 
At  T  =  0,  each  tetrahedron  has  2  up  and  2  down  spins,  so  =  I  for  all  i 
So  ue  have 

F(T-C)=1  (6-66; 

As  T  goes  to  infinity,  all  micro  states  have  the  same  probability  to  appear 
Thus  the  probability  to  obtain  a  tei' r  ib-^'lron  with  2  up  and  2  down  versus  tha  + 
for  all  other  type  of  tetrahedra  i;;  j  to  10,  or  equivalently',  3  to  5.  This 
means  we  have 


■V_  v 


Kwwnn  wun  nn  WT  PUT  wnwn-itr  V'JVW\’sr>.V\.’VW\JVL’VV«^\/Vl'V 


Which  IS  the  reason  why  was  ■•^^''.■r  ter!  ■-■■'  -t/h  f  .  r  -.1’  n :  n  t  ’jp/2  h' 
conf igurati'^ns,  as  discussed  above. 

W 1 1.  T.  n c  'i!  0  f  1  r.  1 1. 1  1 V  n ,  P  i  - .  .  i  ^  i*  .i  ^  \  n  ^  -  a  --a  ^  r  - 

to  1  at  sere  temperature  a'  !  '/anishing  at  h:.:h  t .: -per.' tur  '  re  ;:r.d  cut  y 
F  does  near  th"'  temperature  ■•■her'^  *he  sr'^^:":'  t.^at  ‘  ^  .s  m.a/Ci.m.um, 

calculated  it  in  the  Monte  Carlo  c.ul  eti-'n.  t"'.--  ;  ?  ■  -  ."uo-cn  .,n 

6-34.  The  curve  shows  that  F  '.r.cr  'se  monct  ihica'. ly  as  *;'.e  tem'-cat 
decreases  as  e:<P'=‘cted  ‘^or  other  '■rd-''r  para.m.e*  ers,  ’out  a  .''r.kiric  pr'perto 
the  uctreraely  .slow  decay  above  t'..,  •..,  -  wcH  .  r^'l  i '.  i v  'l , '  .'m.  Ml  value 
"^max-  hard  to  find  a  sin.^ul 's  :  y  .n  F  t  any  '  ■  m. a  *  ve  from. 

gras;>.  Calculations  of  F  in  toth  haac.ng  and  rooling  prcc.^sses  also  shows 
hystereses,  which  is  consistent  vath  wh.i  t  happen*--:  in  si  nil';  .lalculaliohs 
the  energy  U. 

The  long  range  ordering  intrc -^-’ced  above  is  boss--!  on  pr sper ti-as 
small  spin  clusters,  not  on  conventicr  ■■.I  magnetic  rnc  s.  In  some  sense, 
IS  not  a  typical  type  of  magnetic  ordering,  as  -discussed  by  An-d-ers'-n,  but 
ordering  does  have  the  physical  meaning  of  being  assc-'iated  with  the  dror 
entropy  from  2^^  for  the  totally  disordered  -ta^e  to  about  (3/2)*'^/*^  fcr 
"ordered  tetrahedra".  We  interpret  t:  is  result  as  eviden-'e  for  strong  1 
correlations  in  spin  orderin.g  over  a  1  ;rge  rai"  e  i.n  ‘  emperatures,  whi'h 
difficult  to  classify  as  s-rdin.-ry  mag.nelicaticn,  but  apn-ear  in.-t-^^id  ;r.  re  1 
a  spin  glass  type  of  freezing.  We  are  not  sure  .f  F  ..'  the  b-ist  o:  only  ■ 
parameter  at  this  moment;  -.iicother  definition  rniglM  show  a  sharper  i^b. 
transition.  Furt:  -r  calcula  ions  of  the  correlatio".  functions  -.'d  ‘  !ie 


would  be  very  illuminating. 


5.  Conclusion 


Detailed  Monte  I’arlo  C3lculati'':n'  on  truotrateJ  r-p  nilr.  iiave  pro'/i  l--d 
extensive  results,  v/hich  show  a  lark  1  a  sharp  phase  trar.sition  due  to  t'.e 
complete  frustration  of  the  spin  lattice,  but  the  p'  ;""  ility  of  extensive 
local  order.  CcmparisO'n  of  the  si  eoifio  heat  ■-'a  1 'u la t : :  ns  v.-if  tiie  experi.men- 
tal  results,  Fips.  4-9  and  4-10,  in-lic."- 1  is  tliat  the  I'dC’  -0,.  vder  'di-v 
good  agreement  wit;i  our  results  for  the  frustrated  latt::-,  •.vliile  the  CT.'h,'/ll 
and  the  f.vo  corresponding  Zn  samples  appear  to  have  si',.  >  p  'ransitijns.  In 
the  next  section  -.ve  will  see  that  a  tetragonal  dis'or'ion  of  the  lattice, 
reported  in  ZnCr^O,',,  but  not  CdCrTO/^,  leads  un  a  sharp  p'it.'e  transition. 

Comparing  the  susceptibility  from  our  Mcnt.'  Carlo  calculation  with  the  tne 
from  experiments.  Fig.  5-3,  we  can  see  the  results  agree  qualitatively  with 
each  other  except  for  paramagnetic  tails  in  the  data  at  low  temperature  which 
are  especially  pronounced  in  the  powder  samples.  We  will  shov/  in  next  se:ti:n 
that  If  we  change  the  periodic  boundary  cc.-.  .itKn  used  to  .;late  an  in:::  ue 
crystal  to  free  boundary  conditions  appropriate  to  a  finite  sited  grain,  a 
paramagnetic  tail  will  appears  in  tlie  raiculated  suscepLibil i ty. 


1.  I  nt  r  Oil  u  i' t  :oj] 

As  mentioned  in  the  last  sect::;!,  the  Monte  '.e  1 a  t :  nt 

the  fully  frustrated  spinel  do  not  dive  a  sharp  transit:  jn,  v-:..  ;-.  is 
in  some,  but  not  all  of  tlie  eta  fFia.  >  an  !  -'‘-I'',  ■  :r "  1  ■  ' .  :  • 

basically  two  ways  to  reduce  the  frustra  ion  and  to  hat'-  a  ;  has---  t: 

tion.  One  is  to  introduce  second  neares'  neijihibor  in  tp-r  a ,  th-:-  .s 

to  allow  the  lattice  to  undergo  a  crystal  transforrnat  on  w,.:o:-;  roduios  tin- 
symmetry.  We  do  not  consider  the  firri  alternative  at  this  t.me  because  1; 
there  is  little  solid  information  on  the  nature  of  the  seci  .d  n-arest  nei,:hb'“r 
interactions  other  th -.n  they  are  at  leas'  a:-  ■rde-r  c:  m  tg;.:  t  u  d  .i  smaller  than 
the  nearest  neigl'ibor  interactions,  and  ■  neutron  sc,- 1  ten  n,-’  t' 1 1 -a  r. :  e  indi¬ 
cates  that  Zn'.r2Q.^  ond  many  other  spin.  und  -n  go  a  lat*i  .'  ;-  tnanri  ■  r;:  i 
below  the  magnetic  transition  temperatur  .  The  frustration  in  the  syst.-.m  .may 

be  total’y  or  par"  lly  removed,  allo'.v:  ■  ■  :r-d...:'.ry  magne'::  lent:  ra.nge  or.h^r 
to  be  achieved.  The  ."  mplest  and  poss;.''/  mcs''.  important  dis'orti...  s 
tetragonal  one  which  makes  the  latti  :-'r  '-onstan*  c  in  z-'- '  rc'c  t  i-;-n  un. _l 
.1  :  1'  ir.  :-;-y  pl.ane.  '.v'e  consi-Jer  thi.s  i:  d  ,.f  -l.stortnan  .';:■■■■  ;t  is  •  t:n  1  in 
ZnCr20,,  et  al.,  1976).  There  ar-'-  *v;o  Mff>^’"ent  -lases  s'ssibl-^: 

Ac  <  0,  and  Ac  ■>  0;  we  find  the  nfund  s.t.-.'e  prop'^rtK-  !  narr;'  -  'ut  Mn.t  ■ 
Carlo  : a  1 c  u 1 a  1 1 o  n  s  f  c  r  b  o  t  h  o  f  t  h  e  s  e  c  a  ■  s . 

For  convenience  in  ’he  follow. n..’  discussions,  l-^t  : -ipr '-sen*  •  b.-' 

coupling  -ronstant  betvr'een  neares’  neii.’tt-.rs  in  each.  ;  i.in-'  and 

sent  the  -ne  be’w-.-en  r.-- irest  nei):hb'--t.  in  !i‘-eren’  plan-  .-s  in  ;-';g.  •  ’ 

The  HamT-’.-nian  '.'an  b.e  •wti'.'en 


2.  The  Ground  Stjter 

In  the  Ac  <  0  case,  o...'  •.•/pects  th,‘. Jj  •:  J2’  It  means  that  the 
between  nearest  neirjh''ors  in  eren'  planes  should  be  satisfied  f; 
anti-parallel  alifinment,  as  sh.own  in  Its.  6-30.  Because  of  coupli.ng 
chains  in  a  f  I'^ne  throush  chains  the  next  plane,  parallel  s 
required  in  the  whole  plane.  Bo  the  .cr-..-;;!.;!  state  confiqura  1 1  ;n  .s 
ing;  spins  are  parallel  in  the  sam'  -lane  and  an  1 1 -  par  ;  1 1  ,il  in 
Dianes. 

In  tins  case,  one  ''.^n  easily  ub'  am  the  nround  sta'e 
Cjuntin.c  the  bcii'ls,  which  ('ives 


U'TtO  ) 


state,  except  for  an  overall  flip  of;  .11  ', 


■A  '  •^  r'l  I  '  *■>  '^  ^  t  -i  •*  ,“*  H  .'■■ 


equals  to  2,  which  means  the  ground  sta ; 


Sq  -  (kg‘l;.:)/i: 


The  frustrat; -n  has  been  removed  from  i.'e  s".'’; 


spin  IS  uniquely  determined  if  its  r, 


In  the  Ac  >  0  care,  one  has 


the  nearest  n e i g h b o r s  in  each  p  1  a n 


parallel  spin  chains,  as  s h, c- n  ; n  " i g . 


each  other  on  one  chain  moke  a  near---st  u'^ichbo’’  •"'■n  on  an.' 


next  plane)  either  up  or  ccwn.  Thar  i.os  th.e  .'has?  :f  a  "p  o 


determines  wheiher  the  first  spin  n 


the  state  of  the  other  cha.nr,  Th-;  o  ;;  ;,  the  nuur:.;  :  n; 

the  following:  anti-parallel  ordered  "un  cha;:!.',  paralle:  ‘  ■'.•a: 

the  same  plane  and  perpendicular  t  o»  each  ctluer  ;n  a  f  e n  t  p 
mutually  independent  phases. 

Ttie  ground  state  energy  can  be  '  im.'.  ‘  -d  .n  a  way  s;::iil;\r  * 
case.  It  IS  easy  to  see  that  bonds  ■'  upling  'different  nianes  gu 
ener,gy.  All  net  energy  comes  from  bon.l:  .nsi:-  ila’ies.  3  e  r..;ve 


’.'(T-oi  :  '>t- 


The  , 'round  deger.-’ra'ry  co'"'^s  from  the  ; 


tions  of  ph  'es  of  cliains,  or  in  .th 


rirst  spin  on  ej'sh  chain,  wh.c:;  eq  il.' 


the  u  n  d !  s.  t  ^  r  t  :■  d  i  a  t '  .  -e ,  !  t  .  s 


I  n  this  case,  f  :  r  '  r , . 


t  w  '  I  r  p .  n  ;■  i  ; 


e  V  e  n  with 
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crystal  is  cut:  if  we  cut  the  lattice  regularly  aloii;;  surfaces  .:■£  c.i--  u;at 
cells,  each  te trah. eUrc u  on  the  surface  contr ihu'.es  1/?.  "free”  spin,  •.Yf.ile  a 
more  irr?,eular  cut  can  allow  each  tetralieUr'^r!  on  t.;'_‘  su-f  ice  to  cori*:  rib'.: t  e  1 
"free”  spin.  The  ma.  .iir.  i  tucle  of  the  Curie  tail  in  our  ca  1  :u  ’  at :  ■  is  coiuesr  c  n  hr. 
exactin’  to  the  numhor  cl  "ire;"  spins  v;e  produce  fy  th'  h’lindary  nondit 
Therefore,  1  usir.n  a  free  boundary  condition  inst'^ad  ■  f  ;  p-iuclic  .ere,  we 
get  an  additional  paramagnetic  tail  in  t.he  magnetic  susceotilnli  t  y  f  .er  ..11  *,:ie 
cases  previously  cal:u!.ated,  in  bct.h  distorted  '"d  un  1  istorted  lattices. 

5.  Conclusions 

The  study  of  the  tetragonally  transfer.-  'd  spinal  ’-ftice  hr-.-,  le.ad  to  a  r.um.ber 
of  further  insights.  sharp  thermol  -'.amic  t,' ansi  cion  occurs  '.vhen  the  frus¬ 
tration  IS  reduced  in  the  ideal  spinel  lattice;  *he  data  .show  that  samples  of 
both  distorted  and  undistorted  nature  can  apparenr.17  be  made.  The  m. .1  .yn e 1 1  c 
properties  are  relatively  insensitive  to  tlie  Irttice  distortion,  but  rather 
depend  sensitively  .in  the  surface  condition  of  tlie  sinnoC  .jriins,  co'mractin.g 
the  grains  into  a  denser  compo.site  m.i*  rial  .apP‘?'  '  tr  rc-.duce  bit  no't  e:t.:'ol;. 
eliminate  the  free  spins  t  on  the  surface.  A  number  of  intriguing  :.:urs- 
tions  remain  concerning  the  nature  of  the  phase  diagram  a.  a  function  of  th.e 
distortion  or  properties  of  the  composite  in  w'hlch  the  grains  are  sintere.i. 
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ilopefully,  single  crystals  will  lao'como  available',  an 
would  have  a  major  impact  on  ti-e  research.  Small  crystal 
cold  be  uscol  in  all  measuromer.ex,  .except  thei"  al  conductiv 
may  be  too  optimistic  to  expect  larC'C  crystals,  but  such, 
would  allow  measu  r  j.mon  ts  of  anisotropic  properties. 
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phase  diaepam  for  the  distxuted  npir-.  ’  structures  so  ti'.at 
■ansition  temperature  can  be  accuratelv  predicted  f'-'i'  -r'x 
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mAVEL: 

a.  By  tbie  Research  Associate: 
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